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dynamic  programming. 


Roman  Symbols 

AD  Design  tolerance  for  the  quality  characteristics  (tolerance  on  diameter  for 

the  bearing)  (in). 

A.  Dynamic  load  multiplier. 

Xk  Lagrange  multiplier  for  the  &'th  constraint. 

Mf)  Failure  rate  of  a  part  as  a  function  of  its  operation  time, 

p  Lubricant  viscosity  (Reyn). 

p  Oil  density  (0.0311  lbf/in3). 

to  Bearing  speed  (radians/sec). 

0  A  constant  in  lubricant  oil  viscosity-temperature  relation 

v|/(a)  Lund  and  Saible  instability  criterion 

p  Oil  density  (0.0311  lbf/in*). 

a0  Standard    deviation    of  the   quality   characteristic    of  interest    in  the 

manufacturing  process. 
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The  objective  of  this  research  is  to  address  the  need  for  a  comprehensive 
approach  to  designing  quality  products  with  optimum  overall  life-cycle  cost. 
Product  life-cycle  considerations  include  research  and  development,  design 
innovations,  manufacturing,  quality  assurance,  operations,  support,  maintenance  and 
waste  disposal. 

The  emphasis  in  this  research  is  on  developing  an  algorithmic  framework  for 
integrating  the  many  separate  efforts  that  deal  with  each  phase  of  the  product  life- 
cycle,  in  order  to  design,  manufacture,  operate  and  maintain  the  product  with 
minimum  life-cycle  cost.    A  further  objective  is  to  develop  a  framework  for 
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quantitative  determination  of  the  value  added  from  research  and  development,  and 
design  innovations.    This  framework  can  aid  in  decision  making  for  investment  in 
research  and  development  and  acquiring  new  technology  by  performing  cost/benefit 
analysis. 

To  illustrate  the  approach  an  integrated  decision  system  for  the  optimum 
design  of  automotive  journal  bearings  has  been  developed.    The  stages  in  product 
development,  manufacture  and  operation  of  the  system  are  serial  in  nature;  hence,  a 
multilevel  and  dynamic  programming  approach  combined  with  nonlinear 
optimization  techniques  are  used  for  overall  optimization.    Sensitivity  of  the 
optimum  solution  to  problem  parameters  (design  variables  and  design  constraints)  is 
evaluated  in  order  to  perform  cost/benefit  analysis  of  investments  in  R&D  and  new 
concepts. 
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CHAPTER  1 
INTRODUCTION 


1 . 1    Background 

Due  to  the  globally  competitive  nature  of  today's  manufacturing  industries  and 
rising  costs  of  products,  it  is  imperative  to  design  products  which  are  economical  to 
produce  and  operate.  If  the  products  are  designed  only  for  optimal  performance,  their 
manufacturing,  operations  and  support  costs  may  not  be  optimum.  There  is  a  growing 
need  to  design  products  and  systems  in  terms  of  optimum  life-cycle  costs,  rather  than 
just  a  segment  of  its  life-cycle  costs  such  as  production  costs.  Designing  the  products 
and  systems  in  terms  of  optimum  life-cycle  costs  ensures  an  overall  economy  and  not 
the  economy  of  purchase  price  only.  Also,  in  many  industries  it  has  been  recognized 
that  the  greatest  impact  on  total  life-cycle  cost  results  from  the  design  decisions  made 
at  the  early  stages  of  the  design  of  a  product.  For  example,  the  material  specifications 
or  allocation  of  tolerances  at  design  stage  can  have  great  impact  on  the  cost  of 
manufacturing,  quality,  reliability  and  life  of  the  product.  Thus  life-cycle  cost  is  of 
paramount  importance  in  the  decision  making  process  and  it  must  be  considered  as  a 
major  evaluation  criterion  at  the  design  stage.   The  life-cycle  cost  directly  or  indirectly 
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accounts  for  other  parameters  such  as  product  performance,  effectiveness,  size,  weight, 

capacity,  manufacturability,  reliability,  maintainability  and  so  on. 

Before  a  design  engineer  can  accommodate  the  constraints  and  cost  of 
manufacturing,  quality  assurance,  operation  and  support  of  the  product,  information 
must  be  supplied  to  the  designer  in  a  manner  that  can  be  readily  used.  The  design 
engineer,  then,  not  only  can  evaluate  the  performance  of  the  product,  but  also  can 
accommodate  the  constraints  imposed  by  other  functions  (departments)  of  an 
organization.  The  components  of  the  life  cycle  cost  can  be  expressed  as  a  function  of 
various  parameters  and  decision  variables  which  affect  its  life-cycle. 

Often,  the  minimum  life  cycle  cost  may  be  outside  the  design  constraints.  For 
example,  a  designer  may  find  that  a  product  design  is  constrained  by  minimum  surface 
roughness  producible  by  the  existing  processes  and  machines.  In  the  designer's  opinion 
a  better  surface  roughness  might  lower  the  life-cycle  cost  of  the  product  in  terms  of 
decreased  quality  cost  and  increased  useful  life.  In  order  to  improve  the  surface 
roughness  of  the  design,  a  cos^enefit  analysis  should  be  performed  to  analyze  the 
financial  resource  allocation  in  either  buying  new  machines  and  R&D  to  develop  new 
processes  or  analytical  theories.  A  sound  methodology  for  performing  such  a 
cost/benefit  analysis  for  expanding  the  design  space  should  be  available  to  the  designer 
in  order  to  explore  the  potential  gains  in  the  optimum  life-cycle  cost  and  consequently 
to  justify  any  new  investments.  The  development  of  such  a  methodology  is  also  one 
of  the  goals  of  this  research.  The  methodology  integrates  the  scientific  and 
management  aspects  of  the  product  realization. 


1.2   Relevant  Literature 


In  spite  of  the  considerable  interest  and  practical  importance,  the  development 
of  automated  systems  which  integrate  all  aspects  of  the  design,  manufacturing,  quality 
assurance,  and  operation  aspects  along  with  the  cost  management  has  not  received 
appropriate  attention. 

1.2.1    Optimum  Parameter  Design  for  Performance 

Corser  and  Seireg  [Cors85]  developed  an  algorithmic  framework  for  integration 
of  design,  manufacturing  and  inspection  stages  of  the  life  cycle.  An  automated  design 
procedure  was  developed  using  simplified  cost  models  of  manufacturing  and  quality 
control  functions  that  reflect  the  overall  cost  patterns  as  a  function  of  design  parameters. 
A  simple  hydrodynamic  bearing  was  used  to  illustrate  the  approach.  A  dynamic 
programming  method  was  used  to  optimize  the  overall  product  cost. 

Many  other  systems  have  been  developed  to  optimize  the  individual 
departmental  tasks  such  as  design  of  a  mechanical  element.  Seireg  and  Ezzat  [Seir69] 
developed  an  automated  design  system  for  the  selection  of  length,  clearance,  and 
lubricant  viscosity  to  optimize  a  performance  characteristic  of  hydrodynamic  journal 
bearings.  The  objective  was  to  penalize  high  temperatures  and  pressures  in  the 
designed  bearing.  Sutherland  and  Roth  [Suth75]  considered  manufacturing  tolerances 
and  costs  in  the  design  of  mechanisms  for  function  generations  problems. 
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1.2.2   Optimum  Inspection  Schemes 

Many  schemes  have  been  presented  for  optimal  inspection  and  quality  control 
during  production.  After  a  design  is  handed  over  to  the  manufacturing  department, 
these  schemes  deal  with  the  optimal  inspection  strategy  during  production  in  order  to 
conform  to  the  prescribed  design  as  closely  as  possible.  The  objective  here  is  to 
determine  the  inspection  strategy  which  maximizes  the  average  outgoing  quality  and 
minimizes  the  cost  of  inspection.  These  schemes  try  to  achieve  a  balance  between  the 
cost  of  inspection  and  the  cost  of  outgoing  (external)  quality.  But  these  inspection 
schemes  do  not  try  to  influence  the  product  specifications  at  the  design  and 
development  stage. 

A  skip-lot  sampling  plan  initially  calls  for  inspection  of  every  lot  (normal 
inspection)  and  inspection  of  only  a  fraction  of  the  lots  (skip-lot  inspection)  when  lot 
acceptance  demonstrates  the  quality  to  be  good.  This  results  in  reduction  of  average 
number  of  items  inspected  while  maintaining  a  probability  of  acceptance  at  least  as  high 
as  lot-by-lot  inspection  would  yield.  The  important  decision  parameters  here  are  the 
sample  size  inspected,  and  the  acceptance  number.  Other  skip-lot  parameters  are  the 
number  of  good  lots  before  skipping  is  introduced,  and  the  fraction  of  lots  to  be 
inspected  during  skip  lot  sampling.  Hsu  presented  an  optimal  scheme  for  inspection  of 
skip-lot  sampling  plans  for  simple  processes  [Hsu80].  He  also  presented  an  optimal 
scheme  for  skip-lot  inspection  of  multistage  production  processes  [Hsu84].  The 
objective  was  to  minimize  the  average  production,  appraisal,  and  failure  cost  per  unit 
produced. 
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X-bar  (X)  charts  are  used  during  production  for  process  control  where  the 

quality  characteristic  of  interest  is  a  continuous  variable.  In  the  installation  of  an  A' 
control  chart  monitoring  system,  the  user  must  specify  the  four  quantities:  the  sample 
size  n,  the  time  interval  (in  hours)  between  samples  h,  the  average  value  of  quality 
control  characteristic  x ,  and  the  width  of  the  control  limit  k  (quality  control  limits  = 
x±ko).  Many  researchers  have  presented  schemes  for  optimum  economic  design  of  A' 
control  charts.  Montgomery  [Mont82a]  presented  a  scheme  to  choose  the  three 
parameters  {h,n,k)  such  that  expected  loss  per  hour  is  minimized.  Collani  [C0II86] 
presented  a  scheme  to  choose  these  parameters  so  that  expected  loss  per  item  is 
minimized. 

Other  inspection  and  quality  control  optimization  studies  for  different  sampling 
plans  have  been  reported  in  the  literature.  Some  of  the  important  studies  include 
Hassan  and  Knowles  [Hass79].  Gibra  [Gibr78,  Gibr81]  presented  models  for 
economical  optimal  determination  of  the  parameters  of  np-control  charts.  Collani  and 
Sheil  presented  an  approach  for  controlling  the  process  variability  [Coll89].  They  gave 
an  economically  optimal  method  of  controlling  a  process  using  X  chart  and  s  chart  so 
as  to  maximize  the  average  profit  per  item  produced.  The  method  determines  the 
optimal  values  of  sampling  interval  h,  sample  size  n,  and  parameter  k.  The  process  is 
considered  out  of  control  whenever  a  plotted  point  on  the  s  chart  exceeds  the  limit  ka0. 
Montgomery  [Mont2]  presents  an  excellent  overview  and  literature  survey  of  economic 
design  of  control  charts. 
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1.2.3    Statistical  Tolerance  Analysis 

Many  quality  engineering  studies  deal  with  the  techniques  for  statistical 
tolerancing.  Evans  presents  an  excellent  description  of  the  state  of  the  art  in  tolerancing 
[Evan74,  Evan75a,  Evan75b].  The  statistical  tolerancing  is  generally  taken  to  mean  the 
following  two  aspects:  (1)  the  problem  of  ascertaining  the  distribution  of  the  response 
of  a  mechanism,  when  the  component  tolerance  distributions  are  given,  and  (2)  the 
problem  due  to  shifting  and  drifting  of  component  tolerance  distributions.  These 
references  do  not  deal  with  the  problem  of  determining  the  component  tolerances,  given 
the  overall  tolerance  on  the  performance  characteristics  of  the  system. 

An  implementation  aspect  of  tolerancing  is  the  tolerance  allocation,  where  the 
task  is  to  find  the  component  tolerances,  given  the  assembly  tolerance.  Chase  and 
Greenwood  [Chas88]  presented  a  tolerance  allocation  method  which  minimizes  the  cost 
of  production  or  assembly.  This  is  accomplished  by  defining  a  manufacturing  cost 
versus  tolerance  curve  for  each  component  in  the  assembly.  The  optimization  algorithm 
varies  the  tolerance  for  each  component  and  searches  for  the  best  combination  of 
tolerances  that  minimizes  the  manufacturing  cost.  Taguchi's  quality  engineering 
methods  include  optimum  tolerance  design  as  presented  in  the  following  section. 

1 .2.4   Taguchi  Methods  of  Quality  Engineering 

Genichi  Taguchi's  approach  to  quality  involves  engineering  and  statistical 
methods  to  achieve  improvements  in  cost  and  quality  by  optimizing  the  product  design 
and  manufacturing    process.      Taguchi    defines   two   aspects   of  an  overall   quality 
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engineering  system  —  on-line  and  off-line  quality  control  [Tagu78a].     The  quality 

control  activities  at  the  manufacturing  stage,  which  are  conducted  to  keep  the 
manufacturing  process  in  statistical  control  and  to  reduce  the  manufacturing  variations 
(imperfections)  in  the  product,  are  termed  on-line  quality  control  methods.  Quality 
control  methods  such  as  cause  and  effect  diagrams,  process  capability  studies,  process 
quality  control,  and  control  charts  are  referred  to  by  Taguchi  as  on-line  quality  control 
methods.  Off-line  quality  control  methods  are  quality  control  and  activities  conducted 
at  the  product  or  process  design  stage  in  the  product  development  cycle.  The  overall 
aim  of  the  off-line  quality  control  activities  is  to  improve  product  reliability, 
manufacturability,  and  to  reduce  product  development  and  product  operation  costs.  Off- 
line quality  activities  include  design  reviews,  sensitivity  analysis,  prototype  tests, 
accelerated  life  tests,  and  reliability  studies. 

Taguchi  also  formalized  a  three-step  process  for  off-line  quality  engineering 
[Tagu86,  Tagu89]  which  includes  system  design,  parameter  design,  and  tolerance 
design.  System  design  is  a  process  of  applying  scientific  and  engineering  knowledge 
to  produce  a  basic  functional  prototype  design.  The  system  design  defines  the  initial 
settings  of  product  or  process  design  characteristics.  Once  the  system  design  is 
accomplished,  parameter  design  is  undertaken  to  identify  and  establish  the  optimum 
values  of  various  parameters  (nominal  settings)  that  minimize  the  performance  variation 
in  the  system.  The  variations  in  performance  characteristics  from  its  optimal  values  are 
generally  a  function  of  the  design  parameters. 
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After  optimum  parameter  design  for  a  product,  an  important  and  often  critical 

task  faced  by  a  designer   is  tolerance   design  —  determining   the  manufacturing 

tolerances  around  the  nominal  settings  for  various  parameters.    Tolerance  design  is  a 

method  for  determining  tolerances  that  minimize  the  sum  of  product  manufacturing 

cost,  and  cost  of  product  operation  over  its  lifetime.    Tolerances  that  are  too  narrow 

increase   the   manufacturing    cost,   and  tolerances   that  are  too  wide   increase   the 

performance  variation  and,  therefore,  a  product's  lifetime  operation  cost. 

Taguchi's  tolerance  design  methods  are  based  on  minimizing  the  manufacturing 
and  operational  costs.  He  proposed  relationships  between  tolerances,  quality  control 
parameters  and  manufacturing  capability  by  using  quadratic  functions  to  model  cost  as 
a  function  of  tolerances.  In  these  studies,  the  total  cost  includes  inspection  costs  and 
the  cost  of  deviation  of  performance  due  to  variance  of  parameters  from  the  target  or 
design  value. 

It  is  still  a  common  practice  in  the  industry  to  assign  tolerances  by  convention 
rather  than  by  scientific  considerations.  In  the  development  of  the  integrated  decision 
system  in  this  research,  some  of  the  Taguchi  methods  for  modeling  the  quality  control 
function  will  be  applied.  Further  discussion  of  quality  control,  inspection,  and  process 
modeling  will  be  presented  in  Chapter  5. 

1.2.5    Optimizing  the  Manufacturing  Cost 

There  have  been  numerous  studies  on  the  topic  of  optimum  design  for 
manufacture,  optimum  manufacturing  plant  design  and  scheduling  problems.    A  few 
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examples  of  such  studies  of  product  design  for  optimum  manufacturing  cost  are  listed 

here.  Hayes,  Davis  and  Wysk  [Haye81]  presented  an  approach  to  determine  the  number 

of  machines   and  their  operating  rates  for  each  machine   center  in  a  production 

environment.     They  used  a  dynamic  programming  approach  to  solve  the  problem. 

Sutherland  and  Roth  [Suth75]  considered  the  manufacturing  tolerances  and  costs  in  the 

design  of  mechanisms  in  function  generating  problems.     The  textbook  by  Trucks 

[Truc76]  discusses  the  methods  of  economical  production.    Boothroyd  articulated  the 

need  for  developing  and  using  economic  criteria  for  design  for  economical  production 

and  product  assembly.    For  example,  he  proposes  [Boot84]  methods  for  design  of  parts 

so  that  the  assembly  cost  is  minimized.    He  also  presents  [Boot85]  the  models  for 

approximate  cost  estimating  of  turned  parts.    In  his  1990  work  [Boot90],  he  proposes 

models  for  costing  procedures  which  are  intended  to  form  a  basis  of  design  for 

economical  manufacture  (DFM).   Here  he  proposes  the  manufacturing  cost  models  for 

machining  and  injection  molding  processes.    The  idea  is  to  use  these  models  to  choose 

the  appropriate  materials  and  manufacturing  processes  at  the  design  stage. 

1.2.6   Investment  in  Manufacturing,  Quality,  and  New  Technologies 

New  accounting  methods  for  estimating  manufacturing  costs  are  needed  which 
will  enable  us  to  accurately  model  the  cost  of  manufacturing  including  such  intangibles 
as,  flexibility,  quality  and  new  technology.  This  will  also  allow  the  designers  to 
optimize  the  manufacturing  costs  at  the  design  stage  and  make  sound  decisions  on 
investments  in  manufacturing  systems  and  new  technologies.   The  present  management 
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costing  practices  place  too  much  emphasis  on  labor  costs  and  do  not  properly  measure 

intangibles  such  as  production  flexibility  and  product  quality.     Over  the  past  few 

decades  the  cost  of  labor  has  been  steadily  decreasing  as  a  percentage  of  the  total  cost 

of  a  product.     Hence  the  practice  of  measuring  overheads  using  labor  cost  is  not 

justifiable  anymore.  Methods  to  measure  the  value  of  intangibles  such  as  flexibility  and 

quality  have  been  proposed  but  are  not  yet  practiced  widely. 

1.3   Traditional  Product  Realization  Process 

Typically,  when  developing  an  automated  decision  system  for  product  design, 
the  main  objective  is  to  optimize  only  the  performance  characteristic  of  the  product. 
Once  the  optimum  design  parameters  are  found,  the  concerns  of  how  to  manufacture 
the  product  are  addressed.  If  necessary,  the  design  may  involve  one  or  more  iterations 
in  the  design  stage  based  on  any  new  constrains  found  from  the  efforts  to  realize  the 
design  in  manufacturing.  These  constraints  typically  result  from  the  limitations  of  shop 
capacity,  or  capability,  or  possibly  the  lack  of  proper  quality  designed  into  the  product. 
These  iterations  generally  impose  unexpected  additional  costs  and  delays  of  the  product 
delivery.  Later,  as  the  product  is  used  or  operated,  the  customer's  complaints  and 
suggestions  are  fed  back  to  the  design  engineer,  thus  requiring  more  design  iterations 
in  order  to  address  the  operation,  quality,  reliability  and  maintainability  issues. 

The  above  approach  to  product  design  has  been  the  common  method  of  product 
realization  in  most  manufacturing  industries.    Subsequently  there  is  little  direct  input 
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of  manufacturing  concerns,  process  capabilities,  quality  and  maintainability  for  the 
design  engineer  to  incorporate  in  the  design  decisions  such  that  a  universal  optimum  is 
obtained.  Universal  optimum  in  this  context  implies  the  minimum  life-cycle  cost  with 
highest  quality  and  reliability  subject  to  given  operational,  manufacturing  and 
technological  constraints. 

Figure  1 . 1  shows  the  stages  in  product  life  cycle.  The  dashed  lines  indicate  the 
possible  feedback  presented  for  design  iterations.  Such  iterations  are  a  result  of 
unimplementable  or  non-optimum  design  specifications  such  as  shape  or  dimensions 
that  are  impossible  to  manufacture;  tolerances  that  are  unrealistically  tight;  materials 
that  are  costly,  unsuitable  to  task  or  manufacture;  performance  characteristics  that  are 
unattainable  as  designed;  or  quality  problems. 

The  decisions  of  normal  product  development/realization  efforts  at  each  stage 
of  the  product  life  cycle  are  sequential  in  nature.  The  output  of  one  department  (stage) 
is  the  input  to  the  next  department  (stage).  Often  the  decisions  made  in  one  department 
are  constraints  to  the  other  department.  Each  department  has  different  immediate 
concerns  in  realizing  a  product  but  ultimately  has  the  same  goal  —  build  the  best 
product  at  the  lowest  possible  cost.  Each  department  tries  to  minimize  the  cost  of  its 
own  separate  task,  but  one  unit  of  cost  reduction  in  one  department  may  not  necessarily 
lead  to  a  decrease  in  overall  costs.  Thus  individual  departmental  cost  reduction 
decisions  independent  of  other  departments  will  in  general  lead  to  suboptimal  life-cycle 
product  design. 
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1.4   Modeling  the  Product  Life  Cycle 

The  life-cycle  cost  of  a  product  can  be  defined  as  the  cost  incurred  to  society 

from  the  product  conception  stage  to  the  product  disposal  stage.    The  stages  in  the 

product  life  cycle  can  roughly  be  describes  as  follows 

product  planning  and  marketing 

research  and  development 

design  and  prototype 

manufacturing 

quality  assurance 

operation,  support  and  maintenance 

disposal 

In  this  research,  the  overall  life-cycle  cost  is  considered  to  consist  of 

•  manufacturing  cost,  CM 

•  operation  cost,  C0 

•  internal  quality  cost,  CIQ 

•  external  quality  cost,  CEQ 

In  order  to  design  the  products  for  optimum  life  cycle,  costs  in  all  the  stages  of  the  life- 
cycle  as  functions  of  various  design,  manufacturing,  and  quality  control  parameters  can 
be  formulated.  Next,  constrained  optimization  techniques  can  be  used  to  optimize  the 
life-cycle  cost  as  a  function  of  the  parameters  of  design,  manufacturing,  and  quality 
control. 

The  block  diagram  in  Figure  1 .2  represents  the  decision  making  process  for  the 
life  cycle  of  a  product  as  a  set  of  serially  connected  stages.  The  following  discussion 
assumes  that  there  are  /  stages.  The  output  of  one  stage  is  the  input  to  the  next  stage 
of  decision  making  and  can  have  great  influence  on  the  subsequent  stages.    Often  the 
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parameters  of  one  stage  act  as  constraints  to  another  stage,  thus  restricting  the  design 

space.  A  stage  of  the  life  cycle  can  be  represented  by  a  block  diagram  as  in  Figure  1.3. 

The  stage  i  has  a  set  of  input  parameters  (or  state  variables)  S^,  a  set  of  decision 

variables  (or  stage  design  variables)  X:,  a  set  of  transformation   equations  Fip  an 

objective  function  U^  and  a  set  of  output  parameters  St. 

The  stage  model  for  each  stage  is  given  by  its  transformation  equations  that  map 

the  stage  input  parameters  and  stage  decision  variables  to  the  stage  output  parameters. 

The  stage  model  receives  input  parameters,  then  in  conjunction  with  stage  decision 

variables  and  stage  behavior,  it  produces  the  stage  output.   Thus,  the  output  of  stage  i, 

Sj,  is  a  function  of  the  input  parameters  £,.,  and  stage  decision  variables  Xt.    This 

functionality  for  stage  i,  called  transformation  equations,  can  be  expressed  as 

S,  =  F/S,.,,  *,),   for  i  =  1,...,/. 

The  transformation  equations,  objective  functions,  and  constraints  could  be  of  simple 
closed  form,  empirical  charts,  or  they  may  require  numerical  modeling. 

The  objective  function  at  each  stage  is  the  cost  component  associated  with  that 
stage.  For  example,  the  design  stage  objective  function  is  the  cost  of  R&D  and  product 
development  and  the  operation  stage  objective  function  is  the  operation,  support  and 
maintenance  for  the  life  of  the  product.  The  stage  objective  function  can  be  expressed 
as 

U,  =  U{S,h  X). 
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Figure  1.2:  Product  life-cycle  as  serially  connected  stages. 
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Design  Decision  Variables, 

x, 

State  Variables,  SM 

STAGE  i 

Si  =  F,(SM,  Xs) 

State  Variables,  St 

(input) 

(output) 

f 


Objective  Function,   Uj(SM,  X,) 


Figure  1.3:  Stage  i  of  the  product  life-cycle. 


17 
Optionally,  a  stage  decision  could  be  constrained  by  its  given  input  and  required  output. 

These  constraints  for  stage  i  can  be  expressed  as 

g*  =  gAS„X,Sd*0,    k=\,...,K, 

Finally,  the  overall  objective  function  for  all  the  stages  combined  is  the  life-cycle  cost, 
which  can  be  expressed  as 

U  =   U(So,  X„  X2 X,)   =    L  £/,(£,„*,) 

To  integrate  the  four  functions  of  research  and  development,  design, 
manufacturing,  quality  assurance,  operations  and  maintenance  into  a  single  system, 
simplified  models  are  developed  in  the  following  chapters  to  describe  each  stage. 

1.5   Research  Objectives 

The  objective  of  this  research  is  to  address  the  need  for  a  comprehensive 
approach  to  the  overall  life-cycle  cost  optimization  of  products  at  the  design  stage.  The 
focus  of  the  present  work  is  on  developing  an  algorithmic  framework  for  integrating  the 
many  separate  efforts  that  deal  with  each  stage  of  the  product  life  cycle.  A  further 
objective  of  the  research  is  to  develop  a  methodology  for  performing  the  cost/benefit 
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analysis  (sensitivity  analysis)  for  making  decisions  to  invest  in  plants,  equipment,  R&D 

and  new  technologies. 

In   this   research,    algorithmic    and   heuristic    methodologies    for   life-cycle 

parameter  tolerance  design  are  proposed.     The  concept  is  based  on  optimizing  a 

product's   life  cycle  cost  per  hour  of  its  operation.      Life  cycle  cost  consists  of 

manufacturing  cost,  internal  and  external  quality  cost,  and  operating  cost  during  the  life 

of  the  product.    These  three  components  of  the  life  cycle  cost  are  strongly  influenced 

by  not  only  the  parameters  but  also  the  tolerances  for  those  parameters  of  the  product 

design. 

Manufacturing  cost  consists  of  material  and  machining  costs.  Quality  cost  is 
considered  to  consists  of  internal  (in-plant)  cost  of  reworking  and  scrapping,  process 
deviation  diagnosis  cost  and  quality  measurement  cost.  External  quality  cost  (during 
usage  by  consumer)  is  the  incremental  operating  cost  due  to  variations  in  the 
manufactured  quality.  Quality  control  parameters  such  as  control  limits  for  quality 
control  charts  and  frequency  of  quality  measurement  interval  are  obtained  as  a  by- 
product of  optimization  process.  Operating  cost  reflects  the  product  performance  and 
maintenance  parameters  during  the  life  cycle  of  the  product. 

In  order  to  determine  the  investments  in  manufacturing  plant  and  equipment, 
research  and  development,  or  new  technologies,  trends  in  life-cycle  cost  as  a  function 
of  the  design  parameters  are  explored.  Research  and  development  and  new  technology 
can  generate  new  alternatives  at  different  stages  of  the  product  life  cycle  and 
consequently  expand  the  design  space  by  shifting  the  design  constraints.    Such  a  shift 
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in  design  constraints  may  result  in  lower  life-cycle  costs  caused  by  the  benefits  due  to 

improved  performance,  reliability  or  quality  of  the  product.    Reasonable  quantitative 

assumptions  are  made  to  estimate  the  cost  of  R&D  and  new  technologies.    But  in  a 

real-life  situation,  where  reliable  cost  data  are  available,  R&D  cost  can  be  determined 

more  accurately.     After  design  optimization  and  life-cycle  cost  determination,  an 

exploration  of  the  design  space  is  done  in  order  to  determine  the  relative  importance  of 

design  parameters  and  design  constraints. 

Hence  a  further  objective  of  the  research  is  to  develop  a  methodology  for 

sensitivity  analysis  of  optimum  design  with  respect  to  life-cycle  decision  parameters  and 

active  constraints.    This  is  necessary  for  quantitative  determination  of  the  value  added 

from   research   and  development,    design   innovations,    new  technologies    and   new 

investments. 

In  this  research,  journal  bearings  with  hydrodynamic  and  thermo-hydrodynamic 

models  are  taken  as  examples  to  illustrate  the  approach  and  the  developed  algorithms. 

The  journal  bearing  life  cycle  is  setup  as  a  sequential  optimization  problem  as  depicted 

in  Figure  1 .2.  To  illustrate  the  developed  methodology,  parameter  and  tolerance  design 

is  undertaken   for  journal  bearings.     The  clearance   in  a  journal  bearing   strongly 

influences  its  performance  and  cost,  and  a  careful  tolerance  design  for  clearance  is 

critical  for  its  optimum  performance.     The  operating  and  waste  disposal  cost  of  the 

bearing  is  evaluated  by  simulating  the  energy  loss  and  oil  replacement  during  its  life. 

The  research  also  deals  with  the  development  of  an  integrated  knowledge  based 

system  for  the  design  of  journal  bearings.   Here  heuristic  suggested  design  is  validated 
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using  the  algorithmic  procedures.    The  class  of  problems  considered  here  assumes  a 

priori  knowledge  of  the  operating  specifications,  with  the  aim  of  generating  an  optimum 

bearing  configuration. 

1.6   Optimum  Solution  Strategies 

1.6.1    Dynamic  Programming  Strategy 

Though  the  problem  is  serial  in  nature,  it  is  highly  coupled  mathematically,  as 
we  will  see  later.  A  special  method  which  can  be  used  to  decompose  (uncouple)  the 
serial  problem  is  the  use  of  dynamic  programming.  Dynamic  programming  is 
particularly  suitable  for  serial  problems  where  the  decisions  taken  at  one  stage  influence 
only  downstream  stages.  Because  the  optimization  is  carried  out  stage  by  stage, 
dynamic  programming  can  be  considered  a  decomposition  technique.  Following 
Bellman's  principle  of  optimality  [Bell57]  each  decision  made  at  any  stage  takes  into 
account  the  effects  of  that  decision  on  downstream  stages.  The  drawback  of  the 
dynamic  programming  is  that  it  is  very  difficult  to  deal  with  the  global  constraints  in 
the  problem.  However,  it  has  the  advantage  of  not  requiring  the  continuity  of  the 
objective  function.  It  is  the  basic  approach  used  in  this  research  for  achieving  the 
global  optimum. 
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1.6.2  Multi-level  Optimization 

A  general  class  of  problems  amenable  to  decomposition  has  a  hierarchy  of 
variables.  The  decomposition  problem  in  such  cases  can  be  called  multi-level 
optimization  problems.  This  kind  of  a  multi-stage  (serial)  problem  can  also  be 
decomposed  into  a  multi-level  (hierarchical)  problem.  Multi-level  is  the  second 
approach  used  in  this  research  for  achieving  the  global  optimum  and  it  is  discussed  in 
detail  in  Chapter  6. 

1.6.3  Heuristic  Based  Strategy 

The  engineering  design  process  is  generally  quasi-procedural,  wherein  heuristic 
suggested  design  is  validated  using  algorithmic  procedures.  The  class  of  problems 
considered  here  assumes  a  priori  knowledge  of  the  operating  specifications,  with  the 
aim  of  generating  an  optimum  bearing  configuration.  The  task  involves  the  preliminary 
synthesis  of  a  suitable  bearing  configuration,  dictated  by  a  decision  making  capability 
and  the  user  specified  operating  conditions.  The  preliminary  configuration  is  analyzed 
to  arrive  at  the  optimum  solution  from  a  minimal  manufacturing  and  life  cycle  operating 
costs  standpoint.  The  procedure  is  iterative  in  nature,  by  which  alternative 
configurations  are  synthesized  using  heuristics,  and  the  best  configuration  is  identified 
based  on  economic  criteria.  A  novel  feature  of  this  approach  is  that  the  design 
constraint  information  is  utilized  as  a  feedback  to  the  heuristic  design  segment  of  the 
program  to  improve  the  optimum  design. 
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Most  design  decision  systems,  now  available  are  usually  limited  to  problem 

solving  tasks  in  the  deterministic  or  algorithmic  domain.   However,  creative  aspects  of 

design  are  often  ill-structured  and  require  the  judgement  of  the  designer.    Therefore 

attempts  have  been  made  in  the  past  few  years  to  incorporate  qualitative  decision 

making  capabilities  based  on  experience  in  design  procedures  [Dixo83,  Liu87,  Haje89]. 

In  this  research,     an  effort  has  been  made  to  present  a  design  methodology  which 

integrates  knowledge  based  methods  and  algorithmic  analysis/optimization  methods  for 

optimum  life-cycle  design. 

The  following  three  chapters  describe  the  four  components  of  the  life  cycle  cost. 
The  journal  bearing  design  problem  is  formulated  in  the  next  chapter.  Operation  and 
support  stage  is  modeled  in  Chapter  3,  where  a  bearing  simulation  procedure  is 
developed  in  order  to  evaluate  the  bearing  performance  and  calculate  the  cost  of 
operation  and  support.  The  manufacturing  and  quality  assurance  stages  and  associated 
cost  formulations  are  described  in  Chapter  5. 

In  Chapter  6,  the  optimization  strategies  are  described  in  detail  where  the  process 
of  uncoupling  the  problem  and  making  it  into  a  truly  sequential  problem  so  that  it  can 
be  formulated  in  terms  of  dynamic  programming  and  multilevel  optimization  techniques 
is  also  discussed.  This  is  followed  by  the  incorporation  of  heuristic  methods  in  the 
integrated  design  system.  Chapter  7  provides  a  few  examples  of  life-cycle  optimization 
using  the  integrated  design  decision  system.  Chapter  8  presents  a  few  examples  in  cost- 
benefit  analysis  of  alternative  investments  in  order  to  improve  life-cycle  cost.   Finally, 
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Chapter  9  presents  the  summary  of  the  work  and  concludes  with  some  recommendations 

for  further  research. 


CHAPTER  2 
MODELING  THE  DESIGN  STAGE 


The  objective  at  the  design  stage  is  to  specify  the  detailed  system  parameters 
such  that  the  desired  system  performance  characteristics  are  achieved.  After  having 
designed  a  system,  the  design  specifications  can  then  be  passed  onto  manufacturing, 
quality  assurance  and  operations  stages  as  input  to  these  stages.  The  purpose  here  is 
not  to  model  the  design  completely  or  perform  an  exhaustive  analysis  of  all  the  minor 
variables.  Rather,  the  purpose  is  to  select  the  design  parameters  and  performance 
characteristics  that  effect  the  manufacturing,  quality,  and  operation  and  support  costs. 
Another  objective  of  the  design  stage  is  to  investigate  alternative  designs;  applicable 
alternative  theories  for  design  relationships;  and  perform  cost/benefit  analysis  of 
investments  in  research,  development  and  new  technologies.  This  aspect  of  the  design 
stage  is  elaborated  in  Chapter  6.  In  this  chapter,  journal  bearing  relationships  for 
hydrodynamic  and  thermo-hydrodynamic  lubrication,  design  specifications  (stage  inputs) 
and  design  variables  (stage  decisions)  are  described.  The  differences  in  analytical 
approach  between  hydrodynamic  and  thermo-hydrodynamic  models  of  lubrication  are 
discussed  and  algorithms  for  evaluating  the  performance  are  presented. 
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2.1    Journal  Bearing  Design  Considerations 

The  nomenclature  used  to  describe  a  journal  bearing  is  depicted  in  a  schematic 
diagram  in  Figure  2.1.  The  centers  of  the  bearing  and  journal  are  at  Ob  and  Oj 
respectively.  The  radial  clearance  C  is  the  difference  between  the  radii  of  the  bearing 
and  the  journal.  L  is  the  length  of  the  bearing.  Figure  2.2  shows  a  schematic  diagram 
for  the  design  stage,  illustrating  the  input  parameters,  decision  variables  and  output 
variables.  The  main  independent  input  parameters  for  the  journal  bearing  under 
consideration  are 

Bearing  load  (lb.),  W 

•  Journal  rotational  speed  (rps),  N 

•  Journal  diameter  (in.),  D 

The  main  design  decision  variables  to  be  determined  by  he  designer  are 

Bearing  length  to  diameter  ratio,  L/D 

•  Radial  clearance  (in.),  C 
Lubricant  average  viscosity  (Reyn),  p. 

Bearing  material  (strength,  hardness,  machinability,  surface  roughness 
and  wear  characteristics) 

•  Bearing  configuration  (shapes,  holes,  groves) 

•  Important  tolerances  —  specifically  tolerances  on  bearing  and  journal 
diameters  (in.),  AD 

The  constraints  on  journal  bearing  design  are  considered  as  follows: 
minimum  oil  film  thickness, 

h0  >  homm  =  5.0x1 0-5  in,  (2.1) 

maximum  outlet  temperature, 

T0  <  Tmax  =  300  °F,  (2.2) 
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Figure  2.1:  Journal  bearing  geometry  and  nomenclature. 
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Figure  2.2:  Schematic  representation  of  design  stage 
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maximum  pressure, 

^^^.u  =  30xl03psi,  (2.3) 

lower  limit  on  the  viscosity  of  the  oil  used, 

u>umjn  =  0.01xlO-6Reyn,  (2.4) 

upper  limit  on  the  viscosity  of  the  oil  used, 

u  <  umax  =  50.0x1 0"6  Reyn,  (2.5) 

upper  limit  on  length  to  diameter  ratio, 

L/D<(L/D)max  =1.0,  (2.6) 

lower  limit  on  length  to  diameter  ratio, 

L/D  >  {L/D)min  =  0.25,  (2.7) 

a  criterion  ensuring  the  dynamic  stability  of  the  bearing  proposed  by  Lund  and  Saibel 
[Lund67],  which  relates  the  dimensionless  rotor  mass  and  a  function  of  modified 
Somrnerfeld  number. 

"»!  <  *a)  (2.8) 

w 

where 

M  =  rotor  mass  per  bearing  (lb), 

C  =  bearing  radial  clearance  (in), 

co  =  journal  rotational  speed  (rad/sec), 

W0  =  bearing  load  (lbf),  and 

ct  =  modified  Somrnerfeld  number  =  7tS(L/D)2 
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vj/(ct)  is  a  function  of  co,  W0,  L,  D,  C  and  u  as  presented  by  Lund  and  Saible 

[Lund67]  and  is  illustrated  in  Figure  2.3.    In  order  to  compute  \\i,  the  curve  in  Figure 
2.3  was  numerically  modeled  as 


a  = 


3.0a055,       a  <  0.28 
6.88a0094,    0.28  <  a  <  2.9 
7.65,  a  >  2.9 


(2.9) 


2.2   Hydrodvnamic  Theory 


In  the  analysis  of  journal  bearing  using  hydrodynamic  lubrication  theory,  the 
following  physical  quantities  are  assumed  to  be  given. 


Bearing  Speed  N 

Load  W 

Diameter  D 

Clearance  C 

Inlet  oil  temperature  T{ 

Oil  type  (SAE  grade) 

A  function  relating  the  oil  viscosity  u  to  average  temeperature  T  as 

follwos 


u  =  Fn  Exp[60/(r+95)] 


(2.10) 


Here  \i0  is  a  constant  representing  the  base  viscosity  for  an  oil  'o'  —  the  SAE  grade  of 
the  oil;  b0  is  a  constant  for  oil  V;  and  FM  is  a  deterioration  factor  representing  the 
increase  in  viscosity  due  to  contaminants  and  oil  deterioration. 
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Figure  2.3:  Journal  bearing  dynamic  stability  criterion  developed  by  Lund  and 
Saibel. 
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The  design  stage  transformation  relations  are  defined  by  the  equations  relating 

the  bearing  performance  characteristics.      In  this  section,  the  system  of  equations 

describing  the  hydrodynamic  journal  bearings  design  will  be  presented.    The  equations 

which  govern  the  operating  characteristics  are  the  Reynolds  differential  equations  for 

hydrodynamic  journal  bearing  [Reyn86].    For  the  bearing  design,  the  Raimondi  and 

Boyd's  numerical  solution  to  Reynolds  differential  equations  for  a  hydrodynamic 

journal  bearing  [Raim58]  are  among  the  most  widely  known.  Their  numerically  solved 

results  are  given  in  the  form  of  design  charts,  with  Sommerfeld  number  (5)  and  length 

to  diameter  ratio  (L/D)  as  the  main  parameters.    These  numerical  solutions  are  based 

on  the  assumption  that  the  film  formed  in  the  bearing  is  isoviscous  and  is  independent 

of  pressure  and  temperature  variations.   The  film  viscosity  (a)  is  calculated  at  the  mean 
temperature  in  the  bearing.    Based  on  isoviscous  hydrodynamic  considerations,  the 

effective  Sommerfeld  number  S  is  defined  as 


S  =  faLID) 


»aN 


C 


2 


(2.10) 


Here  e  is  the  journal  eccentricity  (distance  between  journal  and  bearing  centers),  e  (e/C) 

is  the  eccentricity  ratio  of  the  bearing.   The  objective  of  the  journal  bearing  analysis  is 

to  find  the  following  physical  quantities  in  the  bearing. 

Average  temperature  in  the  bearing  Tmg 
Average  pressure  in  oil  film  P 
Maximum  pressure  in  oil  film  P^ 
Minimum  oil  film  thickness  h0 
Coefficient  of  friction  / 
Oil  flow  rate  throgh  the  bearing  Q 
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Specifically,  the  curve  fittings  by  Seireg  and  Dandage  [Seir82]  to  Raimondi  and  Boyd's 

design  charts  were  used  in  this  research.    The  curve  fitted  Seireg-Dandage  relations 

describe  various  bearing  parameters  as  a  function  of  length  to  diameter  ratio  (L/D)  and 

Sommerfeld  number  (S).   The  parameters  used  in  this  research  which  these  equations 

describe  are  the  temperature  rise  (AT),  maximum  pressure  (P^),  oil  flow  rate  (Q), 

coefficient  of  friction  (/),  and  minimum  film  thickness  (h0)  in  the  bearing  as  the 

following  nondimensional  quantities 

i.  1/,  JL,  -Q-,  H^l  (2.ii) 

C      CJ      Pmax      RNCL  P 

The  equations  for  the  above  quantities  in  [Seir82]  are  of  the  following  form 


nondimensional  quantity  =  K. 


D 


These  equation  are  defined  with  different  values  of  constants  £,,  K2,  K3,  K4  for  different 
ranges  of  Sommerfeld  number  S,  and  length  to  diameter  ratios  L/D,  for  different 
nondimensional  quantities. 

One  of  the  most  important  assumptions  made  in  the  Raimondi-Boyd  numerical 
analysis  is  that  the  viscosity  of  the  oil  film  is  constant.  As  the  work  is  done  on  the  oil 
during  flow  through  the  bearing,  the  temperature  of  the  oil  rises.  From  the  viscosity, 
temperature  equation  given  above,  the  viscosity  drops  off  with  rising  temperature. 
Since  the  analysis  is  based  on  constant  viscosity,  the  task  is  to  find  the  value  of 
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effective  temperature  and  viscosity  to  be  used  in  the  analysis.    To  do  this,  following 

two  quantities  are  first  defined 

T,  =     Temperature  evaluated  using  viscosity  temperature  relation  given  by 

equation  (2.10). 

T2  =     Temperature   evaluated   using  temperature   rise  A71  given  by  Seireg- 

Dandage  curve  fits  to  Raimondi  and  Boyd's  solution. 

Then  T,  and  T2  are  iteratively  changed  in  order  to  minimize  the  difference  between  the 
two  temperatures  (|T,-T2|).  In  order  to  minimize  the  temperature  difference,  golden 
sections  directional  search  method  is  used.  Now  the  effective  oil  film  temperature  is 
taken  to  be  Tavg  =  (T}+T2)/2.  Based  on  the  r,vg  just  evaluated,  h0,  /  Q,  and  Pmax  are 
evaluated  using  Seireg-Dandage  curve  fits  given  by  relations  of  type  in  equation  (2.11). 

2.3    Thermo-hvdrodvnamic  Theory 

Thermo-hydrodynamic  behavior  model  which  takes  into  account  the  thermal 
effects  in  the  journal  bearing,  and  is  more  accurate  in  representing  the  temperatutres  and 
pressures  in  the  beaing,  was  also  used  for  the  design. 

2.3.1    Difference  From  Hvdrodvnamic  Theory 

In  the  classical  hydrodynamic  theory  presented  by  Reynolds,  an  isoviscous  film 
is  assumed.    This  simplification  is  widely  used  in  bearing  design  since  accounting  for 
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the  effects  of  temperature  variations  along  the  lubricant  film  and  across  its  thickness  are 

avoided.      Many   experimental    observations,    however,    show   that   the    isoviscous 

hydrodynamic  theory,  alone,  does  not  account  for  the  load  carrying  capacity  and  the 

temperature  rise  in  the  fluid  film.  A  thermo-hydrodynamic  phenomenon,  as  explained 

by  Seireg  and  Ezzat  [Seir73]  and  other  researchers,  occurs  in  fluid  film  bearings  which 

influences  the  load  carrying  capacity  of  the  film. 

A  procedure  for  journal  bearing  analysis  based  on  thermo-hydrodynamic  model, 

as  proposed  by  Seireg  and  Dandage  [Seir82],  was  also  used  to  design  the  bearing. 

Seireg  and  Dandage  found  that  only  at  certain  oil  inlet  temperature,  for  any  film 

geometry,  oil  and  speed  does  the  thermo-hydrodynamic  phenomenon  produce  no  effect 

on  the  pressure,  and  the  hydrodynamic  theory  can  be  accurately  applied.  This  condition 

is  assumed  to  represent  thermo-hydrodynamic   equilibrium  in  the  fluid  film.     The 

procedure  for  determining  this  equilibrium  state  (O*)  is  as  follows: 

1 .  For  a  particular  film  geometry,  speed  and  oil,  a  curve  can  be  constructed  for  the 
average  pressure  Pa  versus  average  film  temperature  Ta  based  on  hydrodynamic 
theory  as  illustrated  in  Figure  2.4. 

2.  The  point  representing  the  thermo-hydrodynamic  equilibrium  state  (O*)  is 
empirically  found  to  be  the  point  on  this  curve  where  the  slope  of  the  tangent 


is 


dP  v 

—1  =  £  (2.13) 

dT         V 
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Figure  2.4:  Determination  of  thermo-hydrodynamic  equilibrium  point  O*  from 
pressure  Pa  versus  temperature  Ta  curve  for  hydrodynamic  lubrication  of 
journal  bearing  [Seir73]. 
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where 

V-  volume  of  oil  drawn  into  the  clearance  space  in  cubic  inches  per 

revolution. 
K  =  a  constant  which  is  found  empirically  to  be  a  function  of  (R/C)  as  plotted 

in  Figure  2.5  for  all  the  oils  and  speeds. 

3.  The  temperature  rise  AT*  at  the  O*  condition  can  readily  be  determined  based  on 
the  isoviscous  hydrodynamic  theory.  Consequently,  the  oil  inlet  temperature  T, 
corresponding  to  this  condition  can  be  calculated  from 


T.  =  Ta'-—  (2.14) 


4:    For  the  given  eccentricity,  oil  grade,  and  inlet  temperature  Tt,  the  pressure  speed 
relation  can  be  empirically  expressed  as 


1L  (2.15) 

N' 


Since  the  pressure  distribution  as  predicted  by  the  isoviscous  hydrodynamic 
theory  remains  the  same,  therefore, 
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Figure  2.5:      Determination  of  the  thermo-hydrodynamic   behaviour  constant  K 
from  R/C  ratio  [Seir82]. 


38 


(P    ) 
p* 

max 


_^_  (2.16) 

N' 


This  relationship  is  depicted  in  Figure  2.5  and  compared  to  corresponding 
relationship  predicted  by  isoviscous  considerations  for  the  same  conditions. 

2.3.2   Empirical  Procedure  for  Analysis 

For  analysis  of  thermo-hydrodynamic  behaviour  of  journal  bearings,  Seireg  and 
Dandage  defined  a  modified  Sommerfeld  number  S'  [Seir82],  and  consequently  an 
effective  average  viscosity.  This  accounts  for  the  thermo-hydrodynamic  behavior  of  the 
lubricant  film  and  can  be  directly  used  instead  of  Sommerfeld  number  S  to  evaluate  the 
performance  characteristics  for  any  operating  conditions  using  existing  isoviscous 
analysis  and  data.    They  defined  S*  at  equilibrium  point  O*  as 


S'  =  — 

p: 


R 
C 


(2.17) 


The  Seireg  and  Dandage  [Seir82]  empirical  procedure  for  determining  the 
modified  Sommerfeld  number  S*  is  used.  The  5^  can  be  used  instead  of  the  classical 
Sommerfeld  number  (S)  to  determine  the  bearing  performance  characteristics  using 
available  data  and  methods  based  on  the  isoviscous  hydrodynamic  considerations.  In 
this  procedure  the  inlet  temperature  of  the  oil  (7/,)  in  conjunction  with  temperature  rise 
AT  is  used  to  obtain  the  modified  Sommerfeld  number  by  successive  iterations.  In  this 
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SPEED  N 


Figure  2.6:     Pressure  (P)  and  speed  (N)  relationship  for  a  fixed  geometry  bearing 
[Seir82]. 


40 
thermo-hydrpdynamic   bearing  analysis  algorithm,  the  main  objective  is  to  find  the 

modified  Sommerfeld  number  5*.   The  procedure  is  as  follows. 

(1).    Find  the  bearing  characteristic  constant  K  from  Figure  2.5  using  the  value  of 

RJC.    Then  evaluate  the  parameter  RCL/K  for  the  bearing.    For  computational 

purpose,  the  curve  in  Figure  2.5  was  modeled  as  a  following  relation 


K  = 


l^UxlO'^/C)-4-50109,    if(R/Q>200  (2i8) 

10",    if(R/Q<200 

where  a  =  -7.0631  +7.76371og10(tf/Q-2.25814(log10(*/C))2 


(2).    Calculate  the  value  of  the  parameter 


N_  (R\2 

pi  UJ 


(2.19) 


(3).  Make  an  initial  guess  at  T*  and  Pa\  the  average  temperature  and  pressure  at  the 
O*  condition.  Inlet  temperature  T,  and  average  pressure  Pa  =  W/LD  are  used  as 
initial  estimates. 

(4).    Compute  the  dimensionless  oil  factor 


(2.20) 


(r;+o)2 

corresponding  to  the  current  value  of  Ta\ 
(5).    Compute  the  average  viscosity  corresponding  to  the  current  value  of  T* 


Hi  =  Ae  r:*° 


(6).    Calculate  the  approximation  to  5"  by  using  the  formula 
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(2.21) 


s*  =  n;p;4 


R 
C 


(2.22) 


(7).  If  this  approximation  to  S*  is  sufficiently  close  to  the  previous 
approximation  to  S?,  there  will  no  need  for  further  iteration  and  approximation. 

(8).  Corresponding  to  the  current  value  of  5*  calculate  the  quantity  of  oil  flow 
Q/RNCL  from  Seireg-Dandage  curve  fitted  equations. 

(9).    Estimate  the  new  approximation  to  Pa.  by  using  the  equation 


p:  = 


Q    \(RCV 


RNCL){   K 


(2.23) 


((r;+e)2 

(10).  Corresponding  to  the  current  values  of  S*  and  Pj  find  the  temperature  rise  A7* 

by  using  the  curve  fitted  models  of  type  given  by  equation  (2.11). 
(11).  Revise  the  estimates  for  Ta' 


T*  =  T  +  LaT' 

2 


(12).  Go  to  step  (4). 

(13).  Use  the  current  value  of  S'  for  further  analysis  of  the  bearing. 


(2.24) 
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According  to  Seireg  and  Dandage,  the  bearing  performance  predicted  using  the 

above  procedure  for  thermo-hydrodynamic  model  provides  an  excellent  correlation  with 

the  experimental  data  published  earlier  in  the  literature. 


CHAPTER  3 
MODELING  THE  OPERATION  AND  SUPPORT  STAGE 


The  aim  of  this  chapter  is  to  model  the  physical  and  economic  life  and  the  cost 
of  operations  and  support  during  the  life  of  the  product.  Specifically  the  example  of 
journal  bearing  is  used  to  develop  the  relationships  between  the  operations  and  support 
cost,  and  the  design  and  operating  parameters.  Only  those  costs  which  are  influenced 
by  the  bearing  design  and  manufacturing  parameters  are  considered  in  cost  modeling. 
The  operating  cost  of  the  bearing  is  estimated  to  be  the  cost  of  the  energy  loss  in  the 
bearing  during  its  life.  The  support  cost  is  taken  to  be  the  cost  of  oil  replacements  in 
the  bearing  during  its  life  which  includes  the  disposal  cost.  In  order  to  evaluate  the 
total  operation  and  support  cost,  bearing  performance  characteristics  are  simulated  over 
its  life  as  the  bearing  wears  and  its  performance  deteriorates.  Empirical  wear  models 
dependent  on  operating  conditions,  and  bearing  and  journal  materials  are  used.  The 
bearing  life  is  taken  to  be  minimum  of  its  physical  and  economic  life.  The  physical  life 
is  defined  as  the  point  in  time  when  the  bearing  fails  or  one  or  more  design  constraint 
are  violated  and  cannot  be  brought  inside  the  feasible  design  domain  even  after 
replacing  the  oil.  The  economic  life  is  defined  as  the  point  in  time  when  the  total-cost- 
rate  curve  has  reached  its  minimum  value  and  starts  to  rise  with  time.    A  schematic 


43 


44 
diagram  in  terms  of  input,  output  and  decision  variables  of  the  operation  and  support 

stage  is  shown  in  Figure  3.1. 

3.1    Bearing  Life  Simulation 

In  order  to  evaluate  the  bearing  life,  and  its  operation  cost,  a  simulation  of 
bearing  performance  is  done  as  the  bearing  wears  over  time.  The  underlying  concept 
is  that  with  time  the  bearing  wears  and  the  quality  of  oil  deteriorates  because  of 
oxidation  and  presence  of  contaminants.  The  replacement  of  deteriorated  oil  with  new 
oil  brings  the  performance  characteristics  of  the  bearing  closer  to  their  designed 
optimum  values. 

3.1.1    Wear  Modeling 

The  bearing  clearance  increases  due  to  wear,  hence  the  performance  deteriorates 
from  the  optimum.  It  is  understood  that  in  reality  the  wear  along  the  bearing  surface 
(specially  along  the  tangential  direction)  does  not  occur  uniformly.  Wear  is  more 
pronounced  at  the  bottom  compared  to  the  top  half  of  the  bearing,  leading  to  an 
elliptical  shape  over  time.  But  because  of  the  lack  of  a  theory  for  modelling  elliptical 
bearings,  it  is  assumed  in  this  research  that  the  wear  occurs  uniformly  along  the  bearing 
surface.  Hence  with  time,  the  shape  of  both  the  bearing  and  the  journal  remain  exactly 
circular. 


O    l> 

~     60 
U.  13 


45 


M 

.5 

ma 


u 


•o.S 

So 

I  E 


eo 


Q. 
U 


< 

H 

< 

D£ 

H 

O 

c/3 

t/3 

W 
U 

5s 

Z 

z 

o 

w 

p 

5 

2 

5 

w 

^ 

a. 

O 

t    t    t 


c 

>. 

0 

«— 1 

■*-• 

ro 

§ 

g-      2 

60 

nfig 
ype 

t  of 

to 

iatio 

tior 
due 
var 

II 

c  o 

CO 

«-*  ti 

fee— 

.2    (/)    60 

8  3 

■-  "O 

>   o   c 

earing  des 
aterial,  an 

andard  de 
laracteristi 
anufacturi 

action  def 
ipped  to  c 

CQ  E 

(35  "5  E 

£■5 

u 

eq 
8 

6 

5 

3 
Efl 


5 


i 

5 


60 
03 


u 

■a 


I 


46 
To  predict  the  bearing  wear  rate,  we  propose  a  deterministic  model  which 

combines  the  empirical  relationships  for  bearing  wear  rates  with  sound  quantitative 

assumptions  for  statically  loaded  journal  bearing.  Many  experimental  studies  have  been 

conducted  in  the  past  to  empirically  model  the  journal  bearing  wear.  One  of  the  model 

proposed  by  Watanabe  et  al.  [Wata85]  expresses  the  bearing  wear  rate  as  a  function  of 

contaminant  concentration  in  oil;  contaminant  particle  size;  hardness  of  contaminant, 

bearing  and  shaft  material;  and  relative  motion  between  shaft  and  bearing.   Watanabe' s 

model  for  rate  of  bearing  wear  volume  per  revolution  can  be  expressed  as  follows, 


dW 


dn 


^    <*    [*«#,//„//)  +  KWJH^H,)}  r^  r}2  LD 


(3.1) 


Here  Whb  and  Whs  are  functions  of  bearing,  journal  and  shaft  material  hardness,  and  are 
empirically  determined  to  be 


^JHc,Hb,H)  = 


H 

c 

r 

H 

c 

H". 

H 

s 

H 

c 

H 

c 

H" 

H 

s 

0.25 


15 


H 

-1  >  1 
H 


(3.2) 


H 

c 

~H 


<  1 


Notations  in  the  above  three  expressions  are  as  follows 

dWJdn     =     combined  wear  rate  of  bearing  and  shaft  by  volume  as  modeled 
by  Watanabe  (inVrevolution) 
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(3.3) 


Kb,Ks        -     wear  coefficients  for  bearing  and  shaft  material  respectively 

L  =     bearing  length,  (in) 

D  =     bearing  diameter,  (in) 

r^  =     ratio  of  contaminants  volume  to  total  oil  volume  (including  the  oil 

in  sump  and  pump), 
Hc  ,Hb  ,H=     hardness  of  contaminant  particles,  bearing  and  shaft  material 

respectively, 
TV  =     bearing  speed,  (rps) 

rf  =     ratio    of   minimum    film    thickness    to    average    diameter    of 

contaminants  particles 


rf  = 


d        h 

p  o.min 


(3.4) 


Note  that  in  the  above  equation,  it  has  been  assumed  that  the  average  contaminant 
particle  size  (dp)  in  the  oil  is  approximately  equal  to  the  lower  limit  on  minimum  oil 
film  thickness  (homn)  in  the  bearing. 
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Watanabe  et  al.  did  not  take  into  account  the  effects  of  oil  temperature  change, 

pressure,  and  coefficient  of  friction  [Wata85].    In  order  to  take  the  above  effects  into 

account,  the  following  model  for  the  rate  of  wear  volume  is  proposed 


—±    oc    (fP       TY2  (3-5) 

dn  v     max    °J 


where 

/  =  coefficient  of  friction  for  the  oil  film, 

P^  =  maximum  oil  pressure  in  the  film  (psi), 

T0  =  bearing  outlet  oil  temperature  (°F). 

Hence  the  total  wear  volume  AWV,  resulting  in  An  bearing  revolutions  after  taking  into 
account  the  effect  of  material  properties  in  equation  3.1  and  bearing  performance 
characteristics  in  equation  3.5,  can  be  expressed  as  follows 


dW ,    dW?  „_ 

AWv  =  vl  1  An  (3.6) 

dn       dn 


Here  An,  the  number  of  bearing  revolutions  in  At  hours,  is 


An  =  3600  NAt  (3.7) 
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Finally  the  average  increase  in  bearing  clearance  AC  (in)  due  to  wear  in  time  At 

(seconds)  can  be  derived  as 


AC  =  ^JL  ^i   2600NAt  (3-8) 

dn       dn         nDL 


3.1.2   Modeling  Oil  Deterioration 

Viscosity  of  oil  is  assumed  to  change  with  time  because  of  deterioration  of  oil 
quality  and  generation  of  metallic  particles  due  to  wear  and  the  presence  of  contaminant 
in  the  oil.  The  change  in  viscosity  of  oil  is  reflected  by  addition  of  a  coefficient  FM  in 
viscosity-temperature  relationship  of  the  oil  given  by  equation  2.10  and  is  repeated  here 

u  =  FM  u0  Exp[ty(7-+95)]  (3.9) 

where 

\x0  =  A  constant  representing  the  base  viscosity  for  an  oil  V, 

b0  =  A  constant  for  oil  V, 

F^=         A  'deterioration  factor'  representing  the  change  in  oil  viscosity  due  to 
contaminants  and  deterioration. 

For  a  newly  replaced  oil,  F^  is  unity  and  it  varies  with  the  number  of  bearing 
revolutions,  modeled  as  follows 

F„  -  F^  (1  +  n/n0)m  (3.10) 
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where 


w 


n 


nn 


1, 

number  of  revolution  in  the  bearing  since  last  oil  replacement, 
a  constant,  assumed  to  be  2xl06  revolutions. 


After  combining  the  empirical  model  suggested  by  Watanabe  et  al.  with  the 
effects  of  friction,  temperature  and  pressure,  we  get  the  increase  in  clearance  AC  in  An 
number  of  revolutions  (or,  in  time  interval  At=An/3600N  hours)  as  follows 


WF  =  [KJ^Hc,HbJ1)  +  Kf2{Hc,Hb,H)]  r 


1.3 


h 


(3.11a) 


AC  =  WF  (JP       T)  360(WAf 

v     max       o' 


(3.11b) 


After  many  simulations,  the  following  expression  was  developed  for  contaminats  to  oil 


ratio 


iiLD  ,„   „  . 
r     =  r      +  (C-C  ) 

CO  CO,l  xr         v  /KK 


t 
l+_^ 
t 


(3.11c) 


where, 


rcol  =        initial  ratio  of  contaminats  to  oil  (immediately  after  the  old  oil  is 

replaced), 
L,D  =       bearing  lenght  and  diameter, 
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C  =  current  bearing  clearance, 

Cno  =        bearing  clearance  at  last  oil  replacement  (new  oil), 
V0  -         total  oil  volume  including  the  oil  in  sump  and  pump, 
tno  =         hours  of  operation  since  last  oil  replacement, 
t0  =  a  constant  (assumed  to  be  50.0). 

3.1.3    Modeling  Bearing  Instability 

As  the  bearing  wears  and  the  clearance  increases,  the  bearing  may  tend  towards 
dynamic  instability.  As  the  dynamic  instability  is  approached  in  the  bearing  as  modeled 
by  equation  2.8,  the  dimensionless  rotor  mass  (MC(o2/W)  approaches  \]i(a).  The 
approach  towards  dynamic  instability  is  not  a  sudden  phenomenon,  but  rather  an 
incremental  phenomenon.  The  approach  towards  instabilities  causes  fluctuating  loads 
on  the  bearing,  resulting  in  increased  effective  load.  In  order  to  account  for  this  fact, 
a  factor  called  dynamic  load  multiplier  (X),  which  is  modeled  as  follows 


x. L 


1  - 


(3.12) 


\|<a) 


is  multiplied  to  the  specified  static  load  W0.  The  dynamic  load  multiplier  increases  the 
effective  dynamic  load  (W)  on  the  bearing  as  follows 

W=XW0.  (3.13) 
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The  dynamic  load  on  the  bearing  tends  towards  infinity  as  the  dimensionless  rotor  mass 

(MC(02/Wo)  approaches  the  Lund  and  Saibel  instability  criterion  \\i(a)  [Lund67].   In  all 

the  calculations  for  evaluating  bearing  performance,  the  effective  bearing  load  (W)  is 

used  instead  of  the  specified  static  load  W0.    Only  in  evaluating  the  bearing  stability 

criterion,  used  for  determination  of  bearing  life,  is  the  specified  static  load  W0  used. 

3.1.4   Operation  Cost  Model 

Operation  cost  C0(t)  (a  function  of  operation  time)  of  the  bearing  consists  of  the 
cost  of  energy  loss  due  to  friction  Coe(7)  and  the  cost  of  lubricant  replacement  Cor(/) 
over  the  life  of  the  bearing,  and  is  expressed  as 

Co(0   =   CJi)  +  Cor(0   =   rMt)  +  r0V0 
or 

Co(0   =    l(rebJQpcsAT+r0V0)At  (3.14) 

where 

t  =  bearing  operation  time, 

At  =  operation  time  step  (hrs)  used  during  simulation, 

re  =  cost  rate  of  energy  (0.01  $/kWhr), 

H{f)  =  energy  loss  in  bearing  due  to  friction  (kWhr), 

r0  =  cost  rate  for  oil  (0.04  $/in3), 

V0  =  total  volume  of  the  oil  to  be  replaced  (in3), 
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b  =     energy  conversion  factor,  (1.054  Btu/kWhr), 

J  =     mechanical  equivalent  of  heat  (9336  lbf-in/Btu), 

p  -     oil  density  (0.0311  lbf/in3), 

cs  =     heat  capacity  of  oil  (0.42  Btu/lbf-°F), 

Q  =     oil  flow  rate  (inVsec), 

AT  =     temperature  difference  between  oil  inlet  and  outlet  (°F), 

Total  operation  cost  of  the  bearing  is  evaluated  by  simulating  the  bearing  performance 
over  its  life  and  using  the  summation  in  equation  3.14  over  time. 

3.2   Simulation  Procedure 

A  block  diagram  explaining  the  bearing  simulation  is  given  in  Figure  3.2.  The 
bearing  simulation  is  performed  until  the  end  of  its  useful  physical  or  economic  life, 
whichever  occurs  first.  During  bearing  simulation,  as  the  clearance  increases,  one  or 
more  design  constraints  are  violated,  causing  oil  to  be  replaced.  The  bearing  is  assumed 
to  have  reached  the  end  of  its  physical  life  when  replacing  the  bearing  oil  does  not 
bring  the  bearing  back  into  its  feasible  (valid)  design  domain.  The  bearing  economic 
life  is  defined  as  the  time  of  operation  when  the  cost-rate  (life  cycle  cost  per  hour  of 
operation)  has  reached  its  minimum.  The  bearing  is  simulated  using  the  following 
procedure. 
(1).    The  following  quantities  are  assumed  to  be  given  (input  to  subroutine): 


Figure  3.2: 


(     Start      ) 

y 

Choose  Ci,  L,  n, 

I 


Replace  new  oil.  Reset  n  and 
contaminant  concentration. 


T 


Evaluate  incremental  increase  in 
clearance  AC  in  time  At. 
C  =  C  +AC 


T 


Evaluate  deterioration  in  oil  viscosity  /t,  effective 
dynamic  load  W 


Evaluate  bearing  parameters: 

maximum  pressure,  temperature  rise, 
coefficient  of  friction,  minimum  film  thickness, 
effective  oil  viscosity,  contaminant  concentration. 

Evaluate  the  Constraint  Violations. 
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No 
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Is  the  oil  new  ? 
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Economic  life  of 
bearing  over. 


Evaluate:  Total  Life-Cycle-Cost,  Life,  and 
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f   Stop  bearing  simulation    J 

Bearing  performance  simulation  over  its  life. 
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diameter  D;    L/D  ratio;    initial  clearance  Q    viscosity  evaluation 

coefficient  GM. 
(2).    Initialize  the  following  variables: 

simulation  time  f=10;  life=0;  augmented  Lagrange  function,  ALF=0; 

dynamic  load  W=W0;  cost-rate  Cr=0;  energy  loss  cost  Ce=0;  oil  cost 

Co=0;  life-cycle  cost  LCC=Q; 
(3).    Calculate  the  viscosity  u  at  inlet  temperature  T(  using  equation  3.9. 
(4).    Initialize  the  following  parameters: 

time  to  replace  oil  to=0;   number  of  oil  replacements  n0  =0;   number 

of  revolutions   «=0;     constraint   number  violated   v=0;     viscosity 

evaluation  coefficient  Fl=zGll; 
(5).    Using  the  hydrodynamic  or  thermo-hydrodynamic  analysis  procedure  discussed 
in  Chapter  2,  evaluate  the  bearing  performance.    The  following  parameters  are: 

oil  type  V;    bearing  speed  N;    specified  static  load  W0;    bearing 

diameter  D;  bearing  length  L;    oil  inlet  temperature  Tt;     viscosity 

coefficient  F„. 
The  golden  section  search  procedure  is  used  to  evaluate  effective  oil  temperature 
Tmf,  and  oil  viscosity  u.    Using  the  effective  temperature  Tmg  and  viscosity  u, 
evaluate  the  Sommerfeld  number  S  and  then  using  Seireg-Dandage  equations, 
evaluate: 

coefficient  of  friction  /;    oil  flow  rate  Q;    maximum  pressure  Pmax; 

and  outlet  oil  temperature  T0. 
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(6).    Evaluate  dynamic  load  multiplier  X  using  equation  3.12  and  dynamic  load  W 

using  the  equation  3.13. 
(7).    Evaluate  whether  the  following  bearing  performance   constraints   are  being 


satisfied: 


CI: 

h  >  h 

*o           o.min 

C2: 

P         <  P 

*  max         *  max.u 

C3: 

T  <  T 

•*  o         *  max 

C4: 

H  >  Mmin 

C5: 

H<Hm« 

C6: 

r     <  r 

'  co         '  co, max 

C7: 

(L/D)  >  (Z/D), 

C8: 

(I/D)  <  (I/D), 

C9:       C  >  Cmin+AD 


CIO:      C<Cmax-AD 


CI  1:     Lund  and  Saibel  stability  criterion  of  equation  2.8  and  Figure  2.3. 
Note  that  here,  static  load  W0  is  used  to  evaluate  the  dimensionless 
rotor  mass,  rather  than  dynamic  load  W. 
CI 2:      /|ife  >  /|ifemjn 
C13:      AD  >  AD>min 
If  any  constraint  is  being  violated,  initialize  the  variable  V{  with  that  constraint 
number  and  append  it  to  a  string  Fistr  which  stores  the  history  of  constraint 
violations  during  the  life  simulation.    The  constraint  for  the  life  of  the  bearing 
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(CI 2)  is  evaluated  only  at  the  end  of  life,  when  the  constraint  values  are  also 

calculated   and  stored  in  the  constraint  vector  (G)  after  multiplying   by  a 

constraint  satisfaction  factor  (CSF)  as  described  in  the  overall  optimization 

procedures  given  in  Chapter  6. 

(8).    Calculate  the  cost  of  oil  replacement  C0()  and  cost  of  energy  loss  Coe. 


CHAPTER  4 
MANUFACTURING  AND  QUALITY  ASSURANCE  STAGE 


In  this  chapter,  cost  models  for  manufacturing  and  quality  are  developed. 
Manufacturing  cost  consists  of  material  and  machining  costs.  Quality  cost  is  considered 
to  consist  of  internal  (in-plant)  cost  of  reworking  and  scrapping,  process  deviation 
diagnosis  cost  and  quality  measurement  cost.  External  quality  cost  (during  usage  by 
consumer)  is  taken  to  be  the  incremental  operating  cost  due  to  variations  in  the 
manufactured  quality.  Optimum  values  of  quality  control  parameters  such  as  control 
limits  for  quality  control  charts  and  frequency  of  quality  measurement  interval  are  also 
obtained. 

4.1    Manufacturing  Cost  Model 

After  the  design  stage,  the  manufacturing  engineer  develops  a  process  plan 
which  recommends  the  machine  to  be  used  and  the  operating  parameter  to  use  for  each 
machine,  such  as  speed,  feed,  depth  of  cut,  etc.  The  machine  parameters  chosen  and 
tooling  used  are  usually  based  on  an  optimization  algorithm  to  either  utilize  the 
maximum  production  rate  of  the  machine  or  minimize  the  cost  to  process  the  part.  The 
decisions  of  machine  selections  are  based  on  the  ability  of  the  machine  to  create  the 
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required  geometry,  capacity  of  the  machine  (power,  speed  etc.),  and  the  capability  to 

meet   specifications,    in   conjunction    with   the   results    of  the   machine    parameter 

optimization.  Manufacturing  cost  is  heavily  dependent  on  the  design  specifications,  and 

therefore  on  the  machines  selected  and  the  parameters  used  by  the  process  planner.   A 

schematic  diagram  in  Figure  4. 1  shows  the  modeling  of  manufacturing  stage  in  terms 

of  its  input,  output  and  decision  variables. 

For  the  journal  bearing  example,  the  main  concern  for  the  designer  will  be 
specifying  tolerance  for  diameter,  bulk  and  surface  material  properties  of  the  journal 
bearing  such  as  strength,  hardness,  and  surface  roughness.  Final  processing  steps  such 
as  fine  turning,  grinding,  lapping  and  surface  treatment  such  as  coating  will  depend  on 
designer's  specifications.  The  following  factors  must  be  considered  by  the  process 
engineer  for  selecting  the  process  and  machines: 

Capacity  —  It  is  the  size  of  parts  the  machine  can  handle  such  as  diameter, 
length,  power  and  speed  of  the  machine. 

Capability  —  It  is  the  capability  of  the  machines  and  processes  to  attain 
specifications  such  as  geometry,  dimensional  tolerances,  and  surface  roughness.  Such 
capability  for  a  quality  characteristic  of  interest  x  can  be  expressed  by  a  parameter  such 
as  process  capability  index 

USL  -LSL 

C     =  f -  (3-13) 

p-x  6a 

X 

where 

USLX  =  upper  limit  of  design  specification  for  x, 
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LSLX  -  lower  limit  of  design  specification  for  x, 

ax    =  standard  deviation  of  the  process  for  x. 

Availability  —  It  is  a  factor  varying  between  zero  and  one,  which  shows 
whether  the  machine  is  available  for  use  when  requested,  or  difficult  to  obtain.  This 
factor  can  be  used  to  apply  a  cost  penalty  for  choosing  to  use  a  relatively  unavailable 
machine. 

Cost  —  The  cost  to  process  an  order  on  a  given  machine  is  a  function  of  the  lot 
size,  setup  time,  fixture  and  tooling  requirements,  material  handling  time,  processing 
time,  and  the  base  cost  rate  (cost  per  unit  time  for  the  machine).  Also  the  cost  depends 
on  the  design  specifications  such  as  tolerances,  surface  finish,  and  coatings  etc. 

Very  little  data  is  available  in  the  literature  for  the  cost  of  processing  as  a 
function  of  all  the  variables  mentioned  in  this  section.  This  kind  of  cost  information 
may  be  internally  generated  within  the  manufacturing  department  from  the  machine 
history  databases.  The  cost  factors  thus  derived  can  be  very  useful  to  the  design  and 
process  engineers  in  optimum  design  and  machine/process  selection.  A  few 
manufacturing  handbooks  which  present  the  general  trends  in  manufacturing  cost  as  a 
function  of  tolerance  and  surface  roughness  [ASLE,  Neal73,  Truck76,  Tool83]  were 
consulted.  Some  of  these  relationships  are  shown  in  Figure  4.2  and  Figure  4.3. 
Generally  the  following  functions,  which  relate  the  manufacturing  cost  to  tolerance  or 
surface  roughness,  can  be  found  in  these  handbooks 
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Figure  4.2:      Manufacturing  cost  as  a  function  of  tolerance  and  surface  roughness. 
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Figure  4.3:     Manufacturing  cost  as  a  function  of  tolerances  and  surface  roughness 
found  in  handbooks. 
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Cm  oc  A  +  B/y  (4.2a) 

Cm  oc  A  +  B/y2  (4.2b) 

Cm  oc  A  +  Be^  (4.2c) 

where  y  may  be  tolerance  or  surface  roughness.  Seireg  and  Corser  [Cors85]  modeled 
the  base  manufacturing  cost  for  a  journal  bearing  as  follows 

Cm  =  (AADs  +  BAD  +  Cs  +  D)(39I  +  12)  (4.3) 

where  A,  B,  C,  D  are  constants;  AD  is  tolerance  on  bearing  diameter;  s  is  surface 
roughness;  and  L  is  the  bearing  length.  In  this  thesis,  manufacturing  cost  of  the  bearing 
depends  on  the  bearing  configuration,  material,  length  (L),  diameter  (D),  and  tolerance 
on  diameter  (AD)  and  is  modeled  as, 

Cm  =  aMCm.base  (4.4) 


and 


where 


Cm,base  =  K  +  k,LD2  (4.5) 

a  =  bearing  configuration  cost  factor,  with  the  simple  journal  bearing 

having  a  value  1, 

P  =  material  cost  factor,  with  the  lead  based  babbitt  having  a  value  1 , 

y  =  machining  cost  factor  based  on  the  material  hardness, 

(|>  =  tolerance  cost  factor  =  (&0  +  k/AD°), 

k0-k3     =  cost  coefficients:  0.2,  0.5,  25,  1  respectively. 
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a  =  1/3. 

For  the  purpose  of  this  research,  we  assume  that  bearing  has  a  simple  configuration,  and 

is  made  of  lead  based  babbitt,  hence  a,  P,  and  y  are  taken  as  1. 

4.2   Quality  Cost  Model 

Quality  costs  are  the  costs  incurred  because  poor  quality  may  exist  or  it  does 
exist.  Taguchi  defines  the  quality  costs  [Tagu78a]  as  the  "loss  imparted  to  the  society 
by  the  product  from  the  time  the  product  is  shipped."  Generally  quality  costs  are 
categorized  as  follows. 

Prevention  costs:  These  are  costs  incurred  to  prevent  poor  quality  from  being 
produced.  Example  of  prevention  activities  include  analysis  and  planning  for  quality, 
reliability,  quality  training,  and  development  of  process  controls.  The  fixed  part  of 
these  costs  are  not  needed,  since  the  fixed  cost  do  not  change  in  relation  to  the  product 
parameters  to  be  designed. 

Appraisal  costs:  These  are  costs  that  result  from  the  activities  undertaken  to 
prevent  poor  quality  products  from  being  further  processed  beyond  the  point  at  which 
they  become  non-conforming  or  from  being  delivered  to  customers.  Example  of 
appraisal  activities  include  the  inspection  and  testing  of  raw  materials  and  products 
being  produced. 

Internal  failure  costs:  These  are  costs  associated  with  materials  and  products  that 
fail  to  meet  quality  standards  and  are  identified  before  the  product  or  service  is 
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delivered  to  the  customer.     Examples  of  this  include  the  scrap  and  spoilage,  repair, 

rework,  reinspection  and  downtime  because  of  the  quality  problems. 

External  failure  costs:  These  are  the  costs  incurred  because  poor  quality 
products  are  delivered  to  the  customers.  Examples  include  the  costs  of  handling 
complaints,  warranty  replacements,  product  liability,  and  loss  in  performance  of  the 
product  over  its  life. 

In  the  following  sections,  the  first  three  cost  components  will  be  referred  to  as 
internal  quality  costs  and  the  last  component  is  referred  to  as  external  quality  costs. 
A  schematic  diagram  in  Figure  4.4  illustrates  the  quality  control  stage  in  terms  of  its 
input,  output,  and  decision  variables. 

4.3   Internal  Quality  Cost  Model 

The  cost  of  internal  quality  is  measured  using  the  concept  of  Taguchi's  on-line 
feedback  quality  control,  wherein  measurements  of  product  quality  characteristics  are 
obtained  and  analyzed,  and  the  results  fed  back  to  upstream  processes  for  their 
adjustments.  Such  measurements  and  adjustments  of  the  process  can  be  made  either 
automatically  or  by  a  human  operator.  During  the  production,  it  may  not  be  possible 
to  measure  every  product,  hence  an  attempt  is  made  to  measure  and  inspect  a  fraction 
of  the  products  being  produced.  The  major  assumption  made  in  deriving  the  quality 
costs  in  the  following  sections  is  that  only  discrete  products  are  under  production  and 
the  quality  characteristic  of  importance  is  a  continuous  variable  (e.g.  diameter  or  length) 
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as  opposed  to  an  attribute  (e.g.  go/no-go,  pass/fail).    For  the  journal  bearing  problem, 

the  quality  characteristic  of  interest  is  the  diameter  of  the  bearing  and  the  shaft. 
Prevention  and  appraisal  cost  have  been  assumed  to  consist  of  the  following  two 
components: 

1 .  Cost  of  diagnosing  (inspecting)  the  quality  characteristic  per  unit  of  the  product 
produced, 

CIQ1  =  B/n  (4.6) 

where 
B      =         cost  per  measurement  of  the  quality  characteristic, 
n      =  measurement  interval,  i.e.  measurements  made  at  every  n  product. 

2.  Cost  of  process  adjustment  per  unit  of  product  produced, 

C1Q2  =  A/u  (4.7) 

where 
A      =  cost  per  process  adjustment, 

u      =         predicted    average    number    of   products    between    successive 
adjustment  in  the  process. 

Internal  failure  costs  are  actually  the  internal  (in-plant)  losses  caused  by 
deviations  in  a  products  quality  characteristic  from  its  target  value.  Taguchi  proposed 
the  concept  of  loss  function  for  the  quality  characteristic  of  interest,  which  is  the  "loss 
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a  product  causes  to  society  (in-plant  loss  and  loss  to  customers)  due  to  deviation  in 

product  quality  characteristic   from  its  nominal  specified  value."     In  particular,  a 

quadratic  loss  function  for  internal  failure  costs  is  proposed,  given  by 

ClQ3  =  ka02  (4.8) 

where 

a02       =  standard  deviation  of  the  quality  characteristic  of  interest  for  the 

manufacturing  process, 
k  =  a  constant. 

If  the  design  tolerance  of  the  quality  characteristic  (diameter)  is  AD  and  a  deviation  in 
diameter  D  outside  its  tolerance  limits  D0±AD  causes  an  internal  failure  cost  (due  to 
rejection)  A,  then, 

A  =  kAD2  (4.9) 

or 

k  =  A/AD2.  (4.10) 

Hence,  internal  failure  costs  due  to  a  standard  deviation  aD  in  quality  characteristic  can 
be  given  by 


c,3  =  4-  ct»  <411) 

A/, 
where 
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A   =     in-plant  cost  of  rework  or  scrapping  a  unit  that  falls  outside  the  design 

tolerance  interval  D0±AD,  ($) 

DB  -    target  or  nominal  value  of  bearing  diameter,  (in) 

AD  =    design  tolerance  for  diameter. 

Let  the  control  limits  for  the  production  process  adjustment  be  set  at  D0±b.   When  the 
production  process  is  found  to  be  in  statistical  control,  the  variance  of  D  is  given  by 

o-2  =  62/9  (4.12a) 

if  the  manufacturing  process  follows  a  normal  distribution,  and 

ct2  =  62/3  (4.12b) 

if  the  manufacturing  process  follows  a  uniform  distribution.     When  the  production 
process  is  out  of  statistical  control,  the  variance  is  given  by, 


a2  = 


(w+1) 

- i+T 


^  (4.13) 


where 

n  =       process  measurement  interval  in  units  of  products, 
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x  =       lag  in  the  process,  i.e.  number  of  units  produced  from  the  time  the  an 

out-of-control  process  is  detected  to  the  time  the  production  process  is 

stopped  for  recovery, 

u  =      predicted  average  number  of  products  between  successive  adjustments  in 

the  process. 

Usually  u  is  not  known  in  advance,  hence,  in  order  to  evaluate  the  current  value  of  u, 
a  few  trials  are  needed.  If  at  a  time  when  the  quality  control  limits  are  set  at  ±80,  the 
value  of  u  is  u0,  then  u  at  quality  control  limits  5  is  modeled  as, 


u  =  -J152  (4.14) 

82„ 


The  above  expression  is  derived  from  the  assumption  that  the  average  time  for  a 
randomly-moving  process  to  deviate  by  a  certain  amount  is  proportional  to  the  deviation 
squared. 

The  total  internal  quality  cost  per  product  (CIQ)  is  now  modeled  as  the  sum  of 
the  above  three  costs  —  cost  of  diagnosis  (C/Ql),  cost  of  process  adjustment  (CIQ2),  and 
internal  failure  cost  (C1Q3), 


n  B        C       Ah'1 
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4.4   External  Quality  Cost  Model 

External  quality  cost  is  assumed  to  consist  of  only  one  of  the  components  of  the 
quality  costs  described  in  section  4.2.  It  is  the  component  of  operation  and  maintenance 
cost  caused  by  the  variations  (deviation)  in  manufactured  quality  from  the  design 
specifications.  Variations  in  many  design  parameters  of  the  product  from  their  optimum 
values  cause  the  product  to  perform  suboptimally.  In  the  case  of  the  journal  bearing 
it  may  be  the  clearance  (i.e.  the  bearing  and  journal  diameters);  the  surface  roughness 
and  other  parameters  may  also  cause  external  quality  costs  to  be  affected.  Here  only 
the  variations  in  clearance  are  assumed  to  cause  external  quality  costs.  We  know  that 
operation  cost  (C0)  of  the  bearing  is  a  function  of  the  initial  clearance  C,  (i.e.  for  a  new 
bearing), 

C0  =  f(Q.  (4.16) 

Let  the  actual  initial  clearance  for  the  new  bearing  vary  in  the  range  C/m±Ac  with 
standard  deviation  ac  due  to  manufacturing  variations.  Hence  the  expected  value  of  the 
operation  cost  can  be  approximately  given  by 


1  &A.C) 
E[C0]  *  ACJ  *  ±JLjL 

2  dC] 


<  (4.17) 


Since   it  can  generally   be  assumed   that  C,  is  a  random  variable   having   normal 
distribution,  its  variance  in  terms  of  manufacturing  tolerance  can  be  given  by 
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ac2  =  Ac2/9. 


Hence  the  expected  value  of  the  operation  cost  can  be  expressed  as 


o>2    d2C(C) 
E(C)  =  C  (C/)  +  _5 °1jL 

°       °  2      5C2 


(4.18) 


|c,=c,- 


where  aCi2  is  the  variance  of  initial  clearance  from  the  manufacturing  process.    The 
clearance  in  the  bearing  is  the  difference  between  the  bearing  and  journal  radii 


hence  we  have 


Assuming  that 


C  =  DJ2  -  D/2, 


ac2  =  aDb2/2  +  aDj2/2. 


2  _  2  _  2 

CTDb     -  0"Dj     -   CTD 


(4.19) 


(4.20) 


(4.21) 


and  that  the  diameter  D  is  normally  distributed,  in  the  worst  case  we  can  assume  that 
2Ad=6ctd,  i.e. 

ac2  =  aD2  =  AD2/9.  (4.22) 


Finally    the    external    quality    cost    from    variations    in    bearing    clearance    due   to 
manufacturing  deviations  can  be  expressed  as 
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E(Ca) 


&2d   #Co(C) 


(4.23) 


18 


dC 


CfC,' 


Note  that  the  operation  cost  curve  around  optimum  initial  clearance  C*  must  be  known 
before  external  quality  cost  could  be  calculated. 


4.5   Optimal  Quality  Level 


Now  that  the  internal  and  external  quality  costs  have  been  modeled,  we  can 
determine  the  optimum  quality  control  parameters  by  setting  the  derivative  of  C,Q  with 
respect  to  these  parameters  as  zero.  This  simple  process  gives  the  following  optimum 
values  for  the  quality  control  parameters.  The  optimal  diagnosis  interval  for  the 
production  process  obtained  by  substituting  the  expression  for  u  from  equation  4.14  into 
equation  4.15  and  setting  its  derivative  with  respect  to  n  to  zero.  This  results  in 
optimum  process  measurement  interval  ri  as  follows 


n 


2m  B 


1/2 


(4.24) 


Optimum  quality  control  limit  during  the  production  is  derived  by  differentiating  the 

total  quality  cost  function  in  equation  4.15  with  respect  to  8  and  setting  it  to  zero.   This 
results  in  optimum  quality  control  limits  (D0±8*)  as  follows 
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8*  = 


3C5oA20 
Au 


1/4 


(4.25) 


The  predicted  optimum  manufacturing  process  capability  index  Cp,  defined  as  the  ratio 
of  design  tolerance  (2A)  and  six  sigma  (6a)  of  the  manufacturing  process  is  given  by 


C    = 


2A 


f/iM       16* 

+T 


1/2 


(4.26) 


Substituting  the  expressions  for  u    in  equation  4.24  and  8*  in  equation  4.25  into 
equation  4.15  gives  the  following  expression  for  optimum  internal  quality  cost 


(4.27) 


where  the  two  constants  are  as  follows 


**" 


and 


1/2 

' 
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(4.28) 


K2-± 

192       u 


1 

2 


(4.29) 


The  methods  of  quality  control  presented  in  this  chapter  are  referred  to  as  on-line 
quality  control.  Unlike  the  traditional  practice  of  the  quality  control,  here  the  concept 
of  feedback  in  a  quality  control  system  is  used.    We  determine  the  optimal  control 
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limits  for  the  quality  characteristic  of  interest,  checking  intervals,  and  process  capability 
index  that  minimize  the  total  cost  of  production. 


CHAPTER  5 
THE  DYNAMIC  PROGRAMMING  APPROACH 


The  decisions  of  normal  product  development/realization  efforts  at  each  stage 
of  the  product  life-cycle  are  sequential  in  nature.  The  output  of  one  department  (stage) 
is  the  input  to  next  department  (stage).  Often  the  decisions  made  in  one  department  are 
constraints  to  the  other  department.  Each  department  has  different  immediate  concerns 
in  realizing  a  design,  but  ultimately  have  the  same  goal  —  build  the  best  product  at  the 
lowest  possible  cost.  Each  department  tries  to  minimize  the  cost  of  its  own  separate 
task,  but  one  unit  of  cost  reduction  in  one  department  may  not  necessarily  lead  to 
decrease  in  overall  costs.  Thus  individual  departmental  cost  reduction  decisions 
independent  of  other  departments  will  in  general  lead  to  suboptimal  life-cycle  product 
design. 

5.1    Dynamic  Programming  Problem  Formulation 

The  serial  nature  of  the  problem,  lends  itself  to  the  use  of  dynamic  programming 
solution  strategy.  Dynamic  programming  is  particularly  suitable  where  the  decisions 
taken  at  one  stage  influence  only  'downstream'  stages.  The  dynamic  programming 
method  is  an  efficient  computational  solution  strategy  for  this  kind  of  serial  problems. 
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Because  the  optimization  is  carried  out  stage  by  stage,  dynamic  programming  can  be 
considered  a  decomposition  technique.  Following  Bellman's  principle  of  optimality 
[Bell57]  each  decision  made  at  any  stage  takes  into  account  the  effects  of  that  decision 
on  downstream  stages. 

The  drawback  of  the  dynamic  programming  is  that  it  is  very  difficult  to  deal 
with  the  global  constraints  in  the  problem.  However,  it  has  the  advantage  of  not 
requiring  the  continuity  of  the  objective  function.  Dynamic  programming  is  the  basic 
approach  used  in  this  research  for  achieving  the  global  optimum. 

The  block  diagram  in  Figure  1 .2  represents  the  decision  making  process  for  the 
life  cycle  of  a  product  as  a  set  of  serially  connected  stages.  The  output  of  one  stage  is 
the  input  to  the  next  stage  of  the  decision  making  and  can  have  great  influence  on  the 
subsequent  stages.  Often  the  parameters  of  one  stage  act  as  constraints  to  another  stage, 
thus  restricting  the  design  space.  In  the  block  diagram  in  Figure  1 .2,  there  are  /  stages 
in  the  life  cycle;  each  stage  i  being  represented  by  a  set  of  state  variables  (or  input 
parameters)  Sh„  a  set  of  stage  decision  variables  (or  design  variables)  Xt,  a  set  of 
transformation  functions  F*,  an  objective  function  Ut ,  and  a  set  of  output  parameters 
S, 

The  stage  behavioral  model  for  each  stage  is  described  by  its  transformation 
equations  that  map  the  stage  input  parameters  and  stage  decision  variables  to  the  stage 
output  parameters.  The  stage  i  receives  input  parameters  S,_,  in  conjunction  with  stage 
decision  variables  X,  and  it  produces  the  stage  output  St.  In  general,  this  relationship 
is  described  by  a  set  of  transformation  equations  and  can  be  expressed  as 
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r,<5A/^/)  -  0,     <=lv.,/  (5.1) 

If  the  transformation  function  can  be  explicitly  expressed,  then  the  transformation 
equations  for  stage  i  can  be  written  as 

S,  =  F@„JQ,  (5.2a) 

or  the  relationship  for  f th  output  of  stage  i  can  be  expressed  as 

S^F^.X),   forj=\,..J,  (5.2b) 

where  Ji  is  the  number  of  transformation  equations  (and  output  variables)  for  stage  /'. 
The  objective  function  at  each  stage  is  the  cost  component  associated  with  that  stage. 
For  example,  the  design  stage  objective  function  is  the  cost  of  R&D  and  product 
development,  and  the  operation  stage  objective  function  is  the  cost  of  operation,  support 
and  maintenance  for  the  life  of  the  product.  The  stage  objective  function  can  be 
expressed  as 

U,  =  U,(S,_„  X)  (5.3) 

Optionally,  a  stage  decision  could  be  constrained  by  its  given  input  and  required 
output.    These  constraints  could  be  expressed  as 

Gm  =  GJfiu.  X,  5,)  <  0,   for  m=\,...,M  (5.4) 
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The  transformation  equations,  the  objective  functions,  and  the  constraints  could 
be  of  simple  closed  form,  or  may  require  numerical  simulation,  or  may  be  empirical 
charts.    Finally,  the  overall  objective  function  for  all  the  stages  combined  is  the  life- 
cycle  cost  (LCC),  which  can  be  expressed  as 

U  =  U(S0,  X„...JC„...,  X,),  i=l,..,I  (5.5) 

In  the  case  of  life-cycle  optimization,  the  overal  objective  function  is  additive 

U=l  UfS^Xb        i=l,...,I  (5-6) 

Though  the  problem  of  life-cycle  cost  optimization  has  been  described  as  serial, 

the  stages  are  highly  coupled  in  nature.     A  special  method  which  can  be  used  to 

decompose  (uncouple)  the  serial  problem  is  the  use  of  dynamic  programming.    The 

basic  concept  of  dynamic  programming  is  known  as  Bellman's  principle  of  optimality 

[Bell57]: 

An  optimal  sequence  of  decisions  in  a  multistage  decision  process 
problem  has  the  property  that  whatever  the  initial  stage,  state  and 
decisions  are,  the  remaining  decisions  must  constitute  an  optimal 
sequence  of  decisions  for  the  remaining  problem,  with  the  stage  and  state 
resulting  from  the  first  decision  considered  as  initial  conditions. 

The  following  is  the  general  formulation  of  dynamic  programming  problem.   A  system 
with    /   discrete  stages  can  be  represented  by  a  cascade  block  diagram  as  shown  in 
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Figures  1.2  and  1.3.   The  design  vector  (decision  variables)  and  the  input  vector  (state 

variables)  for  rth  stage  can  be  represented  as, 

X,  ~  {^w—>  Xi,»—>  Xum*  (?•') 

s,-i  =  {Si-u<  •••>  Si-is  ■■•>  Si-u,}>  (5-8) 

where  N,  is  the  total  number  of  design  variables,  Jt  is  the  total  number  of  state  variables 
for  the  z"th  stage. 

The  dynamic  programming  solution  procedure  is  started  with  the  last  stage  and 
applied  at  each  stage  going  backward  until  the  first  stage.  The  optimum  values  of  the 
decision  variables  at  a  stage  depend  on  the  values  of  the  state  variables  at  that  stage  and 
the  effect  of  the  stage  decision  on  subsequent  stages.  We  define  an  optimal  value 
function  for  stage  i,  Ri(SiA),  which  represents  the  optimized  value  of  the  combined 
objective  function  from  stage  i  to  the  last  stage  /.  Evaluation  of  optimal  value  function 
assumes  that  optimum  decisions  will  be  made  in  stage  i  and  all  the  subsequent  stages. 
In  other  words,  R^S^)  is  a  function  denoting  the  optimum  value  of  the  objective 
function  of  a  dynamic  programming  problem  consisting  of  stages  i  to  /  with  the  input 
to  stage  i,  being  given  by  vector  S^.  Therefore  optimum  value  function  at  a  stage  is 
only  a  function  of  the  state  variable  at  that  stage 


RfJSM)  =       Min       £    U#5M,  X)  (5-9) 
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The  selection  of  values  of  the  design  variables  Xi  for  each  stage  /'  and  initial  state  SiA 
is  based  on  the  optimal  objective  function  from  the  final  stage  upto  /' th  stage.  Now  we 
define  a  function  Q„  which  is  the  sum  of  the  objective  function  for  stage  i  and  the 
optimum  value  function  for  stage  i+l,  as  follows 

QiS^X,  5,)  =  UfiSuJQ  +  R^iS.)  (5.10a) 

Recognizing  from  equations  5.1  or  5.2  that  S,  is  a  function  of  SH  and  Xt,  the  above 
relationship  can  be  expressed  as 

QiS^X,)  =  UHS,„XD  +  RM(FkS,M)  (5-lOb) 

which  is  the  sum  of  objective  function  at  stage  i  and  the  optimum  value  function  for 
stage  i+l.  Thus  Q,  is  a  function  of  the  input  state  variables  and  the  decision  variables 
for  that  stage.  The  optimum  value  function  i?,  for  stage  i,  for  a  given  (fixed)  set  of  state 
variables  St,  is  obtained  by  minimizing  Q,  as  a  function  of  decision  variables  Xf.  Hence 
minimum  value  function  Rt  itself  is  a  function  of  state  variables  at  stage  i  and  can  be 
expressed  as 


Rl(S,.l)  =  MmQ,(S,l,X)  (5.11) 
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Now  the  basic  dynamic  programming  recursive  relationship  controlling  the  selection 
of  values  for  design  variables  at  z"th  stage  is  obtained  as  follows 


R{Shl)  =  Min  [UfiS„,  X)]  +  RM(S)  (5-12) 

X 


The  problem  in  equation  5.12  consists  of  choosing  optimum  values  (X',)  for  the 
decision  variables  such  that  the  combined  objective  function  from  stage  i  to  /  is 
optimized.  Thus  the  recursive  dynamic  programming  solution  procedure  given  by 
equation  5.12  proceeds  from  the  last  stage  backwards  and  ends  at  the  first  stage.  At  a 
stage  i  under  consideration,  the  optimum  selection  of  decision  variables  Xt  depends  on 
the  value  of  state  variables  5W  (input  to  stage  /)  and  the  effect  of  optimization  of  the 
stage  on  the  subsequent  stages.  Therefore,  in  order  to  make  optimum  decisions  in 
equation  5.12  at  stage  i,  we  must  know  the  value  of  the  state  variables  S,.,  and  the 
optimum  value  function  for  the  next  stage  (Ri+f)  as  a  function  of  its  input  (S).  Since 
the  value  of  state  variables  S^,  is  not  known  in  advance,  a  reasonable  number  of 
possible  values  within  a  certain  range  must  be  considered.  Hence  a  reasonable  number 
of  dicretizations  of  the  range  of  possible  state  variables  at  the  beginning  of  each  stage 
are  considered  and  the  corresponding  objective  function  evaluated  for  each  stage. 
Specifically,  let  Simm  and  Simax  denote  the  range  of  possible  values  for  state  variables 
at  stage  r+1  as 

Si,min  =   [Sil.mirt    •■    ^ij.mm    ■■    ^Ui.minh  (5.13a) 

Si,max  =   [Sil.max    ■•    ^ij.max-      •    ^Ui.maxh  (5.13b) 
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such  that  the  /  th  variable  lies  within  the  range 

S-.      <S  <S  (5.14) 

Here  a  reasonable  Minimum  set  Simin  and  a  Maximum  set  Sima  for  the  state  variables 
must  be  estimated  in  advance. 

5.2   Dynamic  Programming  Algorithm 

The  dynamic  programming  recursive  relationships  described  above  by  equation 
5.12  is  implemented  through  the  tabular  computations  as  shown  in  Tables  5.1  and  5.2. 
Table  5.1  lists  important  variables  computed  for  stage  i  while  traversing  the  stages 
backwards.  After  having  constructed  a  table  for  each  stage  (total  /  tables),  we  trace 
back  the  path  for  overall  optimum  solution  from  stage  1  forward  to  stage  /.  Table  5.2 
enumerates  the  important  variables  computed  for  stage  /'  while  tracing  back  the  overall 
optimum  solution  and  traversing  the  stages  forward.  Tables  of  type  5.1  and  5.2  will  be 
called  recursion  tables  and  trace-back  tables  respectively. 

5.2.1    Construction  of  Recursion  Tables 

The  following  is  a  description  of  tabular  computations  and  construction  of 
recursion  table  for  stage  i  as  given  in  Table  5.1.    The  description  of  recursion  table  for 
stage  i  is  based  on  the  assumption  that  the  tabular  computations  for  stage  /through  r+1 
have  already  been  performed.    For  stage  i+1,  each  component  of  state  variable  vector 
Si  is  linearly  discretized  within  the  range  described  by  equations  5.13  and  5.14.   If  Ky 
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such  discretizations  are  done  for  the/th  component  of  S,  within  the  range  (Sijmi„,  SiJmax), 

then  the  y"th  discretization  is  given  by 


s         s        +     &■-»     >.min  ,k_X)      k  m  !     r  (5.15) 

ij,k  ij,mm  K  -\  9 

ij 


and  the  total  number  of  the  combinations  P,  of  such  discretization  is  given  by, 


P.  =  ri  K  (516) 


>i     v 


Hence  the  recursion  table  shows  P,  such  discretizations  of  the  state  variable  vector  at 
stage  i+1.    For  stage  i,  the  state  variable  vector  discretizations  S,,,,,,  given  in  column 
(2),  are  used  as  input  to  the  stage  for  evaluating  the  partial  objective  function  U,(Sh 
/  JQ.  The  function  Qt  described  by  equation  5.10  is  used  as  the  full  objective  function 
for  local  stage  optimization.    For  the  /?'th  dicretization,  it  can  be  expressed  as 

Min  Q,,p{S,.llPX)  =  Min  [  UjS^X,)  +  R^Ftf^X))]  (5.17) 

The  crucial  step  now  is  to  evaluate  the  value  of  optimum  value  function  for  stage  z'+l, 
which  is  the  second  term  on  the  right  hand  side  in  equation  5.17,  and  the  column  (6) 
of  the  recursion  table.  After  computing  St  in  column  (5)  (or  F^S^^X^)),  the 
corresponding  optimum  value  function  Ri+/p  is  computed  using  already  computed 
recursion  table  for  stage  i+1.   To  do  this,  state  variable  vectors  S,  given  in  column  (2) 
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and  the  corresponding  optimum  value  function  Ri+I  given  in  column  (7)  of  the  recursion 
table  for  stage  i+l  are  used.  The  recursion  tables  for  stage  r+1  would  only  contain 
discrete  values  of  S,  vectors  in  column  (2),  which,  in  general,  will  not  match  the  desired 
vector  Si+l  in  column  (5)  of  the  recursion  table  for  stage  i.  Thus,  in  order  to  compute 
the  optimum  value  function  for  stage  i+l,  a  general  multi-dimensional  interpolation 
scheme  is  used  and  is  stated  in  the  following  paragraphs. 

5.2.2   Multi-Dimensional  Interpolation 

The  multi-dimensional  interpolation  procedure,  which  was  implemented  for  a 
very  general  case,  uses  quadratic  and  cubic  interpolations.  The  multi-dimensional 
interpolation  problem  can  be  stated  as  follows.  Let  v  be  a  function  of  a  J  dimensional 
vector  u  defined  by  some  unknown  function  <|>, 

v  =  <K«)  (518) 

where 

u  =  {u„u2,  ...  uj}.  (5.19) 

Let  all  discretizations  of  /th  component,  up  be  stored  in  vector  uj  as  follows 

HJ  =  ^ujh..~ujk,..MjK),     for  y=l,..  J  (5.20) 

which  are  stored  in  a  PxJ  matrix  U  as  row  vectors.   The  corresponding  function  values 
(at  nodal  point  p)  vp  are  stored  in  a  vector  V,  such  that 


88 
v,  =  4>(«,),   for  p=\,..P.  (5.21) 

In  the  interpolation  problem  encountered  in  the  general  dynamic  programming  problem, 
matrices  U  and  U,  and  vector  V  are  known.  The  objective  is  to  evaluate  the  unknown 
function  v*=<|)(i/*)  for  a  given  vector  u  using  interpolation  of  surrounding  nodal  points 
in  J  dimensional  space.  In  this  thesis,  both  quadratic  and  cubic  polynomial  multi- 
dimensional interpolations  were  implemented  and  are  described  in  Appendix  A. 

When  the  recursion  table  are  constructed  while  traversing  backwards  through  the 
stages,  the  interpolations  basically  map  the  columns  (2)  and  (7)  of  stage  z'+l  onto  the 
columns  (5)  and  (6)  of  stage  i.  According  to  equation  5.12,  the  optimum  value  function 
for  stage  i,  in  column  (7)  of  the  recursion  table,  is  obtained  by  simply  adding  column 
(4)  and  (6).  J?,  in  column  (7)  also  serves  as  the  objective  function  at  stage  i.  Since  the 
input  state  variables  S0  for  stage  1  are  prescribed  for  the  problem,  no  discretizations  for 
state  variable  vector  are  needed.  Hence  the  recursion  table  for  stage  1  will  always  have 
only  one  row  in  it  (i.e.  K,=J,=\). 

5.2.3    Stage  Optimization 

At  each  stage  of  the  dynamic  programming  solution  approach,  we  encounter  the 
optimization  problem  stated  by  equations  5.10  to  5.12.  The  optimization  is  carried  out 
subject  to  the  M,  inequality  constraints  given  by  equations  5.4.  The  optimization 
problem  for  stage  /,  with  decision  variables  being  Xi  and  the  parameters  being  ShI  „  is 
restated  as  follows 
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minimize  Q./S.^X,)  =    Utf^X)  +  R^Ftf^X))  (5.23) 

subject  to  Gim  =  GaS,,,rX>S)  <  0,   for  m=l,..,M,  (5.24) 

where,  from  equations  5.7  and  5.8,  decision  variable  vector  and  state  variable  vector 
respectively  are 

X,  =  {Xu,  ..  X,„  ..  Xm},  (5.25) 

St-i.p  =  {Si-i,p,i>  ■■  ^i-i.pj>  ■■  Si-i,p.ji}>  (5.26) 

where  N,  is  the  total  number  of  design  variables,  Jt  is  the  total  number  of  state  variables, 
and  M,  is  the  number  of  constraints  for  the  /'th  stage.  The  stage  optimization  problem 
in  equations  5.23  -  5.26  is  stated  for  the  /?'th  combination  of  the  state  variable  vector 
(nodal  point/?),  with  p=\,..Pj.  This  is  a  general  constrained  non-linear  optimization 
problem  and  since  it  is  also  a  part  of  multi-level  optimization  approach,  it  will  be 
discussed  in  detail  in  the  following  chapter. 

5.2.4   Construction  of  Trace-Back  Tables 

After  the  backward  recursive  traversal  procedure  has  been  completed,  and  as  a 
result,  the  recursion  tables  have  been  constructed,  the  overall  optimum  solution  is 
obtained  by  tracing  back  the  existing  recursion  tables  in  the  forward  direction  from  the 
first  stage  through  the  last  stage.  At  stage  1,  the  state  variables  (input  vector  S0)  is 
known  in  advance,  hence  the  recursion  table  contains  only  one  row  (since  there  are  no 
state  variable  discretizations  needed).    At  stage  1,  the  optimum  decision  variables  X' , 
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in  column  (3)  of  recursion  table  (Table  5.1)  are  the  overall  optimum  decisions,  X*„ 

since  they  lead  to  overall  optimum.  The  selection  of  X*,  in  stage  1  leads  to  optimum 
state  variables  S',  at  stage  2.  Subsequently  the  stage  optimization  of  equation  5.10  for 
stage  2,  with  S',  as  input,  will  lead  to  overall  optimum  decisions,  X"2,  and  optimum 
state  variables  S'2  at  stage  3.  In  general,  for  a  stage  i  we  would  know  the  optimum 
state  variables  S*,.,,  and  the  stage  optimization  will  lead  to  overall  optimum  decisions, 
A^,  and  optimum  state  variables  £*,,  As  discussed  earlier,  the  stage  optimization 
problem  in  equations  5.23  -  5.24  would  require  the  optimum  value  function  for  stage 
z'+l  (/?;+/).  In  general,  the  state  variables  5^,  evaluated  (needed  to  evaluate  J?/+/)  will  not 
be  found  in  recursion  table  for  stage  i.  To  obtain  the  optimum  value  function  Ri+„ 
multi-dimensional  interpolation  is  again  invoked.  Note  that  while  tracing  back  the 
optimum  solution,  optimization  is  done  only  for  one  set  of  state  variables  at  each  stage, 
as  against  P,  sets  of  optimization  needed  while  constructing  the  recursion  tables. 

This  process  of  optimization  for  stage  i  is  explained  by  constructing  the  trace- 
back  table  for  each  stage  as  described  by  Table  5.2.  Column  (2)  of  the  trace-back  table 
contains  the  output  state  variable  S^,  from  the  previous  stage  (i-l).  After  stage 
optimization,  the  optimum  decision  variables  X \  are  entered  in  column  (3);  the 
optimum  objective  function  lfj  is  entered  in  column  (4);  the  resulting  output  state 
variables  S*t  are  entered  in  column  (5);  the  corresponding  optimum  value  function  for 
the  next  stage  Rj+„  obtained  by  interpolation,  is  entered  in  column  (6).  Column  (7) 
shows  the  sum  of  column  (4)  and  (6).  By  tracing  back  in  this  manner,  the  over  all 
optimum  solution  to  the  dynamic  programming  problem  is  obtained.    The  optimum 
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value  function  for  stage  i+\  (Rj+I  in  column  (6)  of  Table  5.2)  is  evaluated  by  multi- 
dimensional interpolation  as  described  in  detail  in  Appendix  A. 


Procedure  6.1:  Discretize 

{ 

input  J0 

input  vector  S0[\.J0]; 

for  stage  z-1  to  /  { 

input  number  of  state  variables  J:  for  stage  i; 

input  vectors  S^JLJ^  S^Jl.JJ  for  stage  i; 

input  number  of  discretizations  Lt  at  stage  i; 

forj-l  to  Jj 

create  discretization  points  Su; 

evaluate  Kj> 

I*  Create  discretized  nodal  points  as  follows.  */ 

*=1; 

*«*JM]  =  ConvBase(fc,  L,  J,); 

fory-1  to  Ji 

•*  i.Aj  =  $ ijXkindaU]]' 

for  nodal-point  &=2  to  Kt 

{ 

kindeA\.J]=Add\(L>  kMex); 
for  y'=l  to  Ji 

/*  define  the  nodal  points  matrix  S'  as  follows.  */ 

"•  W  =  $  ijikjndeXiW; 

}; 
}; 
}; 


Procedure  6.2:  DynamicProg 

{ 

call  Discretize; 
for  stage  i=\  to  / 

{ 

for  nodal-point  k=\  to  Kt 

{ 

Optimize_ND_Cnstr  (Functionjg,  i,  S„  M,  xT,  Qp  Qe  X\b  R'lh 

AL,  Q,„  Si-liJt,FiJ,lh  CV); 
I*    Where  the  input  and  output  to  the  procedure  are  as  follows: 
Input: 

function_Qs      Procedure  to  evaluate  Qt. 
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/  Stage. 

•£/.*[  1"«^]  Nodal  point  k  for  state  variable  vector. 

M,  Number  of  decision  variables  at  stage  i. 

A°[1..MJ  Initial  value  of  decision  variable  vector  at  stage  i. 

Qu  Number  of  inequality  constraints  at  stage  /'. 

Qie  Number  of  equality  constraints  at  stage  i. 

Ouput: 

X'jk[\.J]  Resulting  optimum  decision  variable  vector  at  stage 

/  for  £'th  nodal  point. 
R  \k  Resulting  optimum-value-function  for  stage  i  for  £'th 

nodal  point. 
AL  Augmented  Lagrangian  function. 

Q* 

Si-l,k 

FU,,k 

CV  Constraints  violations  array.    The  array  contains  the 

list  of  constraints  (constraint  numbers)  violated  at  the 

optimum  point. 


}; 


Procedure  6.3:  Optimize  NDCnstr  (junction^,  i,  S,  M„  A°,  Q>  Q0  X\b 

R'lMAL,Q,kSi-liMFi_l,iM,CV) 


{ 


Procedure  Input: 

function_Qj  Procedure  to  evaluate  Qt. 

i  Dynamic  programming  stage  number. 

J_  Number  of  parameters  (dimensionality  of  5). 

S[\.J]  Parameter  variable  vector. 

N  Number  of  decision  variables  (dimensionality  of  X). 

A"[1..JV]  Initial  value  of  decision  variable  vector. 

Qi  Number  of  inequality  constraints. 

Qe  Number  of  equality  constraints. 

Procedure  Ouput: 

X*[\..N]  Resulting  optimum  decision  variable  vector  at  stage  /'  for  Fth 

nodal  point. 

R'  Resulting  value  of  optimum  function  (F). 

AL  Augmented  Lagrangian  function  =  R'  +  P';  Where  P*  is  the 

penalty  function  at  the  optimum. 
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Q.  Value  of  optimum-value-function  for  stage  i  (combined  optimum 

objective  function  from  stage  /  to  /). 
nextSj  Output  state  variables  vector  (stage  i+l  input). 

nextQ,  Value   of  optimum-value-function    for   stage    i+\    (combined 

optimum  objective  function  from  stage  r+1  to  I). 
CV  Constraints  violations  array.     The  array  contains  the  list  of 

constraints  (constraint  numbers)  violated  at  the  optimum  point. 

Procedure  Approach: 

The  procedure  uses  a  constrained  non-linear  optimization 
approach  to  optimize  a  given  function  of  decision  variable  vector 
X[l..N\  and  parameters  vectors  S[\.J]  . 

}; 


Procedure  6.4:  ConvBase(kl0,  b,  m) 

{ 

Given  a  number  k,„  in  decimal  base,  convert  it  into  base  b.    It  returns  an  array 

£A[l..m],  of  length  m,  containing  the  converted  number. 

}; 


Procedure  6.5:  Addl(b,  m,  kj 

{ 

Given  a  number  kb  in  base  b  (an  array  of  length  m),  increment  the  number  by  1 

(add  1)  and  return  it  (an  array  fcA[l..m],  of  length  m. 

}; 


CHAPTER  6 
THE  INTEGRATED  DESIGN  DECISION  SYSTEM 


For  the  life-cycle  optimization,  three  different  solution  approaches  are 
investigated  and  implemented.  The  first  two  approaches  tried  are  algorithmic  and  the 
third  approach  combines  an  algorithmic  approach  with  a  heuristic  approach.  The 
algorithmic  optimization  approaches  implemented  are  very  general  in  nature. 

The  first  optimization  approach,  as  presented  in  the  previous  chapter,  is  to  treat 
the  life-cycle  optimization  as  a  dynamic  programming  problem.  A  very  general 
nonlinear  constrained  optimization  module  was  developed  which  can  handle  a  design 
vector  of  arbitrary  size  and  any  number  of  constraints.  The  constrained  optimization 
problem  is  first  converted  into  an  unconstrained  optimization  problem  by  adding  a 
combination  of  interior  and  exterior  penalty  functions  to  the  objective  function,  and  then 
using  augmented  Lagrange  multiplier  method.  The  unconstrained  problem  is  optimized 
using  conjugate  directions  or  variable  metric  methods. 

The  second  optimization  strategy  is  a  multi-level  optimization  approach,  where 
the  design  optimization  is  decomposed  into  many  levels.  Here  the  design  vector  is 
decomposed  into  a  global  design  vector  and  level-specific  design  vectors.  At  each  level 
of  optimization,  calls  to  nonlinear  constrained  optimization  modules  are  made. 


95 


96 
The  third  optimization  approach  is  a  combination  of  heuristic  and  algorithmic 
procedures.  Generally,  the  combination  of  the  two  procedures  in  a  decision  system  can 
provide  the  best  design  tools.  The  heuristic  system  generates  an  initial  design  based  on 
the  application  domain  of  the  bearing  specifications  supplied  by  the  designer.  The 
heuristic  module  can  also  process  the  feedback  about  constraint  violations  and  modify 
the  appropriate  parameters  so  that  the  design  remains  within  the  feasible  domain. 

6.1    Multilevel  Optimization  Approach 

A  general  discussion  of  multilevel  optimization  now  follows.  While  the  well 
known  nonlinear  optimization  techniques  have  proved  effective  where  limited  number 
of  variables  and  constraints  are  involved,  large  number  of  variables  may  require 
excessive  computer  time.  One  remedy  for  this  problem  lies  in  the  resolution  of  the 
large  problem  into  smaller  units.  One  approach,  presented  in  the  last  chapter,  is  the 
dynamic  programming  method,  where  the  problem  is  treated  in  stages  with  optimal 
decisions  being  made  at  each  stage.  Yet  another  approach  to  this  problem  is  based  on 
decomposition  of  the  system  into  a  number  of  smaller  subsystems,  each  with  its  own 
objective  function  and  constraints,  solved  independently  at  the  1  st  level  and  coordinated 
at  the  2nd  level.  General  formulation  of  such  decomposition,  based  on  Kirsch  [Kirs75], 
and  Kirsch  and  Moses  [Kirs79]  is  presented  in  this  sections. 

Consider  the  following  general  optimization  problem  with  decision  variables 
being  the  vector  W 


Minimize 


Subject  to  constraints 


h(W)  =  0 

gW<o 
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(6.1) 

(6.2) 
(6.3) 
(6.4) 


where  W1  are  the  lower  bounds  and  W0  are  the  upper  bound  for  W.  The  multilevel 
optimization  approach  consists  of  decomposition  of  this  general  problem  into 
subproblems,  in  which  so-called  coordination  variables  are  chosen  by  second  level 
controller  to  ensure  solutions  for  the  independent  first-level  system  that  actually 
correspond  to  an  overall  system  optimum.  As  would  be  expected,  success  is  conditional 
on  suitable  conversion  of  the  integrated  problem  to  the  two  level  form,  in  which  the 
first-level  result  must  have  noninteracting  parts  capable  of  separate  treatment.  In  the 
model  coordination  approach,  the  vector  W  is  partitioned  into  subvectors,  T  and  X 


W  = 


Y 
X 


(6.5) 


X  = 


(6.6) 
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in  which   Y  is  the  subvector  of  interactions  or  coordination  variables  between  the 

subsystems,  and  X  is  the  subvector  of  variables  for  subsystem,  with  vector  X„  being  the 

vector  for  n'th  subsystem.    The  partitioning  of  the  decision  variable  vector  W  is  done 

in  such  a  way  so  that  the  equality  and  inequality  constraints  and  the  objective  function 

can  be  rewritten  in  the  form 


h(W)  = 


\{Y,XX) 
h„(W 


(6.7) 


%{W)  = 


g,(iyr,) 


BN(Y,XN) 


(6.8) 


and  the  objective  function  can  be  written  in  the  form 


z  =  nw)  =  j;  fn(Y,xn) 

n=l 


(6.9) 


Based  on  these  assumptions,  the  two-level  optimization  problem  can  now  be  stated  as 
follows. 

First  level  optimization  —  Determine  a  fixed  value  for  Y  through  the  constraints 


f=  r 


(6.10) 
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Then  the  optimization  problem  in  equation  6.9  can  be  decomposed  into  N  independent 

first-level  problems,  each  of  which  is  stated  as  follows.    Find  Xn  such  that: 

Minimize  Z„  =  f„(  Y°,Xn)  (6.11) 

Subject  to  constraints 

hn(F,Xn)  =  0,     n=\,...,N  (6.12) 

aXF^jSO,     n=\,...,N  (6.13) 

XL<X<XU  (6.14) 
If  we  introduce  the  function 

H„(T)  =  min  Z„  (6.15) 

then  the  first-level  problem  consists  of  finding  the  function 

*(h  =  E   #  (*0  ■  E   min  Z  (6-16) 

for  all  ^V„  (n-\,...,N)  satisfying  the  constraints  given  in  equations  6.12  to  6.14. 

Second  Level  Optimization  —  The  task  in  the  second  level  problem  is  to  find 
Y°,  which  minimizes 

H{f)=J:Hn{f)  (6-17) 

Now  the  two  level  optimization  problem  can  be  solved  iteratively  by  the  following 
procedure  (see  Figure  6.1): 
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1 .  Choose  an  initial  value  for  the  coordination  variable,  T°. 

2.  For  a  given  Y°,  solve  the  N  independent  first-level  problems. 

3.  Modify  the  value  of  Y°  so  that  the  value  of  H{Y°)  is  reduced. 

4.  Repeat  steps  2  and  3  until  minimum  H{Y°)  is  achieved. 

The  multilevel  optimization  approach  formulated  in  this  section  is  applied  for  the 
optimum  design  of  journal  bearings  in  the  next  section. 

6.2   Parameter  and  Tolerance  Design  By  Multilevel  Optimization 

Minimizing  the  life  cycle  cost  is  a  coupled,  non-linear  constrained  optimization 
problem.  The  life  cycle  cost  consists  of  manufacturing  cost,  operation  and  maintenance 
cost,  internal  and  external  quality  cost.  The  objective  is  to  minimize  the  total  cost  rate 
(life-cycle  cost  per  hour  of  bearing  operation),  which  can  be  defined  as  follows 


where 


Cr  =  CJt,  =  (CM  +C0  +  CIQ  +  CEQ)lt,  (6. 1 8) 

manufacturing  cost  =  CM  =  C^L,D,AD),  (6.19) 

operation  cost  =  C0  =  C0(L,D,C,fi),  (6.20) 

internal  quality  cost  =  CIQ  =  CIQ(AD),  (6.21) 


external  quality  cost  =  CE    =  C 


d2C 


dC- 


=  C£„(C,AD,C;)    (6.22) 


bearing  life  =  t,. 
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The  optimum  design  is  chosen  subject  to  constraints  on  bearing  design  given  by 

equations  2.1  to  2.8.    The  decision  variables  being  bearing  length  L,  initial  clearance 

Cj,  oil  viscosity  |a,  and  design  tolerance  on  bearing  and  journal  diameter  AD.    It  is  a 

coupled  problem  since  the  external  quality  cost  CEQ  is  a  function  of  double  derivative 

of  operation  cost  function  C0(C). 

The  optimum  bearing  design  problem  is  decoupled  by  dividing  the  optimization 

problem  into  two  levels.    The  schematic  diagram  in  Figure  6.2  presents  the  algorithm 

for  the  two  level  optimum  parameter  and  tolerance  design  for  the  journal  bearing.    In 

the  first  level  we  make  appropriate  assumptions  for  the  values  of  Ap  and  C0"(Cj)  and 

optimize  for  L,  C(,  and  \i.  After  this  optimization  the  operation  cost  curve  as  a  function 

C(  is  available  around  the  optimum  initial  clearance  Cj*.    Thus  C^  can  be  found  out 

more  accurately.     Now  we  optimize  the  total  life  cycle  cost  C|C  as  a  function  of 

tolerance  AD  only.  With  this  optimum  value  of  AD\  we  can  again  go  back  to  first  level 

optimization.    Comparing  with  the  notation  introduced  in  the  general  formulation  of 

multilevel  optimization  problem  (Figure  6.1). 

Y=[C0"]  (6.23) 

X,  =  [L,D,C,fi,AD]  (6.24) 

W  =  [Y,  X,]  =  [C0",L,D,C^,AD]  (6.25) 

The  objective  function  is  evaluated  by  invoking  the  bearing  simulation,  which  not  only 
evaluates  the  operation  cost  but  also  determines  the  life  of  the  bearing. 


Given:  Bearing  parameters;  wear  parameters;  cost 
coeffiecients  for  manufacturing,  quality,  and  operation. 
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I 


•  Assume  initial  values  of  the  decision  variables  L,  Q,  /*,  AD. 

•  Evaluate  internal  quality  cost  as  a  function  of  tolerance,  CIQ(AD). 

•  Guess  the  value  of  internal  quality  cost,  C^q 


3 


•  Optimize  cost  rate:  Cr=  LCC/life 
LCC  =    manufacturing  cost  (Cm)  +  operations  and  maintenance 

cost  (C0)  +  internal  quality  cost  (CIQ)+  external  quality  cost 

(Qq)'> 

•  Subject  to  constraints  on: 

outlet  oil  temperature  T0,  maximum  pressure  PmJ?  ,  minimum  film 
thickness  h0 ,  L/D  ratio,  initial  clearance  Q  ,  oil  viscosity  jh,  bearing 
life,  and  dynamic  stability; 

•  Design  variables:  length  L,  initial  clearance  Q,  viscosity  /*. 


I 


•  For  optimum  journal  bearing  (L",  Q'./i'),  evaluate  2nd  derivative 
of  operation  cost  d*C0  /  dCf  at  optimum  clearance  C0"(C,*); 

•  Optimize  the  LCC  (AD)  =  CM(AD)  +  CIQ(AD)  +  CEQ[AD  ,  C^Q*)] 
with  tolerance  AD  bemg  the  design  variable; 

•  Use  C,Q  and  Qq  values  (calculated  based  on  AD*)  for  next  iteration. 


No 


I 


Have  cost  rate  Cr  and  AD*  converged 


lYes 


Optimum  parameter  and  tolerance  designed. 


Figure  6.2:     Optimum  parameter  and  tolerance  design  for  journal  bearing. 
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6.3   Nonlinear  Constrained  Optimization 

This  section  presents  the  approach  used  to  solve  the  general  multivariable 
nonlinear  constrained  optimization  problem.  This  general  optimization  problem  is 
encountered  in  selecting  the  optimum  decision  variables  while  using  two  different 
approaches  in  this  thesis.  First  is  the  problem  of  stage  optimization  in  the  dynamic 
programming  based  approach  presented  in  the  last  chapter.  Second  is  the  multilevel 
optimization  approach  discussed  in  the  previous  section.  This  general  optimization 
problem  is  solved  by  converting  the  constrained  problem  to  an  unconstrained  one,  using 
the  augmented  Lagrange  multiplier  method.  The  unconstrained  problem  is  then  solved 
employing  a  conjugate  directions  algorithm  or  variable  metric  algorithm,  with  one 
dimensional  search  performed  using  method  of  golden  sections.  A  schematic  diagram 
illustrating  the  relationship  between  various  modules  of  the  integrated  system  is  shown 
in  Figure  6.3. 

A  brief  statement  of  the  general  multivariable  nonlinear  constrained  optimization 
problem  now  follow.  Consider  the  general  optimization  problem  of  equations  6.1  to 
6.4,  restated  as  follows,  with  decision  variables  being  the  vector  X. 

Minimize  Z  =  fiX),  (6.26) 

Subject  to  constraints 

hk(X)  =  0,     k=\,..,K  (6.27) 

gJX)<0,    /-W  (6.28) 
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o  N-dimensional  constrained  optimization 

o  Convert  to  unconstrained  optimization  using 
Augmented  Lagrange  Multiplier  Method 


N-dimensional  unconstrained  problem  solved 
using  gradient  search  methods  (Conjugate 
dierctions  or  variable  metric  methods). 


Direct  directional  search  using  method  of  golden 
sections  and  quadratic  polynomial  interpolations 


Bearing  performance  /  wear  simulation  over  its 
physical/economic  life 


Hydrodynamic 
Lubrication  Theory 


Thermo-hydrodynamic 
Lubrication  Theory 


Figure  6.3:     Relationship  among  various  modules  in  integrated  optimization  system. 
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where  the  inequality  constraints  include  the  set  of  side  constraints  of  type  X1  <  X  <  Xu. 

If  we  had  only  the  equality  constraints,  the  Kuhn-Tucker  conditions  require  that  if  we 

create  the  Lagrangian 


L(X,\)  =  fiX)  +  £   Xkhk(X)  (6.29) 

then  the  minimization  of  the  Lagrangian  subject  to  equality  constraints  of  equation  6.27 
also  provides  the  solution  to  the  optimization  problem  given  by  equations  7.26  and  7.27. 
Here  Xk,  for  k=\,..,K  are  the  lagrangian  multipliers.  Thus  we  can  use  the  exterior 
penalty  function  approach  to  create  a  pseudo  objective  function,  called  augmented 
Lagrangian  function  (ALF)  to  be  minimized  as 

A(X,X,rp)  =AX)  +  E   [  M*(*)  +  rJhfiOf]  (6-30) 

An  efficient  method  for  optimization  is  to  repeatedly  optimize  ALF  with  the  X  and  r 
held  constant.  The  new  set  of  augmented  Lagrange  multipliers  is  calculated  using  the 
following  relation  for  updating  the  lagrange  multiplier  [Vand85] 

K]  =  K  *  2rphk(Xn,      for  k=\,...^  (6.31) 

where  p  denotes  the  iteration  number. 

The  optimization  problem  with  both  equality  and  inequality  constraints  can  also 
be  tackled  by  using  the  augmented  Lagrangian  multiplier  (ALM)  method.  Vanderplaats 
[Vand85]  presents  the  following  algorithmic  approach  to  solve  the  problem.  The 
general  augmented  Lagrangian  for  an  optimization  problem  with  equality  and  inequality 
constraints  is  given  by 


j     r  .        a: 

2 


A{XXrp)  =/X)  +  £    A.V>  +  r^    +  E    M#>  +  ',[*#>! 


t-i 


The  update  formulas  for  the  lagrange  multipliers  are  now  given  by 
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(6.32) 


1^*1 


2r 


max 


gjm, 


2r 


for  j=\,..J 


(6.34) 


*C  =  ^  +  2r/t(A),      for  y=l,...^  (6-35) 

The  general  algorithm  for  optimization  by  augmented  Lagrange  multipliers  is  presented 
in  Procedure  Optimize _ND  Constrained.  The  method  has  been  adopted  to  solve 
general  nonlinear  constrained  optimization  problems  because  of  the  following  features. 

1 .  The  method  is  relatively  insensitive  to  the  value  of  the  penalty  multiplier  rp  and 
the  way  by  which  they  are  updated.  The  rp  does  not  have  to  increase  to  oo  for 
obtaining  final  optimum. 

2.  Precise  satisfaction  of  constraints  g/X)=0  and  hk(X)=0  is  possible. 

3.  Acceleration  towards  optimum  solution  is  obtained  by  updating  the  Lagrange 
multipliers. 

4.  The  starting  point  may  be  either  feasible  or  infeasible. 

5.  At  the  optimum,  the  value  of  the  Lagrange  multipliers  A.*,  *  0  will  automatically 
identify  the  active  constraint  set. 

6.  The  availability  of  A.*,  at  the  optimum  point  and  nearby,  can  facilitate  the 
evaluation  of  sensitivity  of  optimum  solution  to  problem  parameters. 
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The  ALF  (overall  life-cycle  cost  including  the  penalty)  is  evaluated  by  the  Procedure 

Augmented JLCCost.     This  procedure  makes  calls  to  Procedure  MIEQCost  in  order 

to  evaluate  the  manufacturing  cost,  internal  quality  cost  and  external  quality  cost.   The 

Procedure   AugmentedLCCost    also   calls   the   procedure   Augmented  Op  Cost    to 

evaluate  the  operation  cost  of  the  bearing. 


Procedure  6.1:  Augmented _LC  Cost  {from,  N,  J,  K,  X,  Cr,  Vt,  V,  A,  X) 

{ 
Purpose: 

Evaluate  augmented  Lagrange  function  (A)  given  by  equation  7.32  and  the  cost 
rate  (Cr  =  LCC/life). 
Input: 

The  procedure  which  called  this  procedure,  from 
Number  of  decision  variables,  N; 
Number  of  equality  constraints,  J; 
Number  of  inequality  constraints,  K; 
Decision  variable  vector  X,  defined  as  follows: 

Ar/=AD/ADm/n,   X2=C0"(C),  X,=bearing  length  L,  JQ=bearing  diameter  D, 

■X-s~*~M>    -X^C/q,   X7=CEq. 

Output: 

Life-cycle  cost,  C,c; 

Life  of  the  bearing,  t, 

List  of  constraint  violations  at  the  beginning  of  life  at  X0,  V;, 

List  of  constraint  violations  at  the  end  of  life  at  X",  V; 

Augmented  Lagrangian  function,  A; 

Lagrange  multipliers,  X. 

Steps: 

1.      Initialize: 

Initial  value  of  Lagrange  multipliers  X°=0; 

Initial  value  of  penalty  multiplier  rp=\; 

Coefficient  by  which  rp  is  incremented  at  each  iteration,  y=2  to  10; 

Maximum  value  of  rp,  rpma*=100  to  10000; 

x  =  x„, 

J        Jo  ~  JLCoA-oi't 
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2.  If  LCcount  =  0,  then  calculate  the  cost  of  manufacturing,  internal  and  external 
quality  by  using  the  procedure  MIEQCost 

MIEQ_Cost{\,X,CmEQ,  V„  V) 

3.  If  the  current  value  of  tolerance  on  clearance  (AD)  is  out  of  its  range  (ADmi„ 
ADmax)  or  if  it  makes  the  clearance  out  of  its  range  (Cmm,  Cmax),  then  adjust  the 
value  of  the  tolerance  as  follows 


mm 


IfC0-Cm//I<AD,then   AD  =  C0-C, 
IfC^-C^A^then   AD  =  Cmax-C0 
If  AD  <  ADmm,  then   AD  =  ADmin 
If  AD  >  ADma,  then   AD  =  ADmax 

4.  Evaluate  three  different  values  of  initial  clearance,  Cja=C0-Az>  C^C^ 
Cic=C0+AD.  Initialize  3  decision  variable  vectors  Xa,  Xb,  and  Xc,  such  that 
Xa,2=CJC>,m»*  Xb,2=Ci,b/Cimin,  Xc2=CJC  imin.  Now  evaluate  the  bearing 
operation  cost-rate  for  them   as  follows 

Operation  Cost  (N,  X„  Cree  C0t„  Cmla,  Cemr^  ha  Via  Va), 
Operation  Cost  (N,  X¥  Cr>  C0>  C0//fl)  Cenergy<r  tkm  V,¥  Vb), 
OperationCost  (N,  X0  C„  Ca„  Cmla  CmergyM  tlo  Vitf  Vc). 

Evaluate  the  operation  cost  rate  at  these  three  points 

^Or.a  ~  (S'oil.a  *- energy.a' )' " la 

Cor.b  ~  \C0il,b   +   ^energy.b't'hb 
Cor.c  ~  \CoUc  +   (-energy.c)'hc 

Given  3  points  for  initial-clearance  versus  operation-cost-rate,  (C/0,  C0ra),  (Cib, 
C0rb)  and  {Cic,  C0rc),  find  the  second  derivative  of  cost-rate  C0r  with  respect  to 
clearance  C,  (C0r.)  assuming  a  quadratic  polynomial. 

5.  Evaluate  the  gradient  of  CM£(?  with  respect  to  tolerance  AD  and  then  the  search 
direction 

s  =  -  gradient/  |  gradient  |  . 

Do  a  1  dimensional  optimization  with  respect  to  tolerance  AD  in  the  direction  s 
using  the  procedure, 

Optimize!  D  (MIEQCostf),  Augmented  LCCost,   7,  Xb,fs,  cw  a,  X*, 
/,  R,  V,  V) 
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6.  If  the  optimum  tolerance  AD*  makes  (C+AD*)  go  outside  the  specified  range  of 
clearance  {Cpmin,  Cimca),  then  bring  it  back  within  this  range  by  adjusting  AD' 
so  that 

C,,m(„  <  (C,-AD*)    and   Cimax  >  (C,+AD*) 

7.  Evaluate  the  Augmented  Lagrange  function  using  the  procedure 

called  Jrom  =  Augmented  LC  Cost 

ALF  =  Augmented  Op  Cost  (called Jrom,  N,  J,  K,  x,  C„  V,  V) 

8.  Exit. 

} 


Procedure  6.2:   MIEQCost  {called Jrom,  N,  X,  CMEQ,  t„  Vt,  V) 

Purpose: 

Evaluate  manufacturing  cost  (CJ,  internal  quality  cost  (C^),  and  external  quality 

cost  (CEQ). 
Input: 

The  procedure  which  called  this  procedure,  called  Jrom 

Number  of  decision  variables,  jV=1  (or  7,  if  calledfrom  =  Optimize  ID); 

Decision  variables  vector,  X: 

X,=A,/ADmi„,  X2=C0"(C,),  A^bearing  length  L,  ^bearing  diameter  D 

Output: 

Decision  variables  vector,  X: 

X,=A,/ADmm,   X,=C0"(Q,  ^=bearing  length  L,  ^=bearing  diameter  D, 

^5=^Kf>     -X6=C,q,    X~=CEq 

Manufacturing  +  internal  quality  +  external  quality  costs,  CmEQ\ 
Bearing  life,  /; 

List  of  constraint  violations  at  the  beginning  of  life  at  X0,  V;, 
List  of  constraint  violations  at  the  end  of  life  at  X*,  V; 

} 
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Procedure  6.3:   Augmented  Op  Cost  {called Jrom,  N,  X,  C0,  t„  V^  V) 

{ 
Purpose: 

Evaluate  manufacturing  cost  (Cm),  internal  quality  cost  (C/?),  and  external  quality 
cost  (CEQ). 
Input: 

The  procedure  which  called  this  procedure,  called  Jrom 

Number  of  decision  variables,  N=\  (or  7,  if  calledfrom  =  Optimize  ID); 

Decision  variables  vector,  X. 

A/=AD/ADtmiI,   X2=C0"{C),  .Yj=bearing  length  L,  X,=bearing  diameter  D 

Output: 

Decision  variables  vector,  X. 

X1=b.JJLDmt>  X,=C0"(C,),  Zj=bearing  length  L,  ^=bearing  diameter  D, 
-X-s~Cm,    X6—Cjq,   A7—CEq 
Manufacturing  +  internal  quality  +  external  quality  costs,  CmEQ\ 
Bearing  life,  /; 

List  of  constraint  violations  at  the  beginning  of  life  at  X0,  V;, 
List  of  constraint  violations  at  the  end  of  life  at  X*,  V; 
} 


Procedure  6.4:    Optimize _ND  Constrained  {Augmented _LCCQ,  N,  J,  K,  X°,f,  X*,f, 
V„  V,  A,  X) 

{ 
Purpose: 

Perform  constrained  optimization  of  a  function  of  X  (a  JV  dimensional  vector). 

Input: 

Augmented  Life-cycle  cost  function  to  be  optimized,  ALCC; 

Number  of  decision  variables,  jV; 

Number  of  equality  constraints,  J; 

Number  of  inequality  constraints,  K; 

Initial  values  for  decision  variable  vector  X°. 

Output: 

Values  of  the  function  at  X°,  f; 

Optimum  values  for  decision  variable  vector  X*; 

Optimum  cost  function,  /; 

List  of  constraint  violations  at  the  beginning  of  life  at  X0,  V;, 
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List  of  constraint  violations  at  the  end  of  life  at  X*,  V; 
Augmented  Lagrangian  function,  A; 
Lagrange  multipliers,  X. 

Steps: 

1.  Initialize: 

Initial  value  of  Lagrange  multipliers  A°=0; 

Initial  value  of  penalty  multiplier  r=\; 

Coefficient  by  which  rp  is  incremented  at  each  iteration,  y=2  to  10; 

Maximum  value  of  rp,  r/^lOO  to  10000; 

X  =  X0; 

f= fo  =  JlcoJ{0); 

2.  Minimize  augmented  Lagrange  function  (ALF)  A(X,\,rp)  as  an  unconstrained 
function  using  the  conjugate  directions  method  (Procedure  Optimize  ND  CD) 
or  using  the  variable  metric  method  (Procedure  Optimize  NDVM)  as  follows: 

Optimize _ND_CD  (ALCC,  N,  J,  K,  X,  f,  X*,  /,  V„  V,  A), 

or 

Optimize  NDVM  (ALCC,  N,  J,  K,  X,ftX*,/,  V„  V,  A) 

3.  Check  if  the  optimum  solution  has  converged,  using 

|  (f -f)  |  If  <  e;    g/X)  <  0;      |  hk(X')  \   <  8, 
where,  s  and  e,  are  small  numbers. 

4.  Exit,  if  solution  has  converged. 

5.  Update  the  Lagrange  multipliers  using 

Xj  =  \  +  2rp  max[gy(JO,  -X/2rp],  j=\,..J 
h+J  =  Kj  +  2rp  hk(X*),   l^\,..,K 

6.  Update  the  penalty  multiplier  using,    rp  =  yrp. 
If  rp  >  rpmax,  make  rp  =  /•/"". 

7-     /=/ 

8.     Goto  step  2. 
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6.4   Unconstrained  Optimization 

Once  a  constrained  optimization  problem  has  been  converted  into  an 
unconstrained  problem  by  introducing  augmented  Lagrange  function,  it  is  optimized  by 
using 

•  Fletcher  and  Reeves  conjugate  directions  method  [Flet64],  or 

•  variable  metric  methods  —  specifically,   Broydon-Fletcher-Goldfarb-Shanno 
(BFGS)  method  [Huan70]. 

These  two  unconstrained  optimization  methods  were  implemented  and  are  briefly 
described  in  the  following  subsections. 

6.4.1    Conjugate  Directions  Method 

In  this  method,  in  the  initial  search  direction  in  the  first  iteration  is  taken  to  be 
steepest  descent  0=-V/).  On  the  subsequent  iterations  (say  iteration  number  q),  the 
search  direction  is  modified  to  be  a  conjugate  direction  as  follows 

*,  =  -Wy  +  p^  (6.36) 

where  the  scalar  p^  is  defined  as 


p     =      I         q>\  (63?) 

"v     1 12 


IWC-,)I 
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Conceptually  this  method  is  similar  to  Powell's  method,  except  that  here  each  search 

direction  is  conjugate.     Unlike  the  steepest  descent  method,  here  the  information 

gathered  (P)  in  previous  iterations  is  used  to  determine  the  direction  of  search  for  the 
current  iteration.    Theoretically,  this  method  will  minimize  a  quadratic  function  of  N 

variables   in  N     or  fewer  iterations.      In  practice,   it  is  necessary   to     restart  the 

optimization  process  after  every  few  iterations  due  to  the  non-quadratic  nature  of  the 

problems,  or  due  to  numerical  imprecision.    In  order  to  determine  when  to  restart  the 

optimization  process,  two  criteria  are  used.  First,  when  the  one-dimensional  search  fails 

to  improve  the  optimum,  the  process  is  restarted.    Second  when  the  slope  in  a  search 

direction  s  with  parameter  a  is  positive,  i.e. 

dAa)lda  =  V/(A,)-s,  >  0,  (6.38) 

numerical  ill  conditioning  is  indicated  and  sq  is  set  to  -V/A^).  The  conjugate  directions 
algorithm  is  described  in  detail  in  Procedure  Optimize  NDCG. 


Procedure  6.5:    Optimize  NDCG  (function/),  N,  J,  K,  X°,f,  X*,f,  R,  V„  V,  A,  A) 

{ 
Purpose: 

Perform  unconstrained  optimization  of  a  function  ftX)  {X  being  an  N 
dimensional  vector)  using  Fletcher  and  Reeves'  conjugate  gradients  directions 
method. 

Input: 

Function  to  be  optimized  (augmented  Lagrange  function),  /); 
Number  of  decision  variables,  N; 
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Number  of  equality  constraints,  J; 
Number  of  inequality  constraints,  K; 
Initial  values  for  decision  variable  vector,  X°; 
Value  of  the  function  at  X°,  f0. 

Output: 

Optimum  values  for  decision  variable  vector  X"; 

Optimum  value  of  the  function  (augmented  Lagrange  function),  /; 

Optimum  value  of  cost  function,  R; 

List  of  constraint  violations  at  the  beginning  of  life  at  X0,  V;, 

List  of  constraint  violations  at  the  end  of  life  at  A'*,  V; 

Lagrange  multipliers,  X. 

Steps: 

1 .  Evaluate  the  gradient  of  function  X)  using  the  procedure 

Gradient  (function  f(),  N,  X,  f  V/) 
Evaluate  a  =  -V/V/ 

2.  Evaluate  the  search-direction 

s  =  -V/ 

3.  Find  optimum  point  (a*,  X.  and/)  in  the  search  direction  by  1 -dimensional 
golden  sections  search  method  using  the  procedure 

OptimizeJD  (/(),  Optimize _ND,  N,  X,  f  s,  a^  a,  X  /,  R,  V,  V) 

4-    /=/=//=/„; 

5.  If  a*  <  e,  exit. 

6.  Evaluate  the  gradient  of  function  j§  using  the  procedure 

Gradient  (function  f(),  N,  X,  f,  Vf); 

Evaluate  b  =  -V/V/  and  p  =  b/a; 

Evaluate  the  new  search  direction  as,  s  =  -V/+  p\s; 

Assign  a  =  b; 

Evaluate  the  slope  of  the  objective  function  as,  slope  =  sVf. 

7.  If  slope  <  0,  go  to  step  2. 

8.  Find  the  optimum  point  (a*,  X.  and/)  in  the  search  direction  by  1 -dimensional 
golden  sections  search  method  using  the  procedure 
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OptimizeJD  (f(),  Optimize _ND,  N,  X,  f  s,  amax,  a,  X*,  f,  R,  V,  V) 

9.     Check  if  the  optimum  solution  has  converged,  using 

|  (f-f)  |  //<  E,   or,    |  AX  |  /  |  X  |    <  e,   or,  a*  <  e 

where,  s  is  a  small  number. 

Exit,  if  solution  has  converged,  else  goto  step  6. 

} 


6.4.2   Variable  Metric  Method 

The  variable  metric  methods  take  into  account  the  information  gathered  from  the 
previous  iterations,  but  now,  instead  of  carrying  that  information  in  a  single  parameter 
(p),  the  information  is  stored  in  N  parameters.  The  basic  concept  is  to  create  an  array 
which  approximates  the  inverse  of  the  Hessian  matrix  for  quadratic  functions  as  the 
optimization  progresses.  This  method  is  considered  more  reliable  and  efficient 
compared  to  conjugate  directions  method.  Here  the  search  direction  at  iteration  q  is 
defined  as 

sq  =  -HVJ{X)  (6.39) 

In  the  first  iteration,  the  McN  matrix  H  is  taken  to  be  identity  matrix  /,  hence  the  search 
direction  sq  is  simply  the  steepest  descent  direction.  At  the  end  of  iteration  q,  a  new 
H  matrix  is  defined  as 

Hq+I  =Hq  +  Dq  (6.40) 


117 


where  D  is  a  symmetric  update  matrix 


D  =  hppT  +  [HypT  +  piHy)1)/*.  (6.41) 

Vectors  p  and  y  are  defined  as  follows 

P=Xq-  Xq-1 

y  =  VfiXq)  -  VA*W 
and  scalars  a,  x  and  8  are  defined  as  follows 

a=p.y 
r=yTHy 

5  =  (o-+t)/ct2 

The  variable  metric  method  algorithm  is  described  in  detail  in  Procedure  6.6. 


Procedure  6.6:    Optimize _ND_VM  (function /),  N,  J,  K,  X^f,  X,f,  R,  Vt,  V,  A,  X) 

{ 
Purpose: 

Perform    unconstrained    optimization    of  a   function  J{X)    (X  being    an   N 
dimensional   vector)  using  a  variable  metric  method,  or,  Broydon-Fletcher- 
Goldfarb-Shanno  (BFGS)  method. 


Input  and  Output: 

Same  as  in  the  Procedure  Optimize _ND  CG. 

Initializations: 
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J        J  Jp         Jo' 

Steps: 

1.  Initialize  the  pseudo-Hessian  matrix  H[NxN]  to  a  identity  matrix. 
Evaluate  the  gradient  of  function  j\)  using  the  procedure 

Gradient  (function /(9,  N,  X,  f,  V/) 

2.  Evaluate  the  search-direction, 

s  =  -HVf. 

3.  Find  optimum  point  in  the  search  direction  using  1 -dimensional  golden  sections 
search  method  using  the  procedure 

Optimize _1D  (f(),  Optimize _ND,  N,  X,  f  s,  a^  a,  X',  f,  R,  V,  V) 

4.  AX  =  X*  -  X 

5.  Check  if  the  optimum  solution  has  converged,  using 

|(/V)//|    <e;   or,    |  A*  |  /  |  X  \    <e 
where,  s  is  a  small  number.  Exit,  if  solution  has  converged. 

6.  Find  the  new  search-direction  s  using 

Gradient  (function  f(),  N,  X,  /,  V/) 

y  -  v/  -  vf 

c=p.y 
T=yTHy 

(3  =  (a+T)/a2 

Evaluate  the  symmetric  update  matrix  D[N*N]  as  follows 

D  =  fop7  +  [HypT  +  p(Hy)r]/G 
Evaluate  the  pseudo-Hessian  matrix  H  as  follows 

H  =  H  +  D 

7.  Evaluate  the  slope  in  the  search-direction,  5=5.  Vf. 
If  8  >  0,  restart  the  optimization  from  step  1 . 

8.  Update  the  decision  variables  vector  and  the  gradient  vector 
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X  =  X*p  =  X 
v/=v/ 


9.     Goto  step  2. 

} 


6.5   Heuristic  Based  Optimization  Approach 

This  section  deals  with  the  development  of  the  integrated  knowledge  based 
system  for  the  design  of  hydrodynamic  journal  bearings.  The  design  process  is  quasi- 
procedural,  wherein  heuristic  suggested  design  is  validated  using  algorithmic 
procedures.  The  problem  considered  here  assumes  an  a  priori  knowledge  of  the 
operating  specifications,  with  the  aim  of  generating  an  optimum  bearing  configuration. 
The  task  involves  the  preliminary  synthesis  of  a  suitable  bearing  configuration,  dictated 
by  a  decision  making  capability  and  the  user  specified  operating  conditions.  The 
preliminary  configuration  is  analyzed  to  arrive  at  the  optimum  solution  from  a  minimal 
manufacturing  and  life  cycle  operating  costs  standpoint.  The  procedure  is  iterative  in 
nature,  by  which  alternative  configurations  are  synthesized  using  heuristic,  and  the  best 
configuration  is  identified  based  on  economic  criteria  (life  cycle  cost).  A  novel  feature 
of  the  system  is  that  the  design  constraint  information  is  utilized  as  a  feedback  to  the 
heuristic  design  segment  of  the  program  to  improve  the  optimum  design. 

Knowledge   based  methods   now  provide  a  capability   for  representation   of 
qualitative  information,  which  can  be  coupled  with  traditional  procedural  analysis 
methods  to  arrive  at  practical  design  systems.    Hitherto,  such  an  approach  was  not 
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entirely  feasible,  due  to  the  inherent  difficulties  in  the  representation  and  processing  of 

valuable  qualitative  information  that  are  an  essential  element  of  a  successful  design 

practice.    Generic  CAD  systems  available  are  usually  limited  to  problem  solving  tasks 

in  the  algorithmic  or  deterministic  domain.    However,  creative  aspects  in  design  are 

often  ill-structured,  and  require  the  judgement  of  the  designer.  Therefore,  attempts  have 

been  made  in  recent  years  to  incorporate  qualitative  decision  making  capabilities  based 

on  experience  in  design  procedures.    It  is  generally  accepted  that  the  best  approach  is 

one  which  combines  both  heuristic  and  procedural  methods.    This  section  provides  a 

design  methodology  which  integrates  knowledge  based  systems  and  algorithmic  analysis 

methods  for  optimum  life  cycle  design. 

A  schematic  of  the  design  system  is  shown  in  Figure  6.4.    The  system  entails 

two  distinct  modules.     The  first  is  the  heuristic  decision  making  module,  which 

incorporates  an  expert  system  CLIPS  [CLIP88],  and  the  second  is  an  algorithmic 

module  in  'C\  which  includes  the  procedural  analysis  and  optimization. 

6.5.1    Heuristic  Module 

The  basic  framework  for  the  heuristic  decision  system  is  the  CLIPS  rule  based 
environment.  This  utility  invokes  'C  procedural  programs,  making  possible  convenient 
link  between  algorithmic  procedures  and  heuristic.  The  system  provides  for  the 
separation  of  the  knowledge  base,  analysis  procedures  and  domain  independent 
inference  engine.  The  procedure  is  initiated  by  the  input  of  all  the  design  requirements 
to  the  expert  system.     This  system  contains  a  knowledge  base  of  pertinent  design 
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information,  represented  in  the  form  of  rules.    These  rules  are  written  in  the  format 

required  by  the  inference  environment  of  the  expert  system  (CLIPS). 

6.5.1.1   Knowledge  Representation  and  Execution.   The  syntax  of  a  typical  rule 

in  the  knowledge  base  is  shown  below. 


(defrule  <name>  ["comments"] 
(<first  pattern>) 
[(first  pattern>) 


;LHS 


(<n  th  pattern>)] 
=> 

[(<first  action>) 


(<nth  action>)]) 


;RHS 


The  example  in  Figure  6.5  illustrates  a  typical  rule  in  CLIPS  in  the  integrated 
design  decision  system.  The  LHS  consists  of  one  or  more  of  the  patterns  which  form 
the  condition(s)  to  be  satisfied,  before  a  rule  can  be  fired.  An  implicit  'and'  is  present 
in  the  LHS  if  more  than  one  pattern  is  present.  The  RHS  of  the  rule  is  a  list  of  one  or 
more  action(s)  to  be  performed  as  a  result  of  the  firing  of  the  rule.  The  rule  base  also 
includes  meta-rules.  These  are  also  treated  as  rules  in  the  inference  environment,  each 
of  which  is  a  superset  of  two  or  more  rules.  The  advantage  of  using  meta  rules  is  the 
efficiency  of  concatenating  task  specific  rules.  Furthermore,  a  better  control  over 
evaluations  of  rules  is  possible.  In  a  basic  execution  cycle,  the  knowledge  base  is 
examined  to  see  if  the  conditions  of  any  rules  are  satisfied.  This  is  done  by  simple 
matching  of  the  data  with  the  conditions.   If  the  fields  of  the  data  asserted  as  variables, 
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match  the  conditions  of  a  rule,  then  the  rule  is  activated  for  firing.    In  the  case  when 

more  than  one  rule  match,  all  the  matched  rules  are  activated  and  pushed  onto  a  stack. 

The  most  recently  activated  rule  has  the  lowest  priority  and  is  placed  at  the  bottom  of 

the  stack.  The  top  rule  is  selected  and  executed;  as  a  result,  new  rules  may  be  activated 

and  some  previously  activated  rules  may  be  deactivated.    This  cycle  is  recursive  until 

no  further  pattern  matches  are  possible. 


(defrule  rule5  ;name  of  the  rule 

?fl  <-  (parameter  ?p)  ;parameter  variable  p 

?f2  <-  (diameter  ?d)  ;diameter  variable  d 
=> 

(if  (<  ?p  2000.)  ;for  p  less  than  2000 
then 

(bind  ?c  (+  (*  0.00035  ?d)  0.00015)) 

(assert  (clearance  ?c))  ;variable  c  is  a  fact 
else 

(bind  ?c  (+  (*  .0004  ?d)  .0002))       ;c  =  0.0002  +  0.0004  *  d 

(assert  (clearance  ?c))) 

(retract  ?fl  ?f2))  ;remove  facts  fl  and  f2 


Figure  6.5:    Knowledge  representation  in  CLIPS  -  A  sample  rule. 


6.5.1.2  Knowledge  Base  for  Initial  Design.  The  pertinent  information  was 
compiled  from  the  references  used  in  industrial  design  practice  [Full84,  Neal73, 
Tool83].  The  information  supplied  to  the  expert  system  includes  the  bearing  application 
domain,  load,  diameter  of  the  shaft,  and  speed  of  operation.  A  partial  description  of  the 
domain  knowledge  is  included  in  Tables  6.1,  6.2  and  6.3.    Table  6.1  lists  different 
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bearing  configurations  and  corresponding  cost  factors  assigned  to  them.  It  also  lists  the 

factors  for  oil  flow  rates  through  those  bearings.    The  plain  journal  bearings  has  been 

assigned  a  factor  of  1.0  for  both  these  factors.     Table  6.2  lists  different  bearing 

materials  and  associated  material  cost  factors,  harnesses  and  machining  cost  factors. 

These  cost  factors  are  used  to  evaluate  the  costs  in  different  stages  of  the  life-cycle. 

Table  6.3  lists  a  range  of  bearing  applications  and  associated  materials,  tolerances  and 

L/D  ratios  used  in  the  industry  and  as  given  in  the  handbooks.   This  data  in  these  tables 

is  used  in  the  knowledge  base  for  the  selection  of  initial  bearing  design. 

Table  6.1:        Journal  Bearing  Configuration  Data  and  Cost  Factors. 


Configuration 

Configuration 

Relative  Oil 

Cost  Factor  (a) 

Flow  (Qr) 

Axial  Groove  Square  Edge 

1.5 

1.8 

Axial  Groove  Round  Edge 

1.4 

2.0 

Circumferential  Groove  Square  Edge 

1.3 

4.4 

Circumferential  Groove  Round 

Edge 

1.1 

4.9 

Plain 

1.0 

1.0 

Table  6.2:        Journal  Bearing  Material  Data  and  Cost  Factors. 


Material 

Material  Cost 
Factor  ((3) 

Hardness 

(HB) 

Machining  Cost 
Factor  (y) 

Lead-based  Babbitt 

1.0 

21 

1.0 

Tin-based  Babbitt 

1.08 

25 

1.2 

Cadmium  Alloy 

1.45 

35 

1.06 

Copper-Lead 

1.3 

25 

1.02 

Bronze 

0.9 

60 

1.16 

Aluminum  Alloy 

1.2 

45 

1.1 

Silver 

3.0 

35 

1.02 
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Table  6.3:        A  Selected  List  of  Bearing  Applications,  Materials,  Tolerances,  and  L/D 
Ratios. 


Application  Area 

Bearing  Material 

Tolerance 
(in) 

L/D 

Automobile  engines 

Lead  based  babbitt, 

0.0001- 

0.75-1.75 

copper- lead 

0.0007 

Diesel  engines 

Lead-based  babbitt, 

0.0001- 

0.75-2.0 

copper-lead 

0.0007 

Air  compressors 

Lead  based  babbitt 

0.0002- 
0.0004 

1.0-2.0 

Railway  cars 

Cadmium 

0.0005- 
0.0009 

1.8-2.0 

Steam  turbines 

Tin  based  babbitt 

0.0001- 

1.0-2.0 

0.0003 

Generators  &  motors 

Lead  based  babbitt 

0.0009- 
0.0015 

1.0-2.5 

Machine  tools 

Lead  based  babbitt 

0..0001- 

0.0004 

1.5-4.0 

The  task  assigned  to  the  expert  system  includes  the  identification  of  the  initial 
design  parameters.  Specifically,  they  are  the  bearing  material,  configuration  (type), 
length  of  the  bearing,  clearance  between  shaft  and  bearing,  and  lubricant  viscosity.  The 
selection  of  the  material  for  the  bearing  is  chosen  based  on  the  application  domain  of 
the  bearing.  The  chosen  material  has  a  specified  range  of  hardness,  which  defines  the 
material  wear  factor.  The  application  is  categorized  based  on  the  type  of  industrial 
application.  This  categorization  helps  in  the  generation  of  initial  values  for  clearance 
(C),  length  (L),  and  viscosity  of  the  lubricant  (u),  to  be  used  in  the  optimization  phase. 
A  range  of  allowable  clearances,  lengths  and  viscosities  is  available  for  each  application 
area.  A  simple  approach  would  be  to  take  the  average  values  of  the  different  ranges 
as  an  initial  design.    The  initial  bearing  configuration  and  relative  oil  flow  (Qr)  are 
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selected  for  the  given  application  in  accordance  to  tables  given  in  Appendix  A.  The 
configuration  controls  the  relative  oil  flow,  which  is  an  important  component  in 
performance  evaluation.  For  example,  axial  or  circumferential  grooves  and  holes  when 
introduced  in  the  standard  plain  bearing  have  a  significant  effect  on  the  oil  flow.  The 
next  task  is  the  identification  of  the  initial  tolerance  on  the  bearing  diameter.  The 
tolerance  significantly  affects  the  selection  of  the  manufacturing  process  and  bearing 
performance  and  hence,  the  total  cost  of  the  bearing.  The  tolerance  selection  is  based 
on  the  application  domain  and  the  material  selected  previously.  All  the  selected 
parameters  were  introduced  into  the  optimizer  to  achieve  final  design. 

6.5.2   Algorithmic  Module  of  Heuristic  Based  System 

The  design  specifications  which  are  selected  by  employing  heuristic,  are  utilized 
in  this  module  to  attain  an  optimum  design.  The  process  involves  minimization  of  the 
objective  function,  subject  to  the  design  constraints.  Frequent  references  are  made  to 
the  bearing  simulation  module,  to  evaluate  the  bearing  performance  and  operating  cost. 
It  also  makes  references  to  the  heuristic  module  for  evaluating  the  manufacturing  cost. 

Optimization.  The  constrained  minimization  problem  is  converted  to  an 
unconstrained  problem,  using  augmented  Lagrange  multiplier  method  presented  in 
section  sections  6.3  and  6.4.  This  is  solved  employing  a  conjugate  directions  algorithm, 
with  one  dimensional  search  performed  using  golden  sections.  The  objective  function 
evaluation  is  performed  by  invoking  bearing  simulation  module,  and  by  reference  to  the 
heuristic. 
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Bearing  performance  evaluation.     Hydrodynamic   and  thermo-hydrodynamic 

models  of  lubrications  are  utilized  to  evaluate  the  operating  characteristics  of  the  journal 

bearing  such  as  maximum  pressure  (Pmax),  temperature  rise  (AT),  oil  flow  (Q),  minimum 

film  thickness   (h0),  and  coefficient   of  friction   (/).      The   analytical   relationships 

describing  the  journal  bearing  performance  and  constraints  are  given  in  presented  in 

Chapter  2. 

Operating  cost  evaluation.    In  order  to  evaluate  the  operating  cost,  the  bearing 

performance  is  simulated  as  a  function  of  time  as  described  in  Chapter  3.     The 

underlying    assumption    is    that    as    the    bearing    material    wears,    the    clearance 

correspondingly  increases  with  time.    The  time  history  of  the  performance  parameters 

and  total  cost  is  continuously  monitored  and  the  optimum  condition  determined  when 

the  life  cycle  cost/hour  reaches  a  minimum.    That  condition  provides  the  best  solution 

which  can  be  obtained  from  the  initial  design  conditions.    In  other  words,  it  represents 

the  local  optimum  solution  resulting  from  the  output  of  the  heuristic  segment  of  the 

system. 

The  bearing  simulation  is  performed  by  an  incremental  increase  (AC)  in  the 
clearance  due  to  wear.     The  components  of  the  operating  cost  can  be  formulated  as 

sum  of  cost  due  to  energy  loss  and  cost  of  oil  replacements.     As  the  bearing  wear 

progresses,  there  is  a  deterioration  in  its  performance.    Due  to  this,  one  or  more  of  the 

constraints  may  be  violated.    One  of  the  possible  ways  to  rectify  the  problem  is  by 

replacement  of  the  oil  in  the  bearing  which  has  a  cost  associated  with  it.     If  by 

replacing   the   oil,  a  successful    constraint   satisfaction   can   be  achieved,    then  the 
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simulation  is  continued.    Otherwise,  the  bearing  is  presumed  to  have  reached  the  end 

of  its  useful  physical  life.  The  other  criteria  to  end  a  bearing  simulation  is  based  on  its 

economic  life.  This  occurs  when  the  total  life  cycle  cost  per  hour  of  operation  reaches 

a  minimum. 

6.5.3    Feedback  Feature 

A  feedback  feature  is  also  provided  in  the  heuristic  based  system,  which 
facilitates  the  verification  of  changes  in  active  constraints  and  total  life  cycle  cost.  The 
objective  of  this  feature  is  to  provide,  when  possible,  alternative  designs,  which  may 
lead  to  an  improved  optimum  design.  There  are  two  main  parameters  for  which 
alternatives  are  provided  by  the  expert  system,  when  certain  constraints  are  active  or 
violated.  These  parameters  are  bearing  material  and  configuration.  The  heuristic 
module  receives  constraint  information  from  the  algorithmic  module,  and  suggests 
possible  changes.  If  for  example,  the  temperature  constraint  is  active,  one  solution  is 
to  achieve  a  higher  oil  flow,  to  curtail  the  rise  in  temperature.  For  such  an  action, 
heuristic  would  consider  a  different  configuration  or  a  different  bearing  material,  which 
can  withstand  higher  allowable  temperatures.  Out  of  the  two  possibilities,  heuristic 
would  suggest  the  solution  with  a  relatively  lower  manufacturing  cost  factor.  The 
modified  set  of  parameters  are  passed  on  to  the  algorithmic  module  and  the  whole 
process  is  repeated  until  the  life  cycle  cost  for  the  alternative  designs  shows  no  further 
improvement. 
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6.6   Implementation  of  Integrated  Decision  System 

The  dynamic  programming  approach  presented  in  the  last  chapter,  the  multilevel 
optimization  and  heuristic  approaches  presented  in  this  chapter  were  implemented.  The 
algorithmic  programs  were  written  in  C++  language  because  of  its  flexibility  in 
designing  the  data  structure  and  the  programs.  Approximately  6000  lines  of  programs 
were  written.  All  the  relevant  programs  for  dynamic  programming,  constrained  and 
unconstrained  optimization  were  developed,  and  no  commercial  package  was  used. 

The  developed  programs  are  modular  and  designed  to  be  very  general  and 
flexible.  For  example  in  the  dynamic  programming  module,  the  number  of  stages, 
decision  variables  at  each  stage,  number  of  descretizations  at  each  stage  and  the  number 
of  constraints  are  changeable.  When  an  objective  function  or  constraint  need  to  be 
changed,  one  has  to  only  change  the  relevant  modules  which  evaluate  them.  The 
program  modules  which  evaluate  the  objective  function  or  constraints,  are  separate  from 
the  program  modules  which  optimize  these  functions.  Hence  when  the  optimization 
method  need  to  be  changed,  one  only  has  to  modify  the  relevant  modules.  The  number 
of  decision  variables,  constraints  and  objective  functions  are  not  hard  coded  and  can  be 
easily  changed. 

The  most  exacting  and  time  consuming  part  of  the  implementation  was  the 
development  of  the  module  for  simulation  of  the  journal  bearings.  The  development 
was  achieved  after  a  repeated  process  of  trial  and  error  —  where  many  assumptions 
about  the  physics  of  the  bearings  were  made  and  the  results  verified  by  running  the 
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simulations.        Many  times,  the  assumptions   regarding  the  physical   behavior  and 

deterioration  of  the  bearing  performance  were  modified  as  a  result  of  simulations  so  that 

finally  the  bearing  life  simulations  now  bear  a  resemblance   to  the  actual  bearing 

performance  through  its  life.     Literally  thousands  of  such  iterations  of  assumption 

making  were  performed,  leading  to  convergence  into  "sound  assumptions",  which  are 

presented  in  this  thesis. 

6.7   Summary 

The  integrated  design  methodology  presented  in  this  chapter  combines  heuristic 
and  algorithmic  approaches  in  the  integrated  system  for  the  optimum  design  of  machine 
elements  for  minimum  life  cycle  costs.  The  findings  of  the  algorithmic  module  are 
used  by  the  system  in  a  feedback  feature  for  enhancing  the  utility  of  heuristic.  The 
system  is  capable  of  selecting  the  configuration  which  minimizes  the  cost  of  providing 
specified  tasks,  together  with  material  specifications,  design  parameters,  as  well  as 
manufacturing  process  and  tolerances.  Although  a  bearing  design  is  considered  for 
illustration,  the  methodology  offers  sufficient  flexibility  to  be  utilized  in  other  design 
domains. 


CHAPTER  7 
DESIGN  OPTIMIZATION  EXAMPLES 


This  chapter  presents  a  few  examples  of  how  the  integrated  decision  system  is 
used  to  optimize  the  life-cycle  for  journal  bearings.  In  the  first  two  examples  the  same 
bearing,  manufacturing  and  quality  parameters  are  used  and  the  system  performance  is 
based  on  hydrodynamic  and  thermo-hydrodynamic  lubrication  models  respectively. 
Using  the  integrated  design  decision  system,  the  optimum  bearing  designs  are  obtained 
and  the  results  are  then  compared.  In  order  to  compare  the  difference  between  the  two 
analytical  approaches,  the  journal  bearings  optimized  using  hydrodynamic  theory  are 
subsequently  analyzed  using  thermo-hydrodynamic  theory.  An  examples  is  also  given 
for  a  design  optimization  using  the  combined  heuristic  and  algorithmic  procedures. 

7.1    Example  1  —  Hydrodynamic  Model 

In  this  example  the  behavior  of  the  bearing  is  assumed  to  be  governed  by  the 
hydrodynamic  lubrication  theory.  The  task  is  to  design  the  bearing  for  optimum  life- 
cycle  cost  when  the  bearing  length  L,  initial  bearing  clearance  Q,  lubricant  oil  viscosity 
u  and  tolerance  (AD)  are  the  design  decision  variables.  The  integrated  design  system 
is  used  with  the  multi-level  optimization. 
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7.1.1    Design  Specifications 

The  following   bearing  design,   manufacturing    and  quality  control  parameters   are 
specified. 


Design  Parameters 

Bearing  Diameter,  D 

Bearing  Speed,  N 

Load,  W0 

Bearing  material 

Bearing  material  hardness,  Hb 

Journal  material  hardness,  Hs 

Contaminant  hardness,  Hc 

Bearing  wear  coefficient,  Kb 

Journal  wear  coefficient,  Ks 

Rotor  mass,  M 

Oil  inlet  temperature,  T, 

Manufacturing  Parameters 

Cost  factors  a,  P,  y 

Coefficients  of  tolerance  equation 

Coefficients  of  base  cost 


1.75  inch 

3000  rpm 

1000  lb 

Lead-based  Babbitt 

21 

200 

100 

0.016 

0.006 

20  1bm 

110  °F 


1 

ko=0.2,  k,=0.5,  a=l/3 

k2=25,  k3=l 


Quality  Control  Parameters 

In-plant  cost  of  rework/scrapping,  A 
Cost/measurement  of  diameter,  B 
Cost/process  adjustment,  C 
Average  number  of  products  between  successive 

adjustments  interval,  uB 
Quality  control  limits,  ±80 
Process  lag  in  units  of  products,  x 


*-M.base  * 

1.00  $ 
5.0$ 

100  units 
0.0001  in 
5  units 


Operation  and  Support  Parameters 

Oil  cost,  r0 
Energy  cost  rate,  re 
Total  oil  volume,  V0 


0.04  $/in3 
0.01  $/kWhr 
250  in3 
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Sequence  of  Design  Iterations 

Initial  values  for  the  bearing  decision  variables  are  assumed  to  be  as  follows. 

Tolerance  on  bearing  diameter,  AD  0.0004  in, 

Bearing  length,  L  1.75  in, 

Bearing  clearance,  C,  0.001  in, 

Lubricant  oil  viscosity,  u  l.OxlO"6  Reyn, 

Initial  estimate  for  cost  of  external  quality,  CEQ     0.0  $. 

The  multi-level  optimization  algorithms  presented  in  Figure  6.2  and  6.3  are  used  for 

optimization.    The  sequence  of  design  iterations  performed  is  described  below  and  is 

also  presented  in  Table  7.1. 

7.1.2  Iteration  1 

Once  the  value  of  input  parameters  has  been  specified,  the  internal  quality  cost 
CIQ  is  evaluated.  At  level  one  optimization,  the  decision  system  optimizes  the  life-cycle 
cost  rate  (LCC  divided  by  bearing  life)  with  bearing  length  (L),  initial  bearing  clearance 
at  the  beginning  of  life  (CJ,  and  lubricant  viscosity  (n)  treated  as  decision  variables. 
Life-cycle  cost  (LCC)  is  defined  as 


CIC      =  manufacturing  cost  +  operation  cost 

+  internal  quality  cost  +  external  quality  cost 

CM(L,  D,  AD)  +  C0(L,  D,  q,  u)  +  CIQ(AD)  +  CEQ(AD,  CQn(C&.l) 


Bearing  life  is  obtained  by  simulating  the  bearing  performance  over  time,  and 
then  choosing  the  minimum  of  its  physical  and  economic  life.  An  example  of  such  a 
simulation  and  resulting  cost  components  versus  the  operation  time  of  the  bearing  are 
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showed  in  Figure  7.1.  The  optimization  is  carried  out  subject  to  constraints  on 
temperature,  maximum  pressure,  minimum  film  thickness,  contaminant  to  oil  ratio, 
bearing  length  to  diameter  ratio,  bearing  clearance,  lubricant  viscosity,  minimum  life 
of  the  bearing  and  dynamic  stability  of  the  bearing.  These  constraints  can  be  described 
as  follows: 
minimum  oil  film  thickness: 

K  >  homin  =  5.0xl0-5  in, 
maximum  outlet  temperature: 

T0  <  Tmax  =  300  T, 
maximum  pressure: 

P—  *  *U«  -  30xl03  psi, 
lower  limit  on  the  viscosity  of  the  oil  used: 

u  >  umin  =  O.OlxlO"6  Reyn, 
upper  limit  on  the  viscosity  of  the  oil  used: 

u  <  umax  =  50.0x1 0-6  Reyn, 
upper  limit  on  length  to  diameter  ratio: 

L/D<{L/D)max=\.0, 
lower  limit  on  length  to  diameter  ratio: 

L/D  >  (L/D)min  =  0.25, 
and  the  criterion  ensuring  the  dynamic  stability  of  the  bearing  described  by  equation 
2.8. 
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Figure  7.1:     Life-cycle  cost  components  versus  time  of  operation  obtained  by  bearing 
simulation. 
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Optimum    values    of  the    decision    variables    obtained    after   the    first   level 

optimization  are  1=1.67  in,  C^O.0009  in,  u=2.151xl0"7  Reyn.   Other  important  results 

for  the  optimum  design  are  bearing  life=10800  hours,  manufacturing  cost  CM=208.8  $, 

cost  of  energy  loss  COe=303.25  $,  COo=630  $,  and  internal  quality  cost  CEQ=3.4  $. 

Figure  7.2  displays  the  total  cost  and  cost  components  as  a  function  of  initial  clearance. 

The  operation  cost  as  a  function  of  initial  bearing  clearance  Q  around  its  optimum 

value  (Ci*=0.0009  in)  is  automatically  obtained,  and  is  used  to  obtain  the  2nd  derivative 

of  operation  cost  with  respect  to  initial  clearance  at  its  optimum  value  as  follows 

Q>"(Qlo-c/«    =  5.842xl010  $/in2 
Now  all  the  four  components  of  life  cycle  cost  can  be  expressed  in  terms  of  tolerance 
AD  in  the  decision  system  as  in  equation  (7.1).   Figure  7.3  shows  the  life-cycle  cost  as 
a  function  of  tolerance  AD.    Thus  the  optimum  value  of  tolerance  AD  is  obtained  by 
using  a  directional  search  method.    In  the  example,  it  produces  tolerance 

AD        =  0.0001  inch, 
external  quality  cost 

CEQ       =  64.91  $, 
total  life-cycle  cost 

CLC      =  1206.96  $, 
and  cost  rate 

cr         =  QfVlife  =  0.1117  $/hour. 
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Figure  7.3:     Life-cycle  cost  components  versus  bearing  tolerance  on  diameter  (AD). 
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In  this  iteration,  the  first  level  optimization  was  carried  out  to  find  the  optimum 

decision  variables  L,  C  and  u  with  the  tolerance  AD  treated  as  a  fixed  parameter.   Then 

the  second  level  optimization  was  done  to  find  the  optimum  tolerance  ^  with  L,  C,  and 

\i  treated  as  fixed.     Iteration  1  is  based  on  incomplete  cost  information,  since  the 
accurate  values  of  internal  and  external  quality  costs  is  not  available  at  this  stage.   The 

output  of  this  iteration  is  used  to  compute  the  accurate  costs  in  the  2nd  iteration.   Now, 

with  these  intermediate  optimum  values  of  the  parameters  and  the  tolerance,  the  2nd 

iteration  can  be  started. 

7.1.3   Iteration  2 

For  iteration  2,  the  initial  values  of  design  variables  as  a  result  of  iteration  1  are 
1=1.67  in,  C=0.0009  in.,  u=2.151xl0'7  Reyn,  and  AD=0.0001  in.  With  these  initial 
values,  the  first  level  optimization  produces  L=1.67  in,  C=0.00081  in,  u=  2.016xl0"7 
Reyn,  life=11500  hours,  CM=  320.61  $,  energy  cost  C0  =319  $,  oil  cost  COo=760  $, 
C/e=30.71  $.  As  in  iteration  1,  operation  cost  curve  as  a  function  of  initial  clearance 
C,  (a  curve  similar  to  that  in  Figure  7.2)  is  used  to  obtain  the  double  derivative  of 
operation  cost  as  follows 

Co'XC,')  =       5.2295x10'°      at  C>0.00081  in 

AD   =  0.0001  in 

As  in  iteration  1,  a  directional  search  method  is  used  to  obtain  the  optimum  tolerance 
from  costs  versus  AD  curve  (a  curve  similar  to  that  in  Figure  7.3).    It  produces 

CEQ   =  58.11, 


140 
CIC  =  1488.43, 

Cost  rate  -     cT  -  (C£c/r»/c)  =  0.1294  $/hr. 

7.1.4   Optimum  Results 

In  each  iteration,  the  decision  system  automatically  goes  through  a  similar  cycle 

of  the  level-one  and  level-two  optimization  until  the  decision  variables  converge. 

Although,  this  example  has  been  described  using  plots  of  various  quantities,   the 

decision  system  is  completely  automated  and  the  optimization  process  does  not  require 

any  reference  to  these  plots.     But  the  decision  system  can  generate  these  plots  if 

instructed  by  a  user  to  do  so.    In  this  example,  by  the  end  of  the  4th  iteration,  the 

decision  variables  and  objective  function  converge,  hence  the  computation  process  is 

stopped.  Final  values  of  important  life-cycle  parameters  after  the  design  has  converged, 

are  given  in  Table  7.1.  The  optimum  values  for  quality  control  parameters  obtained  are 

as  follows 

Optimum  process  measurement  interval,  m  3.64  «  4  units 

Optimum  quality  control  limits,  8  0.000025  in 

Predicted  average  process  adjustment  interval,  u  9.94  «  10  units 

Optimum  process  capability  index,  Cp  1.23 

The  plots  of  various  physical  quantities  for  the  optimum  bearing,  as  it  is 

simulated  over  its  life,  are  given  in  Figures  7.4  to  7.10.    Only  the  first  260  hours  of 

simulation  is  used  here  for  plotting  so  that  the  graphs  are  clear  and  intelligible.    When 

the  oil  is  replaced  in  the  bearing,  most  physical  parameters  (such  as  oil  viscosity, 

pressure,  temperature,  contaminants,  etc.)  change  suddenly.    Thus  the  discontinuities 
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(kinks)  in  these  simulation  plots  indicate  the  time  when  the  oil  is  replaced  in  the 

bearing.    For  example,  in  Figures  7.4  -  7.10,  the  oil  is  replaced  at  100  and  210  hours. 

Figure  7.4  shows  the  plots  of  the  Sommerfeld  number  (5)  and  the  modified 
Sommerfeld  number  during  the  bearing  life  simulation.  Figure  7.5  shows  the  plots  of 
the  wear  factor  (WF)  and  oil  flow  rate  (Q)  during  simulation.  The  wear  factor  in  the 
bearing  is  given  by  equation  (3.11a)  and  it  increases  with  time  due  to  increased 
contaminants,  oil  temperature,  and  pressure  until  the  time  the  oil  is  replaced.  The 
bearing  oil  flow  decreases  with  time. 

Dimensionless  rotor  mass  {<km=MC<n2IW)  and  stability  criterion  (v|/(a),  a  function 

of  modified  Sommerfeld  number  a),  which  are  the  quantities  used  to  determine  the 
dynamic  instabilities  in  the  bearing,  are  plotted  in  Figure  7.6.    With  time,  the  bearing 

clearance  (Q  increases  due  to  wear  thus  increasing  the  dimensionless  rotor  mass.   For 

clarity,  the  two  quantities  are  plotted  on  a  different  scale. 

Figure  7.7  shows  the  plot  of  average  viscosity  (u.)  and  average  temperature  (Tavg) 

of  the  oil  film  in  the  bearing.    Because  of  the  contaminants  and  sludge  generated  over 

time,  the  viscosity  rises  until  the  oil  is  replaced.    The  temperature  rise  takes  place 

because  of  the  increased  viscosity  and  reduced  oil  flow  rate  through  the  bearing.  A  plot 

of  maximum  pressure  (P^)  in  the  oil  film  is  shown  in  Figure  7.8.    The  maximum 

pressure  decreases  until  the  oil  replacement,  because  the  pressure  distribution  across  the 

oil  film  becomes  flatter.    Figure  7.9  shows  plots  of  minimum  film  thickness  (h0)  and 

coefficient  of  friction  if).  Until  oil  replacement,  both  the  quantities  rise  because  of  the 

increasing  oil  viscosity.   Figure  7.10  shows  the  plots  of  bearing  clearance  (Q  and  ratio 
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of  contaminants  to  oil  (/•„,).    Due  to  bearing  and  journal  wear,  the  clearance  increases 

continuously  and  the  contaminants  increase  until  the  oil  is  replaced. 

7.1.5   Exploration  of  Design  Space 

Figure  7.1  is  a  plot  of  various  cost  simulations  as  the  bearing  goes  through  its 
life.  As  expected,  The  total  cost  (LCQ  keeps  rising  over  the  life.  More  frequent  oil 
changes  at  the  end  of  bearing  life  increase  the  slope  of  oil  cost.  As  expected,  the  cost- 
rate  (LCC  divided  by  the  number  of  hours  of  operation)  keeps  decreasing.  In  this 
particular  example,  the  simulation  is  terminated  because  the  bearing  has  reached  the  end 
of  its  economic  life  where  the  cost-rate  curve  starts  rising. 

Figure  7.2  shows  the  plots  of  various  costs,  cost-rate,  Cfinal  (final  bearing 
clearance  at  the  end  of  life),  and  bearing  life  for  optimum  designs  as  a  function  of  Cinitia| 
(bearing  clearance  at  the  beginning  of  life).  Note  that  the  final  clearance  has  almost  a 
constant  value  of  0.0019  inch  where  the  bearing  approaches  the  dynamic  instabilities. 
Life  keeps  decreasing  as  the  initial  clearance  increases  since  the  increasing  clearance 
reduces  the  allowance  for  bearing  wear.  Similarly,  the  oil  cost  and  energy  cost  decrease 
with  increasing  Cinitial.  It  is  interesting  to  note  that  the  cost-rate  curve  has  a  clear 
minimum  at  Cinitiai=.  00081  approximately,  hence  it  defines  the  optimum  design  clearance 
for  the  bearing. 

Figure  7.3  plots  various  costs  of  optimum  designs  as  a  function  of  bearing 
tolerance  on  diameter  AD.  Again,  as  expected,  the  external  quality  cost  (CEQ)  increases 
with  increasing  tolerance,  since  the  loose  tolerance  causes  the  bearing  performance  to 
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deviate  from  its  optimal  value.   The  cost  of  internal  quality  (CIQ)  drastically  decreases 

with  increasing  tolerances.    This  is  because  the  loose  tolerances  require  less  frequent 

quality  control  inspection,  process  adjustments  and  manufacturing  cost.    The  quantity 

of  interest  here  is  the  total  cost  (manufacturing  cost  +  internal  quality  cost  +  external 

quality  cost)  which  has  a  minimum  at  AD=0001  inch,  confirming  the  hypothesis  that 

there  exists  an  optimum  tolerance. 


Table  7.1:        Results  of  the  sequence  of  iterations  for  parameter  and  manufacturing 
tolerance  optimization  for  example  1 . 


Iterations 

Initial  Values 
Chosen 

Iteration  1 
Results 

Iteration  2 

Results 
(optimum) 

Parameters 

Length,  L  (in.) 

Clearance,  C,  (in.) 

Oil  Viscosity,  \a  (Reyn) 

Tolerance,  AD  (in.) 

Mfg.  Cost,  CM  ($) 

Internal  Quality  Cost,  CIQ  ($) 

External  Quality  Cost,  CEQ  ($) 

Energy  Cost,  C0e  ($) 

Oil  Cost,  C0o  ($) 

C0"(Q,  ($/in2) 

Life  Cycle  Cost,  CIC  ($) 

Life  (hours) 

Cost-Rate,  cr  ($/hour) 

1.75 

0.001 

l.OxlO'6 

0.0004 

0.00 

1.67 
0.0009 

2.15xl0"7 
0.0001 
208.80 
3.40 
64.91 
303.25 
630.00 

5.842xl010 

1206.96 

10800 

0.1117 

1.67 
0.00081 

2.01  xlO"7 
0.0001 
320.61 
30.71 
58.11 
319.00 
760.00 

5.23x10'° 

1488.43 

11500 

0.1294 
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7.2   Example  2  —  Thermo-hydrodynamic  Model 

In  this  example,  thermo-hydrodynamic  theory  is  used  to  model  the  journal 
bearing  with  same  design,  manufacturing  and  quality  control  parameters  as  in  example 
1 .  The  major  differences  between  the  hydrodynamic  lubrication  model  and  the  thermo- 
hydrodynamic  lubrication  model  are  explained  in  Chapter  3. 

In  this  section  the  design  optimization  is  carried  out  for  a  bearing  using  two 
models.  The  optimum  bearing  design  problem  of  example  1  and  2  is  repeated  here  for 
three  different  loads  (250,  1000,  2000  lb)  and  4  different  speeds  (1000,  3000,  5000, 
7000  rpm).  Thus,  twelve  different  combinations  of  loads  and  speeds  are  used  as  design 
specifications  for  design  optimization  using  both  hydrodynamic  and  thermo- 
hydrodynamic  models.  The  result  of  design  optimization  for  these  two  bearing  models 
are  then  compared.  Since  we  already  know  that  the  thermo-hydrodynamic  model  is 
more  accurate,  an  interesting  question  is  that  how  the  optimum  bearings  designed  with 
hydrodynamic  models  actually  perform  compared  to  the  optimum  bearings  designed 
with  thermo-hydrodynamic  models.  In  other  words,  the  objective  here  is  to  find  the 
improvement  in  bearing  performance  and  costs  if  the  bearing  is  designed  using  the 
accurate  design  theory  compared  to  the  bearing  designed  with  inaccurate  design  theory. 
This  also  yields  the  value  (benefit)  of  right  information/theory.  Hence  the  bearings 
already  optimized  using  hydrodynamic  models  are  simulated  using  the  thermo- 
hydrodynamic  model.    The  comparison  of  these  three  cases  now  follows. 
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7.2.1    Optimum  Cost-Rates 

Table  7.2  and  7.3  show  that  the  optimum  cost-rate  for  both  hydrodynamic  and 
thermo-hydrodynamic  bearings  is  very  sensitive  to  bearing  speeds  and  loads.  With 
increasing  bearing  speed  and  load,  higher  pressures  and  temperatures  are  generated, 
resulting  in  increased  energy  loss  and  wear,  thus  increasing  the  cost-rates  of  the 
optimum  designs.  The  bearings  optimized  using  thermo-hydrodynamic  model  generally 
yield  higher  cost-rates. 

The  bearings  designed  with  hydrodynamic  model  are  also  simulated  as  thermo- 
hydrodynamic  bearings  and  the  new  cost-rates  are  given  in  Table  7.4.  Comparing 
Tables  7.2  and  7.3  reveals  that  actual  cost-rates  in  a  bearing  optimized  using 
hydrodynamic  theory  are  much  higher  than  those  predicted  by  hydrodynamic  theory. 
These  optimum  cost  rates  as  a  function  of  bearing  speed  and  load  for  both 
hydrodynamic  and  thermo-hydrodynamic  model  are  plotted  in  Figure  7.11.  Sometimes 
the  bearings  optimized  using  hydrodynamic  model  do  not  even  actually  satisfy  the 
design  constrains.  When  these  bearings  are  simulated  using  the  thermo-hydrodynamic 
model,  their  life  turns  out  to  be  very  short  and  hence  cost-rates  are  very  high.  In  Figure 
7.12,  the  cost  rates  are  compared  for  (a)  bearings  optimized  using  hydrodynamic  model, 
and,  (b)  thermo-hydrodynamic  simulation  of  bearings  optimized  using  hydrodynamic 
model.  The  thermo-hydrodynamic  simulation  of  optimum  bearins  sometimes  results  in 
very  short  life  and  consequently  very  high  cost  rates  —  this  is  evident  in  three  cases 
(250  lb  @  7000  rpm),  (2000  lb  @  1000  rpm),  and  (2000  lb  @  7000  rpm)  in  Figure 
7.12.    In  Figure  7.13  the  plot  of  normalized  cost-rates  as  a  function  of  bearing  speed 
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is  shown.    Here  the  cost-rate  is  normalized  using  the  cost-rate  for  bearing  optimized 

using  the  thermo-hydrodynamic  model,  i.e. 

actual  cost-rate 
normalized  cost-rate    =   


cost-rate  for  the  bearing  optimized 
using  thermo-hydrodynamic  model 


It  is  clear  from  Figures  7.12  and  7.13  that  the  actual  cost-rates  for  the  bearing  optimized 
using  hydrodynamic  bearings  are  much  higher,  specially  at  higher  speeds. 


Table  7.2:       Optimum     Cost-Rates     (10 2    $/hr)    —    Bearings     optimized     using 
hydrodynamic  lubrication  model. 


Load  (lb) 

250 

1000 

2000 

Speed  (rpm) 

1000 

1.44 

3.73 

4.50 

3000 

2.74 

4.63 

7.86 

5000 

4.64 

6.92 

10.39 

7000 

7.57 

10.95 

15.26 

Table  7.3:        Optimum  Cost-Rates  (102  $/hr)  —  Bearings  optimized  using  thermo- 
hydrodynamic  lubrication  model. 


Load  (lb) 

250 

1000 

2000 

Speed  (rpm) 

1000 

1.96 

4.09 

7.01 

3000 

2.79 

5.80 

9.56 

5000 

4.95 

7.79 

11.00 

7000 

12.51 

9.89 

13.93 
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Table  7.4:        Cost- Rates  (102  $/hr)  —  thermo-hydrodynamic   simulation  of  bearing 
designed  for  optimum  hydrodynamic  model. 


Load  (lb) 

2000 

250 

1000 

Speed  (rpm) 

1000 

1.60 

6.30 

00 

3000 

3.30 

10.00 

14.3 

5000 

12.86 

63.78 

157.54 

7000 

00 

14.40 

00 

7.2.2   Optimum  Life 

Table  7.5  shows  that  all  the  optimum  hydrodynamic  bearings  achieve  the 
designed  maximum-life  of  25000  hours.  Except  in  one  case  (250  lb  load,  7000  rpm 
speed),  all  the  thermo-hydrodynamic  bearings  also  achieve  the  same  designed 
maximum-life  of  25000  hours.  When  the  bearings  optimized  using  hydrodynamic 
lubrication  theory  are  simulated  as  thermo-hydrodynamic  bearings  (Table  7.7),  in  most 
cases  the  life  of  the  bearings  turns  out  to  be  much  shorter.  A  life  of  zero  in  Table  7.7 
indicates  that  the  initial  bearing  design  does  not  even  satisfy  some  of  the  design 
constraints.  Figure  7.14  shows  the  normalized  life  of  bearings  for  the  two  models.  It 
is  clear  from  these  figures,  that  specially  at  high  speeds,  the  actual  bearing  life  is  much 
lower  compared  to  what  is  predicted  by  optimum  solution  based  on  hydrodynamic 
theory. 
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Cost  Rate    2 


1000 


3000  5000 

Speed,  N  (rpm) 


7000 


'Optimum  Ihdrodynamic  Bearing 


-•"THD  Simulation  of  Optimum  HD 
Bearing 


_L 


Normalized  Cost  Rate  for  1000  lb  Load 


Cost  Rate 


1000 


3000  5000 

Speed,  N  (rpm) 


7000 


Optimum  Hydrodvnnmic  Bearing 


■   THD  Simulation  or  Optimum  HO 
Bearing 


Figure  7.12:   Normalized  cost  rates  versus  bearing  speed  for  different  bearing  loads. 
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7.2.3    Optimum  Manufacturing  Tolerances 

Figure  7.15  shows  the  plots  of  optimum  tolerances  as  a  function  of  bearing 

speed  (JV)  for  three  different  bearing  loads  (250,  1000  and  3000  lb)  for  the  two  models. 

Table   7.8   and   7.9   list  the   values   of  optimum   tolerances.      Note   that   in  both 

hydrodynamic   and  thermo-hydrodynamic    models,  the  optimum   bearing  tolerances 

generally  decrease  with  increasing  bearing  speeds.    The  higher  bearing  speeds  require 

more  precision,  because  at  higher  speeds,  deviations  from  designed  clearance  affect  the 

bearing   pressures   and  temperatures   more  adversely   than  at  lower  speeds.      The 

conditions  of  low  bearing  load  and  high  speed  (for  example  250  lb,  7000  rpm)  tend 

towards  dynamic  instabilities,  as  modeled  by  Lund  and  Saibel  instability  criterion. 

Hence  low  loads  and  high  speeds  require  low  bearing  clearance  (Q  and  tight  tolerance. 

The  conditions  of  high  load  and  low  speed  (for  example  2000  lb,  1000  rpm)  also 

require  tight  tolerances.    Bearings  optimized  using  thermo-hydrodynamic  model  also 

requires  tighter  tolerances  compared  to  the  hydrodynamic  model. 


Table  7.5:        Life  (hrs)  -  optimum  hydrodynamic  lubrication. 


Load  (lb) 

250 

1000 

2000 

Speed  (rpm) 

1000 

25000 

25000 

25000 

3000 

25000 

25000 

25000 

5000 

25000 

25000 

25000 

7000 

25000 

25000 

25000 

Table  7.6:        Life  (hrs)  -  optimum  thermo-hydrodynamic  lubrication. 
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Load  (lb) 

250 

1000 

2000 

Speed  (rpm) 

1000 

25000 

25000 

25000 

3000 

25000 

25000 

25000 

5000 

25000 

25000 

25000 

7000 

7100 

25000 

25000 

Table  7.7:       Life  (hrs)  -  thermo-hydrodynamic  simulation  of  bearing  designed  for 
optimum  hydrodynamic  model.  (7.3c) 


Load  (lb) 

2000 

250 

1000 

Speed  (rpm) 

1000 

25000 

16000 

0 

3000 

25000 

6500 

6750 

5000 

3200 

370 

210 

7000 

0 

5000 

0 

Table  7.8:        Tolerance  (±,  xlO'3  inch)  -  optimum  hydrodynamic  lubrication. 


Load  (lb) 

250 

1000 

2000 

Speed  (rpm) 

1000 

1.24 

0.64 

0.40 

3000 

0.46 

0.65 

0.46 

5000 

0.40 

0.62 

0.52 

7000 

0.30 

0.43 

0.50 
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Normalized  Life  2 


Normalized  Life  Based  on  Hydrodyanamic  Theory 


1000 


3000  5000 

Speed,  N  (rpm) 


50  1b 


1000  lb 
2000  1b 


7000 


250  1b 


1000  lb 


2000  lb 


Normalized  Life  Based  on  Thermodynamic  Simulation  of 
Optimum  Hydrodynamic  Bearing 


0.75 


Normalized  Life 
(hrs) 


0.25 


1000 


3000  5000 

Speed,  N  (rpm) 


■  250  lb    -°- 1 000  lb  -»-2000  lb 


1000  Ifc 


50  1b 


7000 


Figure  7.14:   Normalized  life  versus  bearing  speed  at  different  bearing  loads. 
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Tolerance,   (in) 


1000 


3000  5000 

Speed,  N  (rpm) 


0  250  lb      ■    1000  lb    ®  3000  lb 


Manufacturing  Tolerance  for  Thermohydrodynamic  Theory 


0.0014 
0.0012 


Tolerance  (in) 


1000 


3000  5000 

Speed,  N  (rpm) 


7000 


0  250  1b 


1000  lb    »  2000  lb 


Figure  7.15:   Optimum  manufacturing  tolerances  versus  bearing  speed  for  different 
loads. 


162 


Table  7.9:        Tolerance  (±,  xlO'3  inch)  -  optimum  thermo-hydrodynamic  lubrication. 


Load  (lb) 

250 

1000 

2000 

Speed  (rpm) 

1000 

0.70 

0.58 

0.28 

3000 

0.62 

0.44 

0.25 

5000 

0.50 

0.41 

0.50 

7000 

0.22 

0.43 

0.55 

7.2.4   Optimum  Design  Parameters 

As  a  result  of  optimization  for  the  12  combinations  of  speeds  and  loads  using 
both  hydrodynamic  and  thermo-hydrodynamic  models,  the  optimum  values  of  design 
parameters  were  obtained.  The  optimum  length  to  diameter  ratio  (L/D)  for 
hydrodynamic  and  thermo-hydrodynamic  models  is  plotted  in  Figures  7.16  and  7.17 
respectively.  The  optimum  initial  clearance  (C,)  for  hydrodynamic  and  thermo- 
hydrodynamic  models  is  plotted  in  Figures  7.18  and  7.19  respectively.  The  optimum 
relative  oil  viscosity  (G^  for  hydrodynamic  and  thermo-hydrodynamic  models  is  plotted 
in  Figures  7.20  and  7.21  respectively. 

7.3    Heuristic  Example 


The  heuristic  based  integrated  design  system  described  in  section  6.5  was  used 
for  solving  a  sample  design  problem  in  this  section.  The  task  assigned  is  the  design  of 
a  suitable  bearing  for  the  following  specifications. 
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7.3.1    Heuristic  Module  Input 

Application  Compressor  main  bearing 

Speed  (rpm)  2000 

Diameter  (in.)  1.75 

Load  (lb)  1750 

The  output  from  the  various  modules  are  given  below. 


7.3.2   First  Iteration  Through  Heuristic 

The  heuristic  suggested  the  following  initial  design  for  the  bearing. 
Material  Lead-based  Babbitt 

Hardness  21 

L/D  Ratio  1.0 

Clearance  (in.)  0.001375 

Tolerance  (in.)  0.0003 

Viscosity  (reyn)  3.3x1 0"6 

Configuration  Plain  bearing 


Based  on  the  initial  design,  the  optimization  module  provided  an  optimum  design 
with  11010  hours  of  life.  At  the  end  of  life  of  the  optimum  bearing,  the  wear 
progresses  to  a  point  where  minimum  film  thickness  constraint  and  maximum  pressure 
constraint  are  violated.    The  parameters  for  optimum  design  are  as  follows. 


Life  (hrs) 

11010 

Iterations 

5 

Length  (in.) 

1.5556 

Clearance  (in.) 

0.001 

Viscosity  (Reyns) 

0.6438x1 0"6 

Oil  replacement  interval 

110  hrs 

Manufacturing  Cost  ($) 

200.03 

Oil  Cost  ($) 

543.75 

Energy  Cost  ($) 

175.32 

Operations  Cost  ($) 

719.07 

Life  Cycle  cost  ($) 

919.10 

Life  Cycle  Cost  per  hour 

of  operation  ($/hr)      0.0835 
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7.3.3    Second  Iteration  Through  Heuristic 

The  active  constraint  information  was  passed  from  the  algorithmic  module  to  the 
heuristic  module.  Since,  The  bearing  designed  above  has  excessive  wear;  the  heuristic 
suggested  tin-based  babbitt  which  is  harder  material  (Hb=35). 

Based  on  the  modified  design,  the  optimization  module  provided  an  optimum 
design  with  14160  hours  of  life.  As  in  the  first  iteration,  the  optimum  bearing  violates 
the  minimum  film  thickness  and  maximum  pressure  constraints  at  the  end  of  its  life. 
The  parameters  for  the  second  optimum  design  are  as  follows. 


Life  (hrs) 

14160 

Iterations 

5 

Length  (in.) 

1.5556 

Clearance  (in.) 

0.001 

Oil  Viscosity  (reyns) 

0.6438X10"6 

Oil  replacement  interval 

110  (hrs) 

Manufacturing  Cost  ($) 

220.32 

Oil  Cost  ($) 

693.75 

Energy  Cost  ($) 

225.41 

Operations  Cost  ($) 

919.16 

Life  Cycle  cost  ($) 

1139.48 

Life  Cycle  Cost  per  hour 

of  operation  ($/hr) 

0.0805 

7.3.4   Third  Iteration  Through  Heuristic 

In  order  to  further  reduce  the  wear  rate,  heuristic  suggested  bronze  which  is  the 
hardest  material  (Hb=60)  in  the  knowledge  base.  Interestingly,  this  time  the  algorithmic 
module  produced  an  optimum  design  with  more  than  50000  hours  of  life.  Up  to  the 
end  of  simulation,  no  constraint  was  violated.    Since  this  design  has  the  minimum  life 
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cycle  cost  per  hour  of  operation,  with  no  active  constraint,  this  is  declared  as  the  final 

design.    Parameters  for  the  final  design  are  as  follows. 


Life  (hrs) 

50000 

Active  Constraints 

none 

Iterations 

8 

Length  (in.) 

1.5556 

Clearance  (in.) 

0.001 

Viscosity  (reyns) 

0.6438x1 0"6 

Oil  replacement  interval 

110  (hrs) 

Manufacturing  Cost  ($) 

290.04 

Oil  Cost  ($) 

2449.60 

Energy  Cost  ($) 

795.70 

Operations  Cost  ($) 

3254.30 

Life  Cycle  cost  ($) 

3535.34 

Life  Cycle  Cost  per  hour 

of  operation  ($/hr), 

0.0707 

It  is  interesting  to  note  that  the  feedback  feature  into  the  heuristic  module 
modified  the  input  to  the  algorithmic  module  and  consequently  produced  a  final  design 
which  is  unconstrained,  and  has  a  lower  life  cycle  cost  per  hour  than  the  previously 
obtained  local  optimum.  The  material  selected  due  to  first  feedback  is  tin  based  babbitt, 
which  has  a  20%  higher  hardness  and  is  10%  more  expensive  in  comparison  with  lead 
based  babbitt.  The  material  selected  due  to  second  feedback  is  bronze,  which  has  a 
114%  higher  hardness  and  is  25%  more  expensive  in  comparison  with  lead  based 
babbitt.  Hence  finally  bronze  is  selected  since  it  gives  the  design  with  the  minimum 
life  cycle  cost  per  hour  of  operation  and  with  no  active  constraints. 
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Optimum  Results 

The  output  from  the  two  modules  during  the  design  process  are  given  in  Table 

7.10.     It  is  interesting  to  note  that  the  feedback  feature  into  the  heuristic  module 

modified  the  input  to  the  algorithmic  module  and  consequently  produced  a  final  design 

which  is  unconstrained,  and  has  a  lower  life  cycle  cost  per  hour  than  the  previously 

obtained  local  optimum.  The  material  selected  due  to  first  feedback  is  tin  based  babbitt, 

which  has  a  20%  higher  hardness  and  is  1 0%  more  expensive  in  comparison  with  lead 

based  babbitt.    The  material  selected  due  to  second  feedback  is  bronze,  which  has  a 

114%  higher  hardness  and  is  25%  more  expensive  in  comparison  with  lead  based 

babbitt.    Hence  finally  bronze  is  selected  since  it  gives  the  design  with  the  minimum 

life  cycle  cost  per  hour  of  operation  and  with  no  active  constraints. 


Table  7.10:      Design  iterations  for  the  heuristic  based  optimization. 
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Bearing  Parameters 

Iteration  1 

Iteration  1 

Iteration  3 

Heuristic  Module  Output: 

Material 

Lead-based 

Tin-based 

Bronze 

babbitt 

babbitt 

Hardness  (BHN) 

21 

35 

60 

L/D  ratio 

1.0 

- 

- 

Clearance  (in.) 

0.001375 

- 

- 

Tolerance  (in.) 

0.0003 

- 

- 

Oil  Viscosity  (Reyn) 

3.3xl06 

- 

- 

Configuration 

plain 

plain 

plain 

bearing 

bearing 

bearing 

Algorithmic  Module  Output: 

Iterations  for  Optimization 

5 

7 

8 

Life  (hrs) 

11010 

14160 

>50000 

Length  (in.) 

1.5556 

1.5556 

1.5556 

Clearance  (in.) 

0.001 

0.001 

0.001 

Viscosity  (Reyn) 

0.6438x1 0"6 

0.6438x1 0"6 

0.643  8xl0"6 

Oil  Replacement  Interval 

110 

110 

110 

(hrs) 

Manufacturing  Cost  ($) 

200.03 

220.32 

290.04 

Oil  Cost  ($) 

543.75 

693.75 

2449.60 

Energy  Cost  ($) 

175.32 

225.41 

795.07 

Operation  Cost  ($) 

719.07 

919.16 

3254.30 

Life  Cycle  Cost  ($) 

919.10 

1139.48 

3535.34 

Life  Cycle  Cost  per  Hour  of 

0.0835 

0.0805 

0.0707 

Operation  ($/hr) 

CHAPTER  8 

MODELS  FOR  COST/BENEFIT  ANALYSIS  OF  R&D 

AND  NEW  TECHNOLOGIES 


The  benefit  of  engineering  or  scientific  research  are  many,  some  of  which  are 
not  measurable  in  purely  economic  terms.  However  one  of  the  roles  of  research  and 
development  in  engineering  is  to  reduce  the  level  of  uncertainty  in  design,  through  the 
improvement  and  refinement  of  information  needed  for  the  design  process.  In  this 
regard,  the  benefits  of  engineering  research  may  be  assessed  in  terms  of  the  savings  in 
design  and  life-cycle  cost  accruing  directly  from  an  investment  in  research.  The 
effectiveness  of  research  and  new  technology  may  be  analyzed  from  the  stand-point  of 
a  tradeoff  between  the  savings  in  life-cycle  cost  versus  the  research  and  development 
investment  required  to  achieve  such  savings. 

In  this  chapter,  the  cost  effectiveness  of  research  in  further  improvement  in  life- 
cycle  cost  is  examined  for  the  journal  bearing  examples.  To  do  this,  investment  in 
research  and  its  subsequent  benefits  reflected  in  the  improved  cost-rate  and  product  life 
are  evaluated.  The  cost  and  benefit  of  developing  the  thermo-hydrodynamic  lubrication 
model  is  estimated.  It  is  shown  that  the  investment  in  R&D  can  both  increase  the 
profits  to  the  manufacturer  and  savings  to  the  user,  though  it  results  in  higher  priced 
bearing.  The  analysis  also  results  in  new  price  for  the  bearing.  A  schematic  diagram 
in  Figure  8.1  illustrates  this  concept. 
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Product  specifications 


I 


Product/process  design  for  optimum  life-cycle  using 
dynamic  programming  or  multi-level  optimization 


I 


Evaluate  sensitivities  of  optimum  solution  to  important 
decision  variables  and  constraints  in  the  life-cycle 


T 


Evaluate  cost  and  benefits  of  investment  in  R&D  and 
new  technologies  by  evaluating  sensitivity  of  optimum 
solution  to  life-cycle  parameters  and  constraints 


T 


Decisions  to  invest  in  R&D  and  new  technologies 


I 


Improved  product/process  design 


Figure  8.1:     Sensitivity  and  cost-benefit  analysis  for  lowering  the  life-cycle  cost. 
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8.1    Cost/Benefit  Analysis  of  Research 

In  this  section  an  example  illustrating  the  use  of  cost/benefit  analysis  for 
decisions  in  investment  in  research  and  development  is  presented.  Suppose  that  the 
thermo-hydrodynamic  model  of  lubrication  did  not  exist,  but  it  is  known  through 
experimentation  that  the  hydrodynamic  model  of  lubrication  does  not  produce  accurate 
results.  It  is  also  known  that  the  deviation  of  experimental  results  from  the  predictions 
made  by  the  hydrodynamic  model  is  due  to  the  thermal  effects  in  the  fluid  film.  Hence, 
a  research  effort  has  to  be  undertaken  in  order  to  develop  the  experimental  setup  and 
investigate  the  thermal  effects  in  the  fluid  film.  The  cost  and  benefit  evaluation  of 
such  an  R&D  effort  for  developing  the  thermo-hydrodynamic  lubrication  theory  is 
presented  in  this  section. 

We  will  now  compare  the  two  models  of  lubrication  for  a  bearing  with  speed 
N=5000  rpm  and  load  ff=1000  lb.  The  rest  of  the  design,  manufacturing,  quality 
assurance,  and  operations  and  support  parameter  are  taken  to  be  same  as  in  example  1 
of  the  last  chapter.  Table  8.1  lists  the  important  parameters  for  the  optimum  bearings 
for  both  hydrodynamic  (HD)  and  thermo-hydrodynamic  (THD)  models. 

From  Table  8.1,  it  can  be  seen  that  the  bearing  designed  using  THD  model  has 
a  life  of  Lthd=25000  hours,  whereas  the  bearing  designed  using  HD  model  actually  has 
Lhd=6250  hours  of  life.  Noting  that  the  life  of  THD  bearing  is  four  times  that  of  HD 
bearing,  by  the  time  an  optimum  THD  bearing  is  replaced  by  the  consumer,  four  HD 
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Table  8.1:       Optimum  HD  and  THD  bearings  for  bearing  speed  N=5000  rpm,  and 
load  W=1000  lb.   (HD=hydrodynamic,  THD=thermo-hydrodynamic). 


Cases 

Optimum 
HD  Bearing 

Optimum 
THD  Bearing 

THD 

Simulation  of 

Optimum  HD 

Bearing 

Parameters 

Manufacturing  Cost  ($), 

183.10 

325.00 

— 

Internal  Quality  Cost  ($), 

2.32 

7.78 

— 

c  +c 

185.42 

332.78 

Oil  Cost  ($),  Coil 
Energy  Cost  ($),  Ce 

565.00 

685.00 

35.00 

975.19 

844.05 

3.97 

External  Quality  Cost 

(J),      CgQ 

5.24 

85.55 

Life  (hr) 

25000 

25000 

6250 

L/D 

0.6948 

.6868 

— 

Initial  Clearance  (in),  Q 

.002594 

.00228 

— 

Oil  Viscosity  (Reyn),  uavg 

Gu=.  31757 

Gu=  1.4364 

Tolerance  (in),  AD 

.001 

0.0005 

— 

Cost-Rate  ($/hr) 

.0692 

0.07789 

0.6378 

bearings  have  been  replaced.  Hence  the  cost/benefit  analysis  is  based  on  the  longer  of 
the  two  lives,  i.e.  25000  hours.  Assuming  that  the  bearing  is  operated  on  average  for 
h=4.28  hours  per  day,  the  lives  of  bearings  in  years  are  1-^=16  year  and  Lhd=4  years. 
The  analysis  is  based  on  the  1st  year  of  sales  of  THD  bearing,  and  1st,  5th,  9th,  and 
13th  year  of  sales  of  HD  bearing. 
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8.2   R&D  Cost 

The  cost  of  research  and  development  (CR&D)  can  be  estimated  as  follows. 


Labor  Cost 

1  theoretical  engineer  100  $/hr 

1  experimental  engineer         100  $/hr 

2  technicians  100  $/hr 
Cost  of  6  months  of  work     300  $/hr 

x  160  hr/month 

x  6  months  =  288,000  $ 

Equipment  cost  estimate  300,000  $ 

Total  equipment  and  labor  cost  estimate  6xl05  $ 


8.3   Manufacturer's  Profit 

When  the  optimum  HD  bearings  are  produced,  future  worth  at  year  v  of 
manufacturer's  profits  in  the  first  year  at  year  v  is  given  by 

^«/  =  (pw-cwj«w(i+iy 

Future  worth  at  year  y  of  manufacturer's  cumulative  profits  over  the  time  horizon  of 
v  years  is  given  by 

Phd,y  =  (Phd  -  Cmhd)nhd  [(1+iY],  if  y  £  Lu 

Pu.,  =  iPu  -  Cmhd)nhd  {(l+iy  +  (\+iyLhd],  if  Lu£  y  <  2L 


he! 
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Pm,  =  (pu  -  Cmhd)nM  [(l+ty  +  (l+iT,M  +  (l+*T2n  if  2LM<  y  <  3LM 

P*y  =  (Pu  -  C„J>nu  [(\+if  +  (\HTLhd  +  (l+iT2LU  +  (Wr3Lhd],     if  3IW<  y  <  4^ 

Future  worth  of  manufacturer's  profits  (due  to  1st,  5th,  9th  and  13th  years  of  sales)  at 
year  Llhd  of  the  profits  over  the  time  horizon  of  Llhd  years  is  given  by 

Phuh  =  (Pu  -  Cmhd)nM  [(\+i)lM  +  (l+i)Llhd-Lhd  +  (l+if^2^  +  (i+i)L'hd-3Lhd) 

Suppose  that  the  manufacturer  invests  in  research  and  development  and  the  thermo- 
hydrodynamic  model  is  developed.  As  a  result,  now  the  bearing  is  designed  based  on 
THD  model,  and  the  future  worth  at  year  y  of  the  manufacturer's  profits  at  year  v  is 
given  by 

Plhd.y  -  [~CMD  +  (pM  -  CmlJn,J  (1+/7,    for  0  <y  <  Llhd. 
where 

Phd       =  sale  price  of  the  bearing  designed  with  HD  model 

Cmhd     =         manufacturing  cost  of  the  bearing  designed  with  HD  model 

nhd       =         Number  of  units  sold  per  year  for  bearing  designed  with  HD 

model 
plhd       =  sale  price  of  the  bearing  designed  with  THD  model 

Cmlhd     =  manufacturing  cost  of  the  bearing  designed  with  THD  model 

nlhd       =  Number  of  units  sold  per  year  for  bearing  designed  with  THD 

model 
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i  =  annual  interest  rate 

(F/A,i,y)  =         Factor  for  future  worth  of  an  annuity  of  y  years  with 

yearly  interest  rate  i  and  is  given  by 


(F/A,  i,  y)=± — '- 


(P/A,i,y)  =  Factor  for  present  worth  of  an  annuity  of  y  years  with 

yearly  interest  rate  i  and  is  given  by 

(FM,  i,  y)=± — '- 

i(i+iy 

Hence,  change  in  manufacturer's  profits  in  year  y,  when  the  bearing  design  basis  is 
changed  from  HD  to  THD  model,  is  given  by 

&*y   =  *thd,y  "  *  hd,y 

8.4   Consumer's  Savings 

When  the  bearing  designed  using  HD  model  is  used,  the  future  worth  at  year  y 
of  initial  cost  to  consumers  can  be  given  by 

Clhd,y  =  (pM  +  pdnu  [(1+01,  for  0  <  y  <  Lhd. 

Clhd,y  =  (phd  +  A)"«  [(l+'T  +  (l+/r/J1,  for   Lhd<  y  <  2Lhd 

CMy  -  ipM  +  p)nM  [(]+iY  +  (\+irLhd  +  (l+iT2IM],  for   2Lhd<  y  <  3Lhd 
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CMy  =  (Phd  +  pfru  [(1+/T  +  (l+irLhd  +  {\W2Lhd  +  (l+/m      for   3Lhd<  y  <  4L 


hd 


and  future  worth  at  year  y  of  cumulative  operation  costs  can  be  given  by 


CoHy  =  ColjihdxS65h  (F/AXy) 


Hence  the  total  cost  to  consumers  at  year  y,  when  HD  model  is  used,  can  be  given  by 


*-hd,y  *-ihd,y   +    ^-ohd.y 


When  the  thermo-hydrodynamic  (THD)  model  is  adopted  as  the  design  basis,  the  new 
initial  costs  to  consumers  are  give  by 

Cuhd.y  =  (Pm  +  PMhd  [(l+Oyl  M   0  <  y  <  Llhd 
and  the  future  worth  at  year  v  of  operation  cost  can  be  given  by 

Co,hd.y  =  Colhdnlhdx365h  (F/AXy),  for  0  <  y  <  Llhd. 
The  total  cost  to  consumers  at  year  y,  when  HD  model  is  used,  can  be  given  by 

Clhd  =   (plhd+p)n,Jl+i)y  +  Colhdnlhdx365h  (F/AXy) 
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Hence  the  savings  to  consumers,  by  changing  the  design  model  from  HD  to  THD,  can 

be  expressed  as 

AS  =  Cm  -  Clhd. 


where 


Pi 


C 


ohd 


c. 


olhd 


Installation  cost  of  the  bearing  at  the  consumers'  site 
Operation  cost  of  HD  bearing  per  hour  of  operation 
Operation  cost  of  THD  bearing  per  hour  of  operation 


8.5    Price  Determination 

In  these  cost  formulations  the  known  quantities  are  CR&D,  Cmhd,  Cmlhd.  The 
unknown  quantity  is  the  price  the  manufacturer  should  charge  for  the  bearing  designed 
with  THD  model  (plhd).  Three  special  cases  are  considered  in  order  to  determine  the 
THD  model  price  pthd. 

Case  (a):  Here  it  is  assumed  that  over  the  time  horizon  of  interest  (L^), 
manufacturer's  profit  rise  is  equal  to  the  savings  to  consumers.  Here  basically,  the 
price  is  set  in  such  a  way  that  manufacturer  passes  half  the  profit  increase  to  the 
consumers.    This  can  be  expressed  as 


AP  =  AS 
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Case  (b):    Here  it  is  assumed  that  the  manufacturer  does  not  pass  any  profit 

increase  to  consumers.    Thus  the  price  is  set  such  that 


AS  =  0 


Case  (c):  Here  it  is  assumed  that  the  manufacturer  passes  all  the  profit  increase 
to  consumers.    Thus  the  price  is  set  such  that 


AP  =  0 


These  three  cases  give  three  different  equations,  which  are  solved  for  three 
different  values  of  the  price  of  the  THD  bearing  model.  These  prices  along  with  the 
manufacturer's  profit  and  consumers'  savings  are  listed  for  two  annual  interest  rates  of 
12%  and  6%  in  Table  8.2  and  8.3  respectively. 

The  calculations  show  a  steep  rise  in  price  of  the  bearing.  For  case  (a)  with  the 
interest  rates  i=\2%  and  6%,  the  price  of  the  bearing  is  determined  to  be  $1171  and 
$1526  respectively.  Even  after  a  price  rise,  both  manufacturer's  profits  and  the 
consumers'  savings  increase.    Thus  the  investment  in  R&D  is  justified. 
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Table  8.2:  Manufacturer's  profit  increase  and  consumers  savings  for  different  prices 
of  bearing  designed  with  thermo-hydrodynamic  model.  Interest  rate, 
Z=12%. 


Price  of  bearing 
designed  with 

thermo- 
hydrodynamic 
bearing  p^,  $ 

Increase  in 

manufacturer's 

profit 

AP,  $ 

Increase  in 

consumers' 

savings 

AS,  $ 

Case  (a):  1171.36 

Case  (b):  1798.03 

Case  (c):  544.69 

1.89xl09 

3.77xl09 

0.0 

1.89xl09 
0.0 

3.77xl09 

Table  8.3:  Manufacturer's  profit  increase  and  consumers  savings  for  different  prices 
of  bearing  designed  with  thermo-hydrodynamic  model.  Interest  rate, 
i=6%. 


Price  of  bearing 
designed  with 

thermo- 
hydrodynamic 
bearing  p^,  $ 

Increase  in 

manufacturer's 

profit 

AP,  $ 

Increase  in 

consumers' 

savings 

AS,  $ 

Case  (a):  1526.35 

Case  (b):  2478.81 

Case  (c):  573.89 

1.21xl09 

2.42x1 09 

0.0 

1.21xl09 
0.0 

2.42x10" 

CHAPTER  9 
SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


The  successful  development  and  implementation  of  the  design  methodology  has 
brought  into  focus  a  variety  of  observations.  This  chapter  starts  with  a  summary  of  the 
research  reported  in  this  thesis.  Some  of  the  observations  that  became  evident  during 
the  course  of  the  research,  and  their  implications  on  future  research  are  then 
enumerated.  Concluding  remarks  include  a  scrutiny  of  the  modeling  of  different  stages 
in  the  life-cycle,  the  solution  strategies,  the  modifications  needed  to  use  the  decision 
system  for  a  different  design  domain,  and  experiences  with  the  design  decision  system 
implementation. 

9. 1    Summary 

In  this  section  the  research  reported  in  this  thesis  is  summarized.  An 
investigation  of  the  concept  of  optimum  product  design  for  life-cycle  cost  was 
undertaken.  Most  of  the  product  life-cycle  cost  is  committed  at  the  design  stage,  hence, 
it  is  appropriate  that  the  life-cycle  considerations  be  included  at  the  design  stage  itself. 
The  idea  is  to  investigate  those  parameters  and  constraints  that  effect  the  life-cycle  cost 
of  the  product  at  the  design  stage  itself. 
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The  life-cycle  was  divided  into  four  stages  —  1.  research,  development  and 

design,  2.  manufacturing,  3.  quality  assurance,  and,  4.  operations  and  support.     A 

journal  bearing  was  chosen  as  the  design  domain  because  it  is  simple  in  configuration 

with   few  decision   variables,    yet   presents   all   the   difficulties,    nonlinearities    and 

constraints  of  a  complex  design  problem.   Its  performance  characteristics  are  also  very 

sensitive  to  the  decision  variables. 

The  life-cycle  cost  was  considered  to  consist  of  R&D  cost,  manufacturing  cost, 
in-plant  (internal)  quality  assurance  cost,  external  quality  cost  (cost  to  consumers),  and 
operations  and  support  cost.  Manufacturing  cost  was  modeled  using  those  design 
parameters  that  affect  the  cost  (length,  diameter,  tolerance,  material). 

The  operation  and  support  cost  was  taken  to  be  the  cost  of  energy  loss  (due  to 
friction)  and  the  oil  replacement  cost  through  the  life  of  the  bearing.  In  order  to 
evaluate  the  operations  and  support  cost,  the  bearing  and  oil  deteriorations  are  simulated 
through  its  life.  The  bearing  life  could  be  either  its  physical  life,  when  it  starts  to 
violate  any  of  the  constraints,  or  its  economic  life,  when  the  cost-rate  (life-cycle  cost 
divided  by  the  number  of  hours  of  operation)  curve  reaches  its  minimum.  Modeling 
of  wear  and  deterioration  of  oil  was  one  of  the  most  sensitive  part  of  the  life-cycle  cost. 

Internal  quality  costs  (in-plant)  were  modeled  assuming  that  the  process  is  being 
controlled  by  either  an  automatic  system  or  a  quality  control  chart  such  as  X  chart. 
External  quality  costs  (to  consumers)  measure  the  cost  of  variability  of  performance 
from  the  optimum  design.  The  measurement  method  of  external  quality  costs  in  this 
research  is  represents  an  improvement  over  the  previous  methods  published  in  literature. 
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For  the  life-cycle  optimization,  three  different  solution  approaches  were 
investigated  and  implemented.  The  first  two  approaches  tried  were  algorithmic  and  the 
last  approach  combines  algorithmic  approach  with  the  heuristic  approach.  The 
algorithmic  optimization  approaches  implemented  are  very  general  in  nature.  The 
decision  system  developed  is  highly  modular  and  could  be  easily  modified  for  a 
different  design  domain. 

The  first  optimization  approach  was  to  treat  the  life-cycle  optimization  as  a 
dynamic  programming  problem.  Dynamic  programming  approach  is  applicable  and  is 
very  efficient  when  a  problem  can  be  thought  of  as  consisting  of  serially  connected 
black  boxes  (stages)  with  output  of  each  stage  being  fed  as  input  to  the  next  stage. 
Each  stage  has  an  objective  function,  input  parameters,  output  parameters,  constraints 
and  a  set  of  decision  variables.  The  objective  is  to  choose  all  the  decision  variables 
subject  to  given  constraints  so  as  to  optimize  the  total  objective  function,  which  is  the 
sum  of  objective  functions  of  all  the  stages.  A  very  general  dynamic  programming 
based  optimization  system  was  developed  where  the  number  of  stages  and  the  number 
of  decision  variables  at  each  stage  can  be  varied. 

A  very  general  nonlinear  constrained  optimization  module  was  developed  which 
can  handle  a  design  vector  of  arbitrary  size  and  any  number  of  constraints.  The 
constrained  optimization  problem  is  first  converted  into  an  unconstrained  optimization 
problem  by  adding  a  combination  of  interior  and  exterior  penalty  functions  to  the 
objective  function,  and  then  by  using  augmented  Lagrange  multiplier  method.    The 
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unconstrained  problem  is  optimized  using  conjugate  directions  or  variable  metric 

methods. 

The  second  optimization  strategy  used  was  multi-level  optimization  approach, 
where  the  design  optimization  is  decomposed  into  many  levels.  Here  the  design  vector 
is  decomposed  into  a  global  design  vector  and  level-specific  design  vectors.  In  this 
research,  the  design  vector  was  decomposed  into  two  levels  —  first  level  being  the 
design  parameters  and  the  second  level  being  the  tolerance.  At  each  level  of 
optimization,  calls  to  nonlinear  constrained  optimization  modules  are  made. 

The  third  optimization  approach  tried  was  a  combination  of  heuristic  and 
algorithmic  procedures.  Ideally,  the  combination  of  the  two  procedures  in  a  decision 
system  can  provide  the  best  design  tools.  The  heuristic  system  generates  initial  design 
based  on  the  application  domain  of  the  bearing  specifications  supplied  by  the  designer. 
The  heuristic  module  can  also  process  the  feedback  about  constraint  violations  and 
modify  the  appropriate  parameters  so  that  the  design  remains  within  the  feasible 
domain.  The  heuristic  module  is  design  domain  specific  and  needs  to  be  completely 
rewritten  for  a  different  design  domain  —  but  it  can  serve  as  a  guideline  for  developing 
the  heuristic  module  for  a  new  design  domain. 

The  decision  system  developed  proved  to  be  a  useful  tool  in  the  optimum  design 
in  the  domain  of  journal  bearing.  Many  examples  of  parameter  and  tolerance  design 
for  a  journal  bearing  for  optimum  life-cycle  cost  were  presented.  It  is  seen  that  a  well 
defined  optimum  for  parameters  and  tolerance  exist.  Values  of  optimum  bearing  length, 
clearance,  viscosity,  and  tolerance  obtained  seem  to  be  very  reasonable. 
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The  research  illustrates  the  importance  of  incorporating  manufacturing,  quality 

assurance,  operations  and  support,  and  cost/benefit  analysis  considerations  at  the  design 

stage  decisions  itself.    Using  the  decision  system,  the  concept  of  identifying  the  need 

for  research  and  development  and  justifying  the  investments  in  R&D  were  investigated. 

The  cost/benefit  analysis  of  alternatives  in  the  life-cycle  identifies  the  parameters  and 

constraints  that  need  to  be  improved.   It  also  provides  the  justifications  for  investments 

in  R&D  and  new  technology. 

9.2   Concluding  Remarks 

In  this  section,  some  of  the  observations  that  became  evident  during  the  course 
of  the  research  are  presented.  Design  for  life-cycle  is  akin  to  the  concept  of  concurrent 
engineering  or  simultaneous  engineering.  In  very  complex  designs/projects,  it  may  not 
be  obvious  that  in  order  to  reduce  total  costs,  what  improvements  are  necessary.  Cost 
modeling,  optimization  and  then  sensitivity  analysis  can  immensely  help  in  determining 
where  to  invest  the  scarce  resources.  Cost/benefit  analysis  of  alternatives  available 
helps  a  decision  maker. 

Modeling  of  bearing  wear  and  deterioration  of  oil  was  the  most  sensitive  and 
time  consuming  part  in  the  development  of  the  decision  system.  An  insight  into  the 
lubrication  and  tribology  practice  helps  make  a  few  intelligent  assumptions,  which  can 
make  the  task  of  bearing  life  simulation  manageable.  The  implemented  decision  system 
calculates  the  objective  function  using  simulations  of  bearing  life.   Performing  the  life 
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simulations  of  the  bearing  was  found  to  be  very  demanding  on  the  computational 

resources.    Each  simulation  requires  an  average  of  10  minutes  on  a  80826  processor 

based  PC  with  a  80287  math  coprocessor.   But  it  was  necessary  for  accurate  prediction 

of  physical  parameters  and  costs  of  the  life-cycle.   Because,  thousands  of  such  bearing 

simulations  were  to  be  conducted,  a  graphics  plotting  module  was  developed  and  was 

incorporated  into  the  decision  system,  thus  providing  an  instant  access  to  the  progress 

of  the  life-cycle  simulation. 

As  a  solution  strategy,  the  dynamic  programming  was  found  to  be 
computationally  slower  than  the  multi-level  optimization  approach.  But  in  large 
problems  with  many  more  stages,  dynamic  programming  could  prove  to  be  more 
efficient  than  multi-level  optimization  approach.  The  investigation  also  revealed  that 
a  combination  of  heuristic  and  algorithmic  procedures  in  a  decision  system  can  provide 
the  best  design  tool.  One  limitation  of  the  heuristic  based  approach  is  that  for  each  new 
design  domain,  the  module  must  be  rewritten. 

The  life-cycle  optimization  problem  for  a  different  design  domain  can  be 
accommodated  by  readily  modifying  the  algorithmic  decision  system.  The  algorithmic 
part  of  the  decision  system  is  modular  in  nature  and  can  be  modified.  The  dynamic 
programming  module  is  also  very  general  and  it  can  incorporate  an  arbitrary  number 
of  stages  and  an  arbitrary  number  of  decision  variables  at  each  stage  —  subject  to 
memory  and  storage  capabilities  of  the  computer. 

In  the  multi-level  optimization  approach,  the  design  domain  optimization  must 
be  decomposed  into  many  levels.    Hence  the  designer  has  to  decompose  the  design 
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vector  into  a  global  design  vector  and  level-specific  design  vectors.  Hence  the  designer 

has  to  incorporate  this  in  the  main  program  of  the  design  module.   He  also  has  to  make 

calls   to   optimization    modules    for   each    level    of  optimization.       The   nonlinear 

optimization  modules  require  to  be  fed  the  normalized  values  of  the  design  vector.   A 

designer  needs  to  replace  the  following  modules  in  the  algorithmic  part  of  the  system 

in  order  to  use  the  decision  system  for  a  different  design  domain. 

Each  design  or  life  cycle  constraints  is  represented  by  a  subroutine  in  the 

system.    Any  number  of  constraints  can  be  incorporated  as  a  design  domain 

requires. 

The  input/output  module  provides  the  input/output  venue  for  the  parameter 

values  and  design  vector. 

•  The  journal  bearing  simulation  module  calculates  the  operation  costs. 

The  internal  quality  cost  evaluation  module  can  be  replaced  if  a  different  quality 
control  method  is  being  used. 

•  The  manufacturing  cost  evaluation  module  can  also  be  changed  to  suit  the 
problem  under  consideration. 

9.3    Recommendations  for  Further  Research 

This  section  outlines  the  recommendations  for  further  research  and  possible 
extensions  to  the  design  decision  system  which  has  been  implemented.  Better 
tribological  models  for  wear  are  needed.    Besides  the  improvements  in  modeling  the 
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manufacturing  cost,  quality  cost,  reliability  cost  discussed  in  the  following  sections, 

additional  factors  or  stages  can  be  included  in  the  life-cycle.    For  examples,  cost  of 

disposal  of  the  replaced  bearing  could  be  incorporated. 

9.3.1    Manufacturing  Cost 

Better  cost  models  for  manufacturing  processes  are  needed.   The  manufacturing 
cost  could  include  the  following  factors: 

•  The  part  parameters  such  as  specified  surface  roughness,  surface  coatings  etc. 

•  The  part  material  characteristics  such  as  machinability,  material  hardness,  any 
heat  treatment  needed  and  material  cost. 

•  Lot  size  of  the  parts. 

•  Process  parameters  (or  machining  parameters)  such  as  cutting  speed,  depth  of 
cut,  feed  rate,  processing  time,  and  tool  material. 

•  Part  handling  parameters  such  as  any  special  jigs  and  fixtures  needed,  setup 
time,  and  material  handling  time. 

•  The  machine  capacity  including  factors  such  as  power  rating  and  speed. 

•  The  machine  capability,  the  ability  of  the  machine  to  attain  design  specifications, 
such  as  geometry,  dimensional  tolerances,  and  surface  roughness. 

The  machine  availability,  a  factor  varying  between  zero  and  one,  which 
represents  whether  a  machine  is  available  for  use  when  requested  or  it  is  highly 
in  demand.  This  factor  then  can  be  used  for  applying  a  penalty  for  choosing  a 
relatively  unavailable  machine. 
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Based  on  these  manufacturing  parameters  and  selected  manufacturing  process, 

one  could  calculate  the  fraction  of  the  products  (or  parts)  being  produced  which  are 

defective.    This  information  is  very  useful  in  evaluating  the  quality  costs  accurately. 

One  could  include  and  make  use  of  a  plant-wise  database  on  material  costs  and 

material  availability.   Similarly  one  could  also  make  use  of  the  databases  on  production 

machine   capacity,   capability  and  availability   in  order  to  make  optimum  process 

plans/scheduling  for  production.     This  data  not  only  can  be  used  for  evaluating  the 

costs,  but  also  it  can  serve  as  constraints  imposed  by  the  manufacturing  function  to  the 

designer.    Better  accounting  practices  can  provide  the  data  on  costs  associated  with 

these  factors.    Already  existing  modules  which  produce  and  optimize  the  production 

schedules  could  be  incorporated.    Cost  of  research  and  development,  innovations  and 

investments  in  new  manufacturing  processes  and  machines  could  also  be  incorporated. 

9.3.2  Quality  Costs 

Depending  on  the  design  domain,  there  can  be  different  quality  control  plans 
used  for  example  attributes  or  variable  control  charts.  Integrated  system  could 
incorporate  attribute  control  charts,  and  process  control  based  on  it. 

9.3.3  Reliability  Cost 

Following  is  a  brief  presentation  of  the  way  reliability  and  maintainability 
consideration  could  be  incorporated  in  life-cycle  analysis.  The  product  reliability  and 
maintainability  cost  is  the  cost  incurred  due  to  scheduled  preventive  maintenance,  cost 
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of  failures  and  the  cost  of  unscheduled  maintenance  following  failures.     It  can  be 

expressed  as 

CR  =  Nff  +  NmCm,  (9.1) 

where  Nf  and  Nm  are,  the  expected  number  of  failures  (or,  unscheduled  corrective 
maintenance)  and  the  number  of  replacements  of  unfailed  parts  (or,  preventive 
maintenance)  during  the  life-cycle  of  the  product,  respectively.  Cm  is  the  cost  of 
replacing  the  enfilade  parts  including  the  loss  due  to  product  (or  system)  unavailability 
during  preventive  maintenance.  Cf  is  the  cost  of  parts  failures  and  subsequent 
replacement  during  corrective  maintenance,  including  the  damage  to  other  parts  and  the 
system.  For  the  journal  bearing  example,  "parts"  may  include  any  of  the  bearing  parts 
or  the  lubricant  oil.   Total  number  of  parts  replacements  is  given  by 

N  =  Nf+Nm.  (9.2) 

Assuming  that  the  total  life-cycle  time  is  very  long  compared  to  the  mean  time  between 
replacements  (MTBR),  we  have 

MTBR  =  t/N.  (9.3) 

If  R(f)  be  the  part  reliability  at  time  t,  and  assuming  that  the  part  is  replaced  after  it  has 
been  in  operation  for  time  T,  then  MTBR  can  be  expressed  as 
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7 


MTBR  =    [/?(/)  dt.  (9.4) 

Therefore  from  the  above  two  equations,  the  number  of  parts  replacements  over  a  long 
time  span  /  is 

N  -  -J—  -  -J—. 

MTBR        J  (9.5) 

[R(t)dt 

Since  only  a  fraction  of  the  parts  will  survive  to  the  next  preventive  maintenance,  the 
number  of  unfailed  parts  replaced  during  the  life-cycle  of  the  product  is  given  by 

N    =  R(T)N  =     //?(r)  . 

r  (9-6) 

\R{t)dt 


Moreover,  since  a  fraction  \-R(T)  of  the  parts  will  fail,  the  number  of  failed  parts 
replacement  is  given  by 


N,  =  [l-R(T)]N  =   t[l   /?(r)]. 

(9.7) 
R{t)dt 


With  expressions  for  A^  and  Nm  determined,  total  reliability  cost  CR  becomes  a  function 
of  replacement  time  T.  Now  we  can  minimize  the  cost  of  reliability  by  making 
dC,/dT=0,  which  can  be  reduced  to 
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C 


R(T)  +  X(T){R(t)dt  -  1  =  m—  (9-8) 

"0  /      « 


which  has  to  be  solved  numerically  if  the  optimum  T  is  to  be  estimated.    Here  R{t)  is 
the  reliability  function  and  X(t)  is  the  failure  rate  defined  as  follows 


=  _J_dm  (9.9) 

R(t)    dt 


If  one  assumes  that  the  failure  rate  follows  the  widely  used  Weibull  distribution,  which 
can  model  a  wide  range  of  failure  modes,  we  have  reliability  function  as 


-(')"  (9.10) 


and  failure  rate  given  by 


e 


m-l 


(9.11) 


Here  9  is  a  scale  parameter  having  dimension  of  time.  Parameter  m  is  a  shape 
parameter  which  depends  on  the  mode  of  failure.  For  failure  rates  increasing  with  time, 
m>l,  and  wear  is  one  of  the  modes  of  failure.  For  failure  rates  constant  in  time,  m=l, 
the  mode  of  failure  is  random  and  the  failure  rate  has  an  exponentially  decreasing 
distribution.  Assuming  that  failure  cost  Cf  is  much  greater  than  the  preventive 
maintenance  cost  Cm,  and  after  performing  a  few  algebraic  manipulations,  the  optimum 
time  of  part  replacement  (MTBR)  is  given  by  [Lewi87] 


MTBR  *  *  9 


1     C 

W-l     C, 


197 


(9.12) 


Besides  estimating  the  optimum  time  to  replace  parts,  it  is  also  important  to  find 
the  optimum  age  of  the  whole  journal  bearing  system.  Many  studies  have  been 
conducted  to  develop  reliability  based  probabilistic  failure  models  based  on  the  wear 
rates  of  various  mechanical  components.  The  references  [Skew85,  Raze86,  Raze87, 
Rhod88,  Lamb88,  Nels89]  suggest  the  reliability  based  design  methods  for  mechanical 
equipment  such  as  clutches,  brakes,  seals,  bearings,  valves  and  gearboxes. 

As  an  alternative  to  the  deterministic  models  used  to  model  the  maintenance 
interval  time  (oil  replacement)  in  the  thesis,  probabilistic  models  of  reliability  and 
failure  rates  could  be  used  to  estimate  the  optimum  maintenance  interval  tm  and  the  life 
t,.  A  number  of  studies  have  been  conducted  in  reliability  based  methods  for  estimating 
the  optimum  maintenance  interval  and  optimum  replacement  interval  for  non-maintained 
and  maintained  systems.  Some  examples  of  these  studies  can  be  found  in  references 
[Jard73,  Bloc86,  Park88,  Regu83,  Geur83]. 

Suppose  we  denote  the  reliability  as  a  function  of  time  of  a  system  without 
maintenance  as  /?(/),  where  t  is  the  operation  time  of  the  system.  If  we  perform 
maintenance  on  the  system  at  intervals  T,  then  for  /  <  T,  maintenance  will  have  no 
effect  on  reliability.    That  is,  the  reliability  of  the  maintained  system  RJj)  is  given  by 

RJf)  =  R{t),   forO<t<T.  (9.13) 
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After  maintenance  at  time  T,  in  the  interval  T  <  t  <  2T,  the  reliability  is  the  product  of 

the  probability  R{t)  that  the  system  survived  to  time  T,  and  the  probability  R(t-T)  that 

a  system  as  good  as  new  at  time  T  will  survive  for  a  time  (7-7)  without  failure,  i.e. 

RJt)  =  R(T)R(t-T),   foxT<t<2T  (9.14) 

The  same  argument  can  be  used  repeatedly  to  obtain  a  general  expression  for  the 
reliability  of  the  maintained  system, 

RM  =  R(T)n  T(t-nT),    for  nT  <  t  <  (n+l)T,    «=0,1,2,...  (9.15) 

Thus,  mean  time  to  failure  (MTTF)  for  a  system  with  preventive  maintenance  can  be 
expressed  as 


RM  (9.16) 


""*  '  \U,)<"  -  Tlffl    -  m 


Total  cost  of  the  system  operation,  maintenance  and  failure  based  replacement  can  be 
estimated  as 

CtC)  =  NJf)CJf)  +  [1-RJ.t)]  C/0  +  C/0  +  CM  (9.17) 
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where 

Nm        =      number  of  maintenance  operations  performed  by  time  /, 

CJj)   =       cost  of  each  preventive  maintenance  of  the  system, 

Cft)   -        cost  of  failure  of  the  system  (includes  down  time  and  damage  to  other 

systems  but  not  the  cost  of  replacement), 
Cfjt)   =        initial   investment   in  the   system   (includes   the   initial   cost,   cost   of 

installation,  and  cost  of  down  time  for  installation), 
C0(f)   =       other  operating  costs  (in  case  of  bearing  costs  such  as  energy  loss  due  to 

friction,  oil  cleaning/filtering  costs). 
Substituting  the  following  expression  for  operation  time  t  in  the  above  equation 

t  =  NmT+  te,    with   te  «  t,   or  te  -►  0  (9.18) 

to  find  the  total  cost  per  unit  time  of  operation, 

C  Jt,T)  =  SA°   +  \\-R(T)'IT}El  +  EE  +  £W  (9.19) 

"A  T         L  J  t  t  t 

Minimization  of  this  quantity  can  give  the  optimum  preventive  maintenance  time  Tm' 
and  optimum  life  //  from  the  expression  for  MTTF  in  equation  9.16. 

9.3.4   Sensitivity  Analysis  of  Optimum  Solution 

After  finding  the  optimum  solution  to  life-cycle  problem  an  additional  step  that 
can  be  highly  valuable  is  the  sensitivity  analysis  of  optimum  solution  to  problem 
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parameters  (design  parameters  and  design  constraints).  Such  a  sensitivity  analysis  of 
the  optimum  solution  is  useful  in  order  to  perform  the  cost/benefit  analysis  of 
investment  in  R&D  and  new  technologies.  An  additional  application  of  sensitivity 
analysis  is  to  provide  a  basis  for  formal  decomposition  of  the  multi-stage  (serial) 
problem  to  a  multilevel  (hierarchical)  problem.  Following  is  a  brief  discussion  of  the 
sensitivity  analysis  problem  and  its  formulation. 

Sobieszczanki-Sobieski  et  al.  [Sobi82,  Sobi83]  give  a  method  of  estimating  the 
sensitivity  derivatives  directly  from  a  set  of  constrained  optimization  problems.  The 
equations  can  be  used  regardless  of  the  type  of  optimization  algorithm  that  is  used  to 
arrive  at  the  optimum  solution.  Sobieszczanki-Sobieski  et  al.  [Sobi82]  also  state  the 
conditions  for  solvability  of  the  sensitivity  equations  developed. 

There  are  at  least  two  ways  of  deriving  the  equations  for  the  unknown 
sensitivities.  One  way  is  to  start  from  the  Lagrange  multiplier  equations  of  the 
constrained  minimum;  the  other  method  begins  with  the  extremum  conditions  of  the 
penalty  function.  In  either  case,  the  same  general  functional  relationships  are 
recognized.  The  first  method  as  described  by  Sobieszczanki-Sobieski  et  al.  is  now 
briefly  presented  as  applicable  to  the  life-cycle  cost  optimization  problem. 

Let  the  objective  function  be 

F(X„X2,...,X„  p„p2,...,pd  (9.20) 

and  the  constraints  active  at  the  optimum  point  be 

gj(X„X2,...,X„  p„p2,:,pk),   for  y'=l  to  m,  (9.21) 
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and  implicitly  at  the  optimum  point  X'(p,,p2,..pk).     Here  X„X2,..Xn  are  the  design 

variables,  and  P1.P2.-Pk  aie  me  parameters  of  the  problem.     There  is  no  need  to 

distinguish  among  the  equality  and  inequality  constraints  in  the  optimization  problem, 

because  all  the  constraints  active  (let  it  be  a  set  of  constraints  J)  at  the  constrained 

optimum  point  may  be  regarded  as  the  equality  constraints.    The  Lagrange  multiplier 

equations  satisfied  at  the  constrained  minimum  point  are 

VF(X)  +  £j  XjVg/X)  =  0,   for;  e  J  (9.22) 

gj(X)  =  0,  and  X,  =  0,   for  j  e  J  (9.23) 

It  is  understood  here  that  X  and  X  correspond  to  their  values  at  the  optimum  point  X* 
and  X*,  respectively.  Differentiating  the  first  equation  with  respect  to  a  problem 
parameter  p,  we  have 


*=i 


dX3Xk   +  £f    J  dXdXk 


exk  +  ^  dgpp  exj  +  ffpw  ^  ^    &gpo    _Q 

dp       ^T     dX.     dp         dXdp       fa    J  dXdp 

for  j  6  J  (9.24) 


Also  differentiating  the  constraint  equation  with  respect  to  parameter  p 


V^5  +  M=0  (9.25) 

tr  dxi  dp      dp 


forj  €  J 
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We  can  now  determine  the  value  of  the  terms  dX/dp,  and  dXJdp  for  i=\,..n  andy'GJ, 

by  writing  the  above  two  equations  in  the  matrix  form 


nxn         nxj 


dXIdp 
dX/dp 


'nxl 


Jx\ 


=  0 


(9.26) 


Here  the  Lagrange  multipliers  X*  are  either  determined  as  part  of  optimization  process, 
or  they  can  be  determined  from 

X*  =  -  [BT]-'[BT]  VF(X)  (9.27) 


A  fundamental  requirement  here  is  that  active  constraints  be  linearly  independent  of 
each  other.  We  can  now  calculate  the  derivative  of  the  objective  function  at  the 
optimum  point  as  follows 


dF(X) 


dp  dp  dp 


(9.28) 


Thus  the  sensitivity  of  the  objective  function  and  design  variables  with  respect  to 
parameter  p  can  be  determined.  The  information  now  available  can  be  used  to  estimate 
when  a  currently  active  constraint  in  J  will  cease  to  be  active  or  when  some  new 
constraint  will  become  active.  An  equality  constraint  (in  J)  which  is  presently  active 
will  become  inactive  when  its  Lagrange  multiplier  goes  to  zero.  Therefore  due  to  a 
change  Ap  in  the  parameter  p,  Lagrange  multiplier  changes  by  {dXJdp)Ap,  hence 


Xj  +  (dX/dp)Ap  =  0,   or,   Ap  =  -  X/(8X/dp) 


(9.29) 
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Similarly,  a  currently  inactive  constraint  will  become  critical  when  its  value  goes  to 

zero,  so 

gjW  +  \yg/X)-dXdp]Ap  =  0,   or,   Ap  =  -g/X)  I  [Vgj(X)-dX/dp]        (9.30) 

Thus  from  the  above  formulations,  the  sensitivities  of  optimum  design  variable  X*  and 
sensitivity  of  optimum  objective  function  F*  with  respect  to  a  parameter  p  can  be 
obtained.  The  sensitivity  analysis  can  thus  be  used  for  further  improving  the  life-cycle 
cost  of  the  products. 


APPENDIX  A 
MULTIDIMENSIONAL  INTERPOLATION 


In  the  dynamic  programming  approach,  in  order  to  compute  the  optimum  value 

function  for  a  stage,  a  general  multi-dimensional   interpolation  scheme  is  used  as 

described  in  this  appendix.  In  solving  the  problem  of  interpolating  a  function  <|>  which 
is  a  function  of  a  N  dimensional  vector,  it  is  assumed  that  a  grid  of  nodes  in  N 

dimensional  space  has  been  laid  with  the  node  coordinates  and  the  function  value  at 

those  node  points  given.    More  specifically,  let  v  be  a  function  of  a  N  dimensional 

vector  u  defined  by  some  unknown  function  <)>, 

v  =  <)>(«)  (A.1) 

where 

u  =  {u„u2,...,uN}.  (A.2) 

Furthermore  a  reasonable  number  of  discretizations  of  the  range  of  possible  values  for 
variable  un  and  the  corresponding  values  of  v  are  given.  Let  vectors  umin  and  umax 
denote  the  range  of  possible  values  of  vector  //  components  as, 

Umin  ~    \Ul,mirf->Un,mirt>->UN,mm}->  (A.3a) 
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and 


such  that 


"mm  \"l,ntav->Mn,mav-'"N,max)' 


u        <  u   <  u 
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(A.3b) 


(A.4) 


If  L„  equally  spaced  discretizations  are  done  for  the  component  u„  within  the  range 
("wAJ-  men  me  ^'m  discretization  of  u„  can  be  given  by 


«     =  u        +   ""■max"t/"'min  (M),      I  -■  \,...,L_;    n  -    1...JV.  (A.5) 


L  -1 


Let  all  discretizations  of  77'th  component  un  be  stored  in  vector  u„  as  follows 


u„  =  l«n/,.. .,uHit,...,uIKJ,     for  «=1,..JV 


(A.6) 


which  are  stored  in  a  Z,x7V  matrix  U  as  row  vectors  as  follows 


V 


«, 


« 


«. 


(A.7) 


The  total  number  of  the  sets  of  combinations  generated  by  such  discretizations  un  is 
given  by 


N 

k  =  n  l 

n 
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(A.8) 


This  means  that  the  discretizations  generate  a  total  of  K  nodal  points  in  N  dimensional 
space.  Furthermore  let,  these  K  nodal  points  uk  in  N  dimensional  space  be  stored  in  a 
KxN  matrix  U  with  rows  containing  the  nodal  points  as  follows 


U  = 


u. 


u, 


«l 


where 


uk  =  {"/*••  "„*•■  «Af»}»    for   fcl,-,K,    and  n=l,...JJ. 


(A.9) 


(A.10) 


The  corresponding  function  values  (at  nodal  point  k)  vk  are  stored  in  a  vector  V,  such 
that 


v*  =  <K«t).    for  k=\,...,K. 


(A.11) 


In  the  interpolation  problem  encountered  in  the  general  dynamic  programming  problem, 
matrices  U  and  U,  and  vector  V  are  assumed  to  be  given.  The  objective  is  to  evaluate 
the  unknown  function  v=$(u)  for  a  given  vector  u  using  interpolation  of  surrounding 
nodal   points  in  N  dimensional    space.      In  this  thesis,   both  quadratic   and  cubic 
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polynomial    interpolations    were    implemented.       We    will    denote    the    degree    of 

interpolation  (the  maximum  value  of  exponent  of  uj  by  p  —  value  of  p  being  2  for 

quadratic  interpolation  and  3  for  cubic  interpolation. 

Quadratic  Interpolation 

In  the  quadratic  interpolation,  if  the  vector  u  is  of  two  dimensions,  then  we  can 
define  a  following  interpolation  function 

<K")  =  a<>o  +  a,gu,  +  a0,u2  +  a„u,u2  +  a20u,2  +  a02u22 

+  a21u,2u2  +  al2u,u2  +  a22u,2u22  (A.  12) 

The  powers  of  the  u,  and  u2  in  quadratic  2  dimensional  case  given  by  this  equation  can 
be  put  in  a  tabular  form  as  in  Table  A.l.  Similarly,  when  vector  u  is  3  dimensional, 
the  coefficients  of  the  interpolation  polynomial  and  corresponding  exponents  can  be 
shown  in  a  tabular  form  as  in  Table  A.2  and  the  function  <)>(«)  is  given  by 

27 


♦(«)    -E'ki"!     "2     Ml"  (A'13) 

/=1 
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Table  A.l:       Exponents  of  polynomial  terms  for  2  dimensional  quadratic  interpolation. 


Polynomial 
Term  Number  i 

Exponents  of 
u,  u2 

1 

00 

2 

0  1 

3 

02 

4 

1  0 

5 

1  1 

6 

1  2 

7 

20 

8 

2  1 

9 

22 

By  induction,  for  quadratic  interpolation,  we  can  conclude  the  following  points 
about  the  interpolation  polynomial  <j>  which  is  a  function  of  N  dimensional  vector  u: 
The  polynomial  expression  for  <|>  for  quadratic  interpolation  has  IN  =  3N  terms. 
The  f  th  term  of  the  polynomial  has  a  coefficient  'a/  of  the  following  form 


a  =  a 

•  J,Ji-JN 


(A.  14) 

The   z"th  term   exponents    for  N  components    of  vector   u   are  j,,...,j„-..,jN 

respectively. 

The  set  [f'i—Jir..JN]  forms  a  number  in  base  TV  (binary,  tertiary  etc.)  which  is 

which  is  related  to  the  term  number  i  as  follows, 


i  =  Uf-Jm-JlAhm  N  +  1>      fOr  '-1,-,/ 
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Table  A.2:       Exponents  of  polynomial  for  quadratic  interpolation  in  3  dimensions. 


Polynomial 
Term  Number 

Coefficient 
of  z"'th  Term 

Exponents 
for 

i 

u,  u2  u3 

U1J2J3) 

1 

aooo 

000 

2 

aooi 

00  1 

3 

a<x>2 

002 

4 

aoio 

0  1  0 

5 

<*01l 

0  1  1 

6 

aOI2 

0  1  2 

7 

°020 

020 

8 

a02i 

02  1 

9 

<*022 

022 

10 

aioo 

1  00 

11 

awi 

1  0  1 

12 

<*I02 

1  02 

13 

<*iio 

1  1  0 

14 

0/7/ 

1  1  1 

15 

a, 12 

1  1  2 

16 

<*120 

1  20 

17 

<*121 

1  2  1 

18 

<*,22 

1  22 

19 

a200 

200 

20 

<*201 

20  1 

21 

a202 

202 

22 

<*210 

2  1  0 

23 

<*211 

2  1  1 

24 

<*212 

2  1  2 

25 

<*220 

2  20 

26 

a22l 

22  1 

27 

<*222 

222 
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Consider  a  hexahedron  or  a  brick  in  N  dimensional  space  (a  iV-cube)  which  has 
2N  vertices.  If  it  has  nodes  on  all  the  vertices,  on  mid-edges,  mid  faces  and  at  mid 
-volume,  then  the  total  number  of  such  nodes  is  given  by  7^=3".  In  one  dimension  the 
nodes  are  I ,=3,  in  2  dimensions  I2=9  and  in  3  dimensions  I3=27.  In  the  interpolation 
problem,  the  next  step  is  to  find  a  JV-cube  from  a  given  set  of  discretizations  in  U  such 
that  the  point  u  to  be  interpolated  lies  inside  this  JV-cube.  This  is  done  by  searching 
for  a  point  u°  through  the  given  discretizations  in  U  such  that  the  w'th  component  of 
u°  is  largest  available  discretization  point  but  smaller  than  the  n'th  component  of  u, 
i.e. 

«*«  *    »*fa.  *  *-  (A.\5a) 

for  I  €  [1.JLJ,    l„  e  [1,..I„],  and  l*l0  (A.  15b) 

This  is  repeated  for  all  the  N  components  (n=\,..N)  of  U  and  the  indices  tM  for  u^lno 
are  found.  We  will  denote  the  N-cube  vertex  («„,„,  n=\,..N)  thus  found,  by  u"  and  all 
the  indices  by  vector  {£„}.  This  completes  the  identification  of  all  the  indices 
({0i,.ioN)  for  one  vertex  of  the  cube.  The  coordinates  of  the  rest  of  the  vertices  of  JV- 
cube  are  found  by  incrementing  the  indices  one  at  a  time.  If  the  f  th  vertex  of  iV-cube 
«',  is  defined  as 

«',  =  {u'iV...,u'in,...,u'^}  (A.16) 
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then  N  indices  in  U  of  the  z"th  vertex  u ',  are  given  by 


u'     =  U  ,  (A.17) 


where  indices  lin  are  evaluated  using 

Hence  the  indices  in  1/  for  subsequent  vertex  /+1  are  found  by  just  algebraically  adding 
1  to  the  number  in  base  N  defined  by  indices  [I ,]base N  =  [lu,...,  t^...,  liN].  For  the  vertex 
i  of  the  A'-cube,  the  corresponding  row  k  in  matrix  U,  is  simply  the  decimal  (base  1 0) 
value  of  [t^ft,  i.e. 

Kbase  10  —     !*• /,/•••>  * i,rr-->  ' i.Nlbase  N  (A.  I") 

Finally,  for  quadratic  interpolation  in  N  dimensional  space,  the  polynomial  for 
function  <j>  is  of  the  following  form, 


♦<•>-      E      <v.A  <  «?.-•»»  (A20) 


Here,  note  that  the  number  of  polynomial  terms  IN  for  <|>  is  same  as  the  number  of  nodes 
defined  by  the  vertices  of  the  A'-cube  just  considered.  Hence  the  polynomial  of 
equation  (A.20)  can  be  used  for  interpolating  the  given  function  values  v',  at  the 
vertices   u' ,  (i=\,..IN)  of  the  A'-cube.      This  is  true  since  it  leads  to  a  set  of  IN 
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simultaneous  linear  equations  in  IN  unknown  coefficients  a,'s.    Let  the  multiplier's  at 

the  nodal  points  u  \  be  stored  as  an  element  of  a  INxIN  matrix  U' 


ulf  -  n  <•  (A-21) 


n-1 


where  V  denotes  the  polynomial  term  number  in  for  the  equation  at  fth  node  point.  Let 
the  corresponding  function  values  at  each  node  of  the  TV-cube  v't  be  stored  in  the 
column  vector  V,  Furthermore,  let  the  coefficients  a,  of  the  polynomials  for  function 
§  be  stored  in  column  vector  A.  Then  the  substitution  of  IN  nodal  coordinates  u',  and 
corresponding  values  v',  into  equation  (A.  15)  results  in  the  following  set  of  IN 
simultaneous  linear  equations  in  1N  unknowns  coefficients 

{V}  =  \U'\  {A}  (A.22) 

This  equation  is  solved  for  the  polynomial  coefficients  as  follows 

{A}  =  [U"l\  {V}  (A.23) 

Thus  the  value  of  the  function  §  at  the  given  point  {«'}  is  obtained  as  follows 

v*  =  {u}-{A}  (A.24) 

This  completes  the  discussion  of  quadratic  interpolation  in  jV  dimensional  space. 


APPENDIX  B 
SIMULATION  EXAMPLES 


In  this  appendix,  two  examples  of  the  journal  bearing  simulations  are  presented. 
The  first  simulation  is  performed  with  hydrodynamic  model  and  the  second  with 
thermo-hydrodynamic  model.  The  resulting  performance  parameter  during  the 
simulation  are  plotted  for  the  life  of  the  bearing. 

B.l    Simulation  1  —  Hydrodynamic  Model 

For  the  simulation  with  hydrodynamic  model,  the  following  bearing  design  and 

operation  parameters  are  specified. 

B.l.l    Input  Parameters 

Design  Parameters 

Bearing  Diameter,  D  1.32  inch 

Bearing  Speed,  TV  5000  rpm 

Load,  W„  1000  lb 

Bearing  material  Lead-based  Babbitt 

Bearing  material  hardness,  Hb  21 

Journal  material  hardness,  Hs  200 

Contaminant  hardness,  Hc  100 

Bearing  wear  coefficient,  Kb  0.02 

Journal  wear  coefficient,  Ks  0.01 

Rotor  mass,  M  5.0  lbm 

Oil  inlet  temperature,  T,  110  °F 
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Fixed  cost  of  manufacturing,  internal  quality 

and  external  quality  (CMie) 
Tolerance  on  bearing  diameter,  AD 
Bearing  length,  L 
Bearing  clearance,  C, 

Lubricant  oil  viscosity  at  temperature  T„  u 
(Relative  oil  viscosity  (G^) 

Simulation  parameters 

Time  increment  for  simulation 


300.0  $ 
0.0004  in, 
1.75  in, 
0.001  in, 
1.0x10*  Reyn, 
0.15) 


10.0  hour 


Operation  and  Support  Parameters 
Oil  cost  rate,  r0 
Energy  cost  rate,  rt 
Total  oil  volume,  V0 
Initial  contaminants/oil  ratio,  rcoi 

Coefficients  for  evaluating  oil  viscosity  deterioration: 


n0  (equation  3.10) 
t0  (equation  3.11c) 


0.02  $/in3 
0.01  $/kWhr 
250  in3 
l.OxlO"8 


l.OxlO"8 

1.2xl07  revolutions 

50.0  hours 


Design  Constraints 

The  following  constraints  are  used  for  the  determination  of  (i)  the  point  for  oil 
replacements  and,  (b)  the  life  of  the  bearing. 


Lower  limit  on  minimum  film  thickness,  hol 
Upper  limit  on  maximum  oil  pressure,  Pmax 
Maximum  oil  temperature,  T^ 
Lower  limit  on  oil  viscosity,  \xmin 
Upper  limit  on  oil  viscosity,  \xmax 
Upper  limit  on  contaminant/oil  ratio,  rco 
Lower  limit  on  L/D  ratio,  (L/D)min 
Upper  limit  on  L/D  ratio,  (L/D)max 
Lower  limit  on  clearance,  Cmin 
Upper  limit  on  clearance,  Cmax 
Lund  and  Saible  instability  constraint 

B.1.2    Simulation  Results 


5.0x1 0"5  in 
30000  psi 
300  °F 

l.OxlO'7  Reyn 
5.0x1 0"5  Reyn 
l.OxlO"7 
0.25 
1.00 

0.00075  in 
0.02  in 


The  results  of  the  simulations  are  plotted  in  from  Figure  B.l  through  B.16. 
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Figure  B.l:    Oil,  energy,  and  total  cost  versus  hours  of  operation  —  simulation  with 
hydrodynamic  model. 
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Figure  B.2:     Cost  rate  (cr)  versus  hours  of  operation  —  simulation  with  hydrodynamic 
model. 
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Figure  B.3:     Clearance    (C)    versus    hours    of    operation    —    simulation     with 
hydrodynamic  model. 
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Figure  B.4:     Average  oil  temperature  (7^)  versus  hours  of  operation  —  simulation 
with  hydrodynamic  model. 
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Figure  B.5:     Maximum  oil  pressure  (P^  versus  hours  of  operation  —  simulation 
with  hydrodynamic  model. 


Figure  B.6:     Coefficient  of  friction  (/)  versus  hours  of  operation  —  simulation  with 
hydrodynamic  model. 
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Figure  B.7:     Oil    viscosity    (u)    versus    hours    of  operation    —   simulation    with 
hydrodynamic  model. 
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Figure  B.8:     Oil    flow   rate   (Q)   versus    hours   of  operation   —   simulation    with 
hydrodynamic  model. 
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Figure  B.9:     Minimum  oil  film  thickness  (h0)  versus  hours  of  operation  —  simulation 
with  hydrodynamic  model. 
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Figure  B.10:  Contaminants  to  oil  ratio  (rOT)  versus  hours  of  operation  —  simulation 
with  hydrodynamic  model. 
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Figure  B.ll:  Sommerfeld  number  (5)  versus  hours  of  operation  —  simulation  with 
hydrodynamic  model. 
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Figure  B.12:  Modified    Sommerfeld    number    (ST)   versus    hours    of  operation 
simulation  with  hydrodynamic  model. 
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Figure  B.13:  Lund  and  Saible  stability  parameter  (y(o))  versus  hours  of  operation 
simulation  with  hydrodynamic  model. 
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Figure  B.14:  Dimensionless    rotor  mass   (MCu2/W)  versus  hours  of  operation 
simulation  with  hydrodynamic  model. 
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Figure  B.15:  Wear   factor   (WF,)  versus   hours  of  operation  —  simulation   with 
hydrodynamic  model. 
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Figure  B.16:  Wear   factor   (WF2)  versus   hours   of  operation   —  simulation   with 
hydrodynamic  model. 
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B.2   Simulation  2  —  Thermo-hvdrodvnamic  Model 

B.2.1    Input  Parameters 

The  above  example  was  also  run  with  thermo-hydrodynamic  model  with  the 

following  changes. 

Initial  clearance,  C,  0.0015  in 

Relative  viscosity,  G^  0.4 

(Oil  viscosity  u  at  temperature  T,  2.688x10" ) 

B.2.2    Simulation  Results 

The  results  of  the  simulations  are  plotted  in  from  Figure  B.17  through  B.32. 
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Figure  B.17:  Oil,  energy,  and  total  cost  versus  hours  of  operation 
thermo-hydrodynamic  model. 
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Figure  B.18:  Cost  rate  (c,)  versus  hours  of  operation  —  simulation  with  thermo- 
hydrodynamic  model. 
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Figure  B.19:  Clearance  (Q  versus  hours  of  operation  —  simulation  with  thermo- 
hydrodynamic  model. 
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Figure  B.20:  Average  oil  temperature  (7^)  versus  hours  of  operation 
with  thermo-hydrodynamic  model. 
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Figure  B.21:  Maximum  oil  pressure  {P„J  versus  hours  of  operation  —  simulation 
with  hydrodynamic  model. 
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Figure  B.22:  Coefficient  of  friction  (/)  versus  hours  of  operation 
thermo-hydrodynamic  model. 
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Figure  B.23:  Oil  viscosity  (u)  versus  hours  of  operation  —  simulation  with  thermo- 
hydrodynamic  model. 
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Figure  B.24:  Oil  flow  rate  (Q)  versus  hours  of  operation  —  simulation  with  thermo- 
hydrodynamic  model. 
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Figure  B.25:  Minimum  oil  film  thickness  (/»„)  versus  hours  of  operation  —  simulation 
with  thermo-hydrodynamic  model. 
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Figure  B.26:  Contaminants  to  oil  ratio  (rro)  versus  hours  of  operation 
with  thermo-hydrodynamic  model. 


Figure  B.27:  Sommerfeld  number  (5)  versus  hours  of  operation 
thermo-hydrodynamic  model. 
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Figure  B.28:  Modified    Sommerfeld    number    (5*)    versus    hours    of  operation 
simulation  with  thermo-hydrodynamic  model. 
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Figure  B.29:  Lund  and  Saible  stability  parameter  (y(o))  versus  hours  of  operation 
simulation  with  thermo-hydrodynamic  model. 
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Figure  B.30:  Dimensionless    rotor  mass   (MCoT/W)  versus   hours   of  operation 
simulation  with  thermo-hydrodynamic  model. 
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Figure  B.31:  Wear  factor  (WF,)  versus  hours  of  operation  —  simulation  with  thermo- 
hydrodynamic  model. 
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Figure  B.32:   Wear    factor    {WF2)   versus    hours   of  operation    —   simulation    with 
hydrodynamic  model. 
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